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We show that cosmological acceleration (Dark Energy effect) falls into to a special class of quantum 
phenomena that occur on a macroscopic scale. Dark Energy is the third macroscopic quantum effect, 
the first after the discovery of superfluidity and superconductivity. Dark Energy is a phenomenon 
of quantum gravity on the scale of the Universe as a whole at any stage of its evolution, including 
the contemporary Universe. The effect is a direct consequence of the zero rest mass of gravitons, 
conformal non-invariance of the graviton fleld, and one-loop finiteness of quantum gravity, i.e. it 
is a direct consequence of first principles only. Therefore, no hypothetical fields or "new physics" 
are needed to explain the Dark Energy effect. This macroscopic effect of one-loop quantum gravity 
takes place in the empty isotropic non-stationary Universe as well as in such a Universe filled by 
a non-relativistic matter or/and radiation. The effect is due to graviton-ghost condensates arising 
from the interference of quantum coherent states. Each of coherent states is a state of gravitons 
and ghosts of a wavelength of the order of the horizon scale and of different occupation numbers. 
The state vector of the Universe is a coherent superposition of vectors of different occupation num- 
bers. One-loop approximation of quantum gravity is believed to be applicable to the contemporary 
Universe because of its remoteness from the Planck epoch. To substantiate the reliability of macro- 
scopic quantum effects, the formalism of one-loop quantum gravity is discussed in detail. The 
theory is constructed as follows: Faddeev - Popov - De Witt gauged path integral — > factorization 
of classical and quantum variables, allowing the existence of a self-consistent system of equations for 
gravitons, ghosts and macroscopic geometry — > transition to the one-loop approximation, taking 
into account that contributions of ghost fields to observables cannot be eliminated in any way — > 
choice of ghost sector, satisfying the condition of one-loop finiteness of the theory off the mass shell. 
The Bogolyubov-Born-Green-Kirckwood-Yvon (BBGKY) chain for the spectral function of gravi- 
tons renormalized by ghosts is used to build a self-consistent theory of gravitons in the isotropic 
Universe. It is the first use of this technique in quantum gravity calculations. We found three exact 
solutions of the equations, consisting of BBGKY chain and macroscopic Einstein's equations. It was 
found that these solutions describe virtual graviton and ghost condensates as well as condensates of 
instanton fluctuations. All exact solutions, originally found by the BBGKY formalism, are repro- 
duced at the level of exact solutions for fleld operators and state vectors. It was found that exact 
solutions correspond to various condensates with different graviton-ghost compositions. Each exact 
solution corresponds to a certain phase state of graviton-ghost substratum. Quantum-gravity phase 
transitions are introduced. In the formalism of the BBGKY chain, the generalized self-consistent 
theory of gravitons is presented, taking into account the contribution of non-relativistic matter in 
the formation of a common self-consistent gravitational fleld. In the framework of this theory, it 
is shown that the era of non-relativistic matter dominance must be replaced by an era of domi- 
nance of graviton-ghost condensate. Pre-asymptotic state of Dark Energy is a condensate of virtual 
gravitons and ghosts with a constant conformal wavelength. The asymptotic state predicted by the 
theory is a self-polarized graviton-ghost condensate of constant physical wavelength in the De Sitter 
space. The Dark Energy phenomenon of such a nature is presented in the form of AGCDM model 
that interpolates the exact solutions of equations of one-loop quantum gravity. The proposed theory 
is compared with existing observational data on Dark Energy extracted from the Hubble diagram 
for supernovae SNIa. We show that AGCDM model has advantages over ACDM model if criteria 
for the statistical probability are in use. Result of processing of observational data suggests that 
the graviton-ghost condensate is an adequate variable component of Dark Energy. We show that 
its role was significant during the era of large-scale structure formation in the Universe. 

PACS numbers; Dark energy 95.36.+x; Quantum gravity 04.60.-m 
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I. INTRODUCTION 

Macroscopic quantum effects are quantum phenomena that occur on a macroscopic scale. To date, there are two 
known macroscopic quantum effects: superfluidity at the scale of liquid helium vessel and superconductivity at the 
scale of superconducting circuits of electrical current. These effects have been thoroughly studied experimentally and 
theoretically understood. A key role in these effects is played by coherent quantum condensates of micro-objects with 
the De Broglie wavelength of the order of macroscopic size of the system. The third macroscopic quantum effect under 
discussion in this paper is condensation of gravitons and ghosts in the self-consistent field of the expanding Universe. 
Hypotheses on the possibility of graviton condensate formation in the Universe proposed by Hu [Ij and Antoniadis- 
Mazur-Mottola 0] in a general form. A description of these effects by an adequate mathematical formalism is the 
problem at the present time. 

We show that condensation of gravitons and ghosts is a consequence of quantum interference of states forming 
the coherent superposition. In this superposition, quantum fields have a certain wavelength, and with different 
amplitudes of probability they are in states corresponding to different occupation numbers of gravitons and ghosts. 
Intrinsic properties of the theory automatically lead to a characteristic wavelength of gravitons and ghosts in the 
condensate. This wavelength is always of the order of a distance to the horizon of events^. 

In this fact, a common feature of macroscopic quantum effects is manifested: such effects are always formed by 
quantum micro-objects, whose wavelengths are of the order of macroscopic values. With this in mind, we can say 
that macroscopic quantum gravity effects exist across the Universe as a whole. The results of this work suggest 
that the existence of the graviton-ghost condensate is directly responsible for the Dark Energy effect, i.e. for the 
observational data of the acceleration of the expansion of the Universe [1, 01 . Most significantly, a graviton-ghost 
condensate formation is direct consequence of the first principles of the theory of gravity and quantum field theory, 
so that no hypothetical fields are needed to explain the Dark Energy effect. 

Quantum theory of gravity is a non-renormalized theory and for this reason it is impossible to calculate effects 
with an arbitrary accuracy in any order of the theory of perturbations. The program combining gravity with other 
physical interactions within the framework of supergravity or superstrings theory assumes the ultimate formulation 
of the theory containing no divergences. Today we do not have such a theory; nevertheless, we can hope to obtain 
physically meaningful results. Here are the reasons for this assumption. 

First, in all discussed options for the future theory, Einstein's theory of gravity is contained as a low energy limit. 
Second, from all physical fields, which will appear in a future theory (according to present understanding), only the 
quantum component of gravitational field (graviton field) has a unique combination of zero rest mass and conformal 
non-invariance properties. Third, physically meaningful effects of quantum gravity can be identified and quantified 
in one-loop approximation. Fourth, as was been shown by t'Hooft and Veltman [5], the one-loop quantum gravity 
with ghost sector and without fields of matter is finite. For the property of one-loop finiteness, proven in j5:] on the 
graviton mass shell, we add the following key assertion. All one-loop calculations in quantum gravity must be done 
in such a way that the feature of one-loop finiteness (lack of divergences in terms of observahles) must automatically 
he implemented not only on the graviton mass shell but also outside it. 

Let us emphasize the following important fact. Because of conformal non-invariance and zero rest mass of gravitons, 
no conditions exist in the Universe to place gravitons on the mass shell precisely. Therefore, in the absence of one-loop 
finiteness, divergences arise in observables. To eliminate them, the Lagrangian of Einstein's theory must be modified, 
by amending the definition of gravitons. In other words, in the absence of one-loop finiteness, gravitons generate 



Everywhere in this paper we discuss quantum states of gravitons and ghosts that are self-consistent with the evolution of macroscopic 
geometry of the Universe. In the mathematical formalism of the theory, the ghosts play a role of a second physical subsystem, the average 
contributions of which to the macroscopic Einstein equations appear on an equal basis with the average contribution of gravitons. At 
first glance, it may seem that the status of the ghosts as the second subsystem is in a contradiction with the well-known fact that the 
Faddeev— Popov ghosts are not physical particles. However the paradox, is in the fact that we have no contradiction with the standard 
concepts of quantum theory of gauge fields but rather full agreement with these. The Faddeev— Popov ghosts are indeed not physical 
particles in a quantum-field sense, that is, they are not particles that are in the asymptotic states whose energy and momentum are 
connected by a definite relation. Such ghosts are nowhere to be found on the pages of our work. We discuss only virtual gravitons and 
virtual ghosts that exist in the area of interaction. As to virtual ghosts, they cannot be eliminated in principle due to lack of ghost- 
free gauges in quantum gravity. In the strict mathematical sense, the non-stationary Universe as a whole is a region of interaction, 
and, formally speaking, there are no real gravitons and ghosts in it. Approximate representations of real particles, of course, can be 
introduced for shortwave quantum modes. In our work, quantum states of shortwave ghosts are not introduced and consequently are 
not discussed. Furthermore, macroscopic quantum effects, which are discussed in our work, are formed by the most virtual modes 
of all virtual modes. These modes are selected by the equality XH = 1, where A is the wavelength, H is the Hubble function. The 
same equality also characterizes the intensity of interaction of the virtual modes with the classical gravitational field, i.e. it reflects the 
essentially non-perturbative nature of the effects. An approximate transition to real, weakly interacting particles, situated on the mass 
shell is impossible for these modes, in principle (see also the footnote ^ on p. 6). 
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divergences, contrary to their own definition. Such a situation does not make any sense, so the one-loop finiteness off 
the mass shell is a prerequisite for internal consistency of the theory. 

These four conditions provide for the rehabihty of theory predictions. Indeed, the existence of quantum component 
of the gravitational field leaves no doubt. Zero rest mass of this component means no threshold for quantum processes 
of graviton vacuum polarization and graviton creation by external or self-consistent macroscopic gravitational field. 
The combination of zero rest mass and conformal non-invariance of graviton field leads to the fact that these processes 
are occurring even in the isotropic Universe at any stage of its evolution, including the contemporary Universe. Vacuum 
polarization and particle creation belong to effects predicted by the theory already in one-loop approximation. In 
this approximation, calculations of quantum gravitational processes involving gravitons are not accompanied by the 
emergence of divergences. Thus, the one-loop finiteness of quantum gravity allows uniquely describe mathematically 
graviton contributions to the macroscopic observables. Other one-loop effects in the isotropic Universe are suppressed 
either because of conformal invariance of non-gravitational quantum fields, or (in the modern Universe) by non-zero 
rest mass particles, forming effective thresholds for quantum gravitational processes in the macroscopic self-consistent 
field. 

Effects of vacuum polarization and particle creation in the sector of matter fields of J = 0, 1/2, 1 spin were well 
studied in the 1970's by many authors (see [6] and references therein). The theory of classic gravitational waves in 
the isotropic Universe was formulated by Lifshitz in 1946 7]. Grishchuk @] considered a number of cosmological 
applications of this theory that are result of conformal non-invariance of gravitational waves. Isaacson ^ has formu- 
lated the task of self-consistent description of gravitational waves and background geometry. The model of Universe 
consisting of short gravitational waves was described for the first time in [lOl| . The energy-momentum tensor of classic 
gravitational waves of super long wavelengths was constructed in pTl Il2j|. The canonic quantization of gravitational 
field was done in [l^.Tl5|. The local speed of creation of shortwave gravitons was calculated in [l^. In all papers 
listed above, the ghost sector of graviton theory was not taken into account. One-loop quantum gravity in the form 
of the theory of gravitons defined on the background spacetime was described by De Witt [l3| ■ Calculating methods 
of this theory were discussed by Hawking fl^ . For the first time, the approach based on the fact that the ghost sector 
of graviton theory is determined by the condition of one-loop finiteness off the mass shell is presented in the present 
paper. 

One-loop finiteness provides the simplicity and elegance of a mathematical theory that allows, in turn, discovering 
a number of new approximate and exact solutions of its equations. This paper is focused on three exact solutions 
corresponding to three different quantum states of graviton-ghost subsystem in the space of the non-stationary 
isotropic Universe with self-consistent geometry. The first of these solutions describes a coherent condensate of 
virtual gravitons and ghosts; the second solution describes a coherent condensate of instanton fluctuations. The third 
solution describes the self-polarized condensate in the De Sitter space. This solution allows interpretation in terms of 
virtual particles as well as in terms of instanton fluctuations. All three solutions are directly related to the physical 
nature of Dark Energy. 

The principal nature of macroscopic quantum gravity effects, the need for strict proof of their inevitability and re- 
liability impose stringent requirements for constructing a mathematical algorithm of the theory. In our view, existing 
versions of the theory of gravitons in macroscopic spacetime with self-consistent geometry do not meet these require- 
ments. In this connection, note the following fact. Because of conformal non-invariance, the trace of the graviton 
energy-momentum tensor is not zero simply by definition of graviton field. Naturally, the information on macroscopic 
quantum effects (that are the subject of study in this paper) is contained in this trace. After the first publication of 
preliminary results of our research we feel that a number of problems under discussion needs a much more detail 
description. Due to some superficial similarities, effects of graviton and ghost condensation in the De Sitter space are 
perceived sometime as another method of description of conformal anomalies, which are calculated by a traditional 
method of regularization and renormalization. (Calculations of such anomalies see in, e.g., [20l. [2ll. [23. l23l [2^. |25| .) 
We would like to emphasize that such analogies have neither physical nor mathematical basis. Conformal anomalies 
describe the effect of reconstruction of the spectrum of zero oscillations of quantum fields. Their contributions to the 
energy-momentum tensor are independent of state vector of the quantum field. They are parameterized by numerical 
coefficients of the order of unity that are factors of quadratic form of the curvature. Conformal anomalies that are 
microscopic quantum effects are able to contribute to the macroscopic evolution of the Universe only if its parameters 
are close to Planck ones. In contrast, the graviton and ghost condensation is a macroscopic quantum effect in the 
contemporary Universe. It is quantitatively described by macroscopic parameters that govern the structure of a state 
vector. These parameters are averaged numbers of quanta in coherent superpositions. Note also that in all papers 
known to us, calculations of anomalies were conducted in the framework of models with no one-loop finiteness off the 
graviton and ghost mass shell. It is shown in Appendices IXII Al IXII Bl that such models are internally inconsistent. 
The rigorous theory of one-loop quantum gravity presented here shows that the effect of conformal anomalies on 
gravitons must be zero (see Appendix IXII C]) . In the empty Universe where the only gravitons exist, the De Sitter 
space can be formed by only a graviton-ghost condensate. 
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Sections |TT] and IIIII are devoted to the derivation of the equations of the theory with a discussion of all the math- 
ematical details. In Section |TT1 we start with exact quantum theory of gravity, presented in terms of path integral 
of Faddeev-Popov [l^ and De Witt [23] ■ Key ideas of this Section are the following, (i) The necessity to gauge 
the full metric (before its separation into the background and fluctuations) and the inevitability of appearance of a 
ghost sector in the exact path integral and operator Einstein's equations fSections III Al and III Bp : (ii) The principal 
necessity to use normal coordinates (exponential parameterization) in a mathematically rigorous procedure for the 
separation of classical and quantum variables is discussed in Sections III CI and IIIDI (iii) The derivation of differential 
identities, providing the consistency of classical and quantum equations performed jointly in any order of the theory of 
perturbations is given Section [II El Rigorously derived equations of gauged one-loop quantum gravity are presented 
in Section HTFI 

The status of properties of ghost sector generated by gauge is crucial to properly assess the structure of the theory 
and its physical content. Let us immediately emphasize that the standard presentation on the ghost status in the 
theory of S'-matrix can not be exported to the theory of gravitons in the macroscopic spacetime with self-consistent 
geometry. Two internal mathematical properties of the quantum theory of gravity make such export fundamentally 
impossible. First, there are no gauges that completely eliminate the diffeomorphism group degeneracy in the theory 
of gravity. This means that among the objects of the quantum theory of fields inevitably arise ghosts interacting 
with macroscopic gravity. Secondly, gravitons and ghosts cannot be in principle situated precisely on the mass shell 
because of their conformal non-invariance and zero rest mass. This is because there are no asymptotic states, in 
which interaction of quantum fields with macroscopic gravity could be neglected. Restructuring of vacuum graviton 
and ghost modes with a wavelength of the order of the distance to the horizon of events takes place at all stages 
of cosmological evolution, including the contemporary Universe. Ghost trivial vacuum, understood as the quantum 
state with zero occupation numbers for all modes, simply is absent from physically realizable states. Therefore, direct 
participation of ghosts in the formation of macroscopic observables is inevitable^. 

Section IIIII is devoted to general discussion of equations of the theory of gravitons in the isotropic Universe. It 
focuses on three issues: (i) The gauge invariant procedure for eliminating gauge non-invariant modes by conditions 
imposed on the state vector (Section IIII A^ : (ii) Construction of the state vector of a general form as a product 
of normalized superpositions (Section IIII C|) ; (iii) The allocation of class of legitimate gauges that are invariant with 
respect to transformations of the symmetry group of the background spacetime while providing the one-loop finiteness 
of macroscopic observables fSections IIIIDI and IIIIE]) . The main conclusion is that the quantum ghost fields are 
inevitable and unavoidable components of the quantum gravitational field. As noted above, one-loop finiteness is seen 
by us as a universal property of quantum gravity, which extends off the mass shell. The requirement of compensation 
of divergences in terms of macroscopic observables, resulting from one-loop finiteness, uniquely captures the dynamic 
properties of quantum ghost fields in the isotropic Universe. 

The treatment set out in Sections |TT] and IIIII in essence, is a sequence of transformations of equations defined by 
original gauged path integral. Basically, these transformations are of mathematical identity nature. There are only 
three elements of the theory, missing in the original integral: 

(i) The hypothesis of the existence of classical spacetime with the deterministic but self- consistent geometry; 

(ii) A transition to the one-loop approximation in the self-consistent system of classical and quantum equations; 

(iii) A class of gauges that automatically provides the one-loop finiteness of self- consistent theory of gravitons in 



^ Once again, we emphasize that the equal participation of virtual gravitons and ghosts in the formation of macroscopic observables in 
the non-stationary Universe does not contradict the generally accepted concepts of the quantum theory of gauge fields. On the contrary 
it follows directly from the mathematical structure of this theory. In order to clear up this issue once and for all, recall some details 
of the theory of S— matrix. In constructing this theory, all space— time is divided into regions of asymptotic states and the region of 
effective interaction. Note that this decomposition is carried out by means of, generally speaking, an artificial procedure of turning on 
and off the interaction adiabatically. (For obvious reasons, the problem of self— consistent description of gravitons and ghosts in the 
non-stationary Universe with \H = 1 by means of an analogue of such procedure cannot be considered a priori.) Then, after splitting 
the space-time into two regions, it is assumed that the asymptotic states are ghost-free. In the most elegant way, this selection rule is 
implemented in the BRST formalism, which shows that the BRST invariant states turn out to be gauge-invariant automatically. The 
virtual ghosts, however, remain in the area of interaction, and this points to the fact that virtual gravitons and ghosts are parts of 
the Feynman diagrams on an equal footing. In the self-consistent theory of gravitons in the non— stationary Universe, virtual ghosts of 
equal weight as the gravitons, appear at the same place where they appear in the theory of 5— matrix, i.e. at the same place as they 
were introduced by Feynman, i.e. in the region of interaction. Of course, the fact that in the real non-stationary Universe, both the 
observer and virtual particles with XH = 1 are in the area of interaction, is highly nontrivial. It is quite possible that this property of 
the real world is manifested in the effect of dark energy. An active and irremovable participation of virtual ghosts in the formation of 
macroscopic properties of the real Universe poses the question of their physical nature. Today, we can only say with certainty that the 
mathematical inevitability of ghosts provides the one-loop finiteness off the mass shell, i.e. the mathematical consistency of one-loop 
quantum gravity without fields of matter. Some hypothetical ideas about the nature of the ghosts are briefly discussed in the final 
Section |X] 
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the isotropic Universe. 

Conceptually discussing of the first two additional elements was not necessary. Their introduction to the formalism 
is to the factorization of measure of path integral and expansion of equations of the theory in a series of powers of the 
graviton field. The existence of appropriate correct mathematical procedures is not in doubt. Agreement on choosing 
of a gauge is also mathematically consistent. Moreover, a gauge is necessary for the strict definition of path integral 
as a mathematical object. 

The existence of class of gauges, automatically providing one-loop finiteness off the mass shell is itself a nontrivial 
property of the theory. The assertion that only such gauges can be used in the self-consistent theory of gravitons in 
the isotropic Universe is actually a condition for the internal consistency of the theory. In Appendix IXII Bl we present 
the formal proof that an alternative formulation of the theory (with no one-loop finiteness) does not exist. 

From the requirement of the one-loop finiteness it follows that the quantum component of the gravitational field 
is of a heterogeneous graviton-ghost structure. In further Sections of work, it appears that this new element of the 
theory prejudge its physical content. 

Sections IIV Al and IIV Bl contain approximate solutions to obtain quantum ensembles of short and long gravitational 
waves. In Section IIV CI it is shown that approximate solutions obtained can be used to construct scenarios for the 
evolution of the early Universe. In one such scenario, the Universe is filled with ultra-relativistic gas of short-wave 
gravitons and with a condensate of super-long wavelengths, which is dominated by ghosts. The evolution of this 
Universe is oscillating in nature. 

At the heart of cosmological applications of one-loop quantum gravity is the Bogolyubov-Born-Green-Kirckwood- 
Yvon (BBGKY) chain (or hierarchy) for the spectral function of gravitons, renormalized by ghosts. We present the first 
use of this technique in quantum gravity calculations. Each equation of the BBGKY chain connects the expressions for 
neighboring moments of the spectral function. In Section [VAI the BBGKY chain is derived by identical mathematical 
procedures from graviton and ghost operator equations. Among these procedures is averaging of bilinear forms of 
field operators over the state vector of the general form, whose mathematical structure is given in Section IlII CI The 
need to work with state vectors of the general form is dictated by the instability of the trivial graviton-ghost vacuum 
fSection IIIIF[) . Evaluation of mathematical correctness of procedures for BBGKY structure is entirely a question 
of the existence of moments of the spectral function as mathematical objects. A positive answer to this question is 
guaranteed by one-loop finiteness fSection llllPp . The set of moments of the spectral function contains information 
on the dynamics of operators as well as on the properties of the quantum state over which the averaging is done. The 
set of solutions of BBGKY chain contains all possible self-consistent solutions of operator equation, averaged over all 
possible quantum ensembles. 

A nontrivial fact is that in the one-loop quantum gravity BBGKY chain can formally be introduced at an axiomatic 
level. Theory of gravitons provided by BBGKY chain, conceptually and mathematically corresponds to the axiomatic 
quantum field theory in the Wightman formulation (see Chapter 8 in the monograph [28]). Here, as in Wightman, the 
full information on the quantum field is contained in an infinite sequence of averaged correlation functions. Definitions 
of these functions clearly relate to the symmetry properties of manifold on one this field is defined. O77.ce the BBGKY 
chain is set up, the existence of finite solutions for the observables is provided by inherent mathematical properties of 
equations of the chain. This means that the phenomenology of BBGKY chain is more general than field operators, 
state vectors and graviton-ghost compensation of divergences that were used in its derivation. 

Exact solutions of the equations, consisting of BBGKY chain and macroscopic Einstein's equations are obtained in 
Sections IV Bl and IV CI Two solutions given in lVBI describe heterogeneous graviton-ghost condensates, consisting of 
three subsystems. Two of these are condensates of spatially homogeneous modes with the equations of state p = — e/3 
and p = e. The third subsystem is a condensate of quasi-resonant modes with a constant conformal wavelength 
corresponding to the variable physical wavelength of the order of the distance to the horizon of events. The equations 
of state of condensates of quasi-resonant modes differ from p — e/3 by logarithmic terms, through which the first 
solution is p > — £/3, while the second is p ^ —e/3. Furthermore, the solutions differ by the sign of the energy 
density of condensates of spatially homogenous modes. The third solution describes a homogeneous condensate of 
quasi-resonant modes with a constant physical wavelength. The equation of state of this condensate is p = —e and 
its self-consistent geometry is the De Sitter space. The three solutions are interpreted as three different phase states 
of graviton-ghost system. The problem of quantum-gravity phase transitions is discussed in Sectior fVDI 

Solutions obtained in Section|V]in terms of moments of the spectral function, are reproduced in Sections I VII and I VIII 
at the level of dynamics of operators and state vectors. A microscopic theory provides details to clarify the structure 
of graviton-ghost condensates and clearly demonstrates the effects of quantum interference of coherent states. In 
Section IVl Al it is shown that the condensate of quasi-resonant modes with the equation of state p > —e/3 consists of 
virtual gravitons and ghosts. In Section [VI Bl a similar interpretation is proposed for the condensate in the De Sitter 
space, but it became necessary to extend the mathematical definition of the moments of the spectral function. 

New properties of the theory, whose existence was not anticipated in advance, are studied in Section IVIII In 
Section IVII Al wc find that the self-consistent theory of gravitons and ghosts is invariant with respect to the Wick 
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turn. In this Section we also construct the formalism of quantum theory in the imaginary time and discuss the physical 
interpretation of this theory. The subjects of the study are correlated fluctuations arising in the process of tunnelling 
between degenerate states of graviton-ghost systems, divided by classically impenetrable barriers. The level of these 
fluctuations is evaluated by instanton solutions (as in Quantum Chromodynamics). In Section [VII Bl it is shown that 
the condensate of quasi-resonant modes with the equation of state p < — e/3 is of purely instanton nature. In Section 
IVIICI thc instanton condensate theory is formulated for the De Sitter space. 

Potential use of the results obtained to construct scenarios of cosmological evolution was briefly discussed in Sections 
IIVI — I VIII to obtain approximate and exact solutions. A main application of the theory of macroscopic effects of 
quantum gravity is to explain the physics of Dark Energy. As a carrier of Dark Energy, quantum gravity offers 
graviton-ghost condensate of quasi-resonant modes. In Section IVIIII the theory of gravitons and ghosts is formulated 
in a self-consistent field, in the formation of which heavy particles (baryons and particles of Dark Matter (neutralino 
as a example)) are involved. In Section lVlII Al it is shown that the cosmological solution describing the evolution of the 
scale factor, non-relativistic matter and graviton-ghost subsystem, automatically satisfies the chain integral identities 
contained in the equations of this theory. From these identities it follows that the stage of the Universe evolution, which 
is dominated by non-relativistic matter, must inevitably be replaced by a stage at which substantive contribution to 
the energy density of the cosmological substratum will produce graviton-ghost condensate of quasi-resonant modes 
of a constant conformal wavelength. This stage of the Universe evolution we consider as the pre-asymptotic one. 
Conversion of pre-asymptotic condensate to a condensate of constant physical wavelength occurs in the process of 
quantum-gravity phase transition. 

In Section IVIIIBI the exact solutions of equations of the one-loop quantum gravity are used to construct model of 
Dark Energy, intended for interpretation of experimental data. The model is based on simple interpolation formulae, 
describing the energy density and pressure of graviton-ghost condensates at various stages of cosmological evolution. 
For the cosmological model based on the physical considerations presented here, we suggest the abbreviation " AGCDM 
model." Here "CDM" means Cold Dark Matter; "G" means a graviton-ghost condensate in the role of Dark Energy; 
symbol A indicates an asymptotic state of the Universe in which the energy density of cosmological substratum goes 
to a constant value. The following cosmological Einstein equations correspond to the AGCDM model 

3H^^K[eg + ^), 6Q = -x[sg + 3pg + ^) , (I.l) 



where 



,, = A. + ai.^i. p. = -A.-ii,„5£. ,1.2) 

are energy density and pressure of graviton-ghost medium; M/a^ is the density of non-relativistic matter; H = d/ais 
Hubble function; Q = a/a is acceleration of the Universe expansion. The symbol Aqo indicates the asymptotic value of 
total energy density of graviton-ghost condensate and equilibrium vacuum of non-gravitational physical fields. Note 
that Aoo ^ 0, even if the non-gravitational vacuum energy ("standard" Einstein's cosmological constant) vanishes by 
some compensation mechanisms of the type of supersymmetry. 

As can be seen from (jl.2[) . the theory of graviton-ghost condensates in an extremely simplified version predicted 
by a three-parameter model of Dark Energy. A quantitative comparison of theory with the observational data is 
conducted in Section IIXI Section IIX Al devoted to a brief review of the physics of Dark Energy problems. Synopsis 
of observational data on supernovae SNIa is given in Section flX Bl Our focus is on the choice between ACDM and 
AGCDM models based on general principles of fundamental physics without involving hypothetical elements. In 
Section IIX CI a comparative analysis of the results of processing the Hubble diagram for SNIa by formulas for both 
models is conducted. Here we show that the AGCDM model not only explains quantitatively the effect of Dark 
Energy, but also reproduces a number of specific details of observational data. 

The results and problems of the theory are briefiy discussed in the Conclusion (Section|X|. Appendix IXII is devoted 
to the cosmological constant problem within the framework of its interpretation as the energy density of equilibrium 
vacuum of non-gravitational physical fields. An internal inconsistency of the theory which lacks one-loop finiteness off 
the mass shell is proved in Appendices IXII AIXII Bl In IXII C| it is shown by the method of dimensional transmutation 
that the one-loop self-consistent theory of gravitons and ghosts is not only finite but is also free of anomalies. 

A system of units is used, in which the speed of light is c = 1, Planck constant is h = 197.327 MeV-fm; Einstein's 
gravity constant is >c = SnG = 8n ■ 1.324 • 10-*^ MeV'^-fm. 



II. BASIC EQUATIONS 

According to De Witt [l^ , one of formulations of one-loop quantum gravity (with no fields of matter) is reduced 
to the zero rest mass quantum field theory with spin J = 2, defined for the background spacetime with classic metric. 
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The graviton dynamics is defined by the interaction between quantum field and classic gravity, and the background 
space geometry, in turn, is formed by the energy-momentum tensor (EMT) of gravitons. 

In the current Section we describe how to get the self-consistent system of equations, consisting of quantum operator 
equations for gravitons and ghosts and classic C-number Einstein equations for macroscopic metrics with averaged 
EMT of gravitons and ghosts on the right hand side. The theory is formulated without any constrains on the graviton 
wavelength that allows the use of the theory for the description of quantum gravity effects at the long wavelength 
region of the specter. The equations of the theory (except the gauge condition) are represented in 4D form which is 
general covariant with respect to the transformation of the macroscopic metric. 

The mathematically consistent system of 4D quantum and classic equations with no restrictions with respect to 
graviton wavelengths is obtained by a regular method for the first time. The case of a gauged path integral with ghost 
sector is seen as a source object of the theory. Important elements of the method are exponential parameterization of 
the operator of the density of the contravariant metric; factorization of path integral measure; consequent integration 
over quantum and classic components of the gravitational field. Mutual compliance of quantum and classic equations, 
expressed in terms of fulfilling of the conservation of averaged EMT at the operator equations of motion is provided 
by the virtue of the theory construction method. 

A. Path Integral and Faddeev— Popov Ghosts 

Formally, the exact scheme of quantum gravity is based on the amplitude of transition, represented by path integral 
MM: 



(II.l) 



(out|in) = J exp Q JiLgrav + LA)d*x^ (detM\) (^SiAk^grf" - B'}j dfi 



J exp (^^ J {Lgrav + La + Lghost)d^x^ n [Ji^i^kV^g''' - B'U dfidfie , 



where 



Lgrav + -7;—V~^g''''Rik " a/-^A 



is the density of gravitational Lagrangian, with cosmological constant included; Lghost is the density of ghost La- 
grangian, explicit form of which is defined by localization of detM'^,; Ak is gauge operator, B^{x) is the given field; 
M\ is an operator of equation for infinitesimal parameters of transformations for the residual degeneracy 77' = (5a;* ; 

X [ i<^k J 

is the gauge invariant measure of path integration over gravitational variables; dfie is the measure of integration over 
ghost variables. Operator M\ is of standard definition: 



M\r^' ^ ^fc(V-5r) = 0, (II.3) 
where 

Sy/~gg''' = -dii^g'W) + V^g^'dir^" + ^gg^'dir,^ (n.4) 

is variation of metrics under the action of infinitesimal transformations of the group of diffeomorphisms. According 
to (jll.ip . the allowed gauges are constrained by the condition of existence of the inverse operator {M\)~^. 

The equation pi.ip explicitly manifests the fact that the source path integral is defined as a mathematical object 
only after the gauge has been imposed. In the theory of gravity, there are no gauges completely eliminating the 
degeneracy with respect to the transformations (|II.4p . Therefore, the sector of nontrivial ghost fields, interacting with 
gravity, is necessarily present in the path integral. This aspect of the quantum gravity is important for understanding 
of its mathematical structure, which is fixed before any approximations are introduced. By that reason, in this Section 
we discuss the equations of the theory, by explicitly defining the concrete gauge. 
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In cosmological applications of the quantum gravity it is convenient to use synchronous gauges of type: 

V^9°° = Vl : V^5°" = . (II.5) 

For that gauge 

Ak = (1, 0, 0, 0) , B' = (V^, 0, 0, 0), (II.6) 

where 7 = 7(x) is the metric determinant of the basic 3D space of constant curvature (for the plane cosmological 
model 7=1). More general approach to the choice of gauge used for cosmological problems is discussed in Section 

nni 

The construction of the ghost sector, i.e. finding of the Lagrangian density Lghost, is reduced to two operations. 
First, det M*j, is represented in the form, factorized over independent degrees of freedom for ghosts, and then the 
localization of the obtained expression is conducted. Substitution of (|II.6p and pi.4p to pi.3p gives the following 
system of equations 

+ V^^ = o, yz-g"^a,^o ^ ^ ^^^^^ 

According to pi.7p . with respect to variables 77", the operator-matrix M*^, reads 

d 



(II., 



(Note matrix-operator is obtained in the form pi.Sp without the substitution of transformation parameters if 
Leutwiller measure dfiL = gg^^djl is used. The measure discussion see, e.g. 0].) Functional determinant of matrix- 
operator det M*^ is represented in the form of the determinant of matrix M*^,, every element of which is a functional 
determinant of differential operator. As it is follows from pl.8|) . 

detAf\. - (det d.^gg'^'dk) x fdet -|) x fdct |-) . (II.9) 



dtj \ at. 

One can see that the first multiplier in pi.9|) is 4-invariant determinant of the operator of the zero rest mass Klein- 
Gordon-Fock equation, and two other multipliers do not depend on gravitational variables. 

Localization of determinant pi.9|) by representing it in a form of path integral over the ghost fields is a trivial 
operation. As it follows from pi.9p . the class of synchronous gauges contains three dynamically independent ghost 
fields 9, iy5, X, two of each — (p, x do not interact with gravity. For the obvious reason, the trivial ghosts ip, x 
excluded from the theory. The Lagrangian density of nontrivial ghosts coincides exactly with Lagrangian density of 
complex Klein-Gordon-Fock fields (taking into account the Grassman character of fields 0, 9): 

Lghost = -^^^§'"8,9 ■ dk9. (11.10) 

The normalization multiplier —l/ix in pi.lO[) is chosen for the convenience. The integral measure over ghost fields 
has a simple form: 

X 

The calculations above comply with both general requirements to the construction of ghost sector. First, path 
integration should be carried out only over the dynamically independent ghost fields. Second, in the ghost sector, it is 
necessary to extract and then to take into account only the nontrivial ghost fields, i.e. those interacting with gravity. 

The extraction of dynamically independent nontrivial ghost fields can be done not only by factorization of functional 
determinant (as it made in pi.9|) ). but by means of researching the equations for the ghosts as well. It is well known 
[30I I , that from the definition of det M\ it follows that the ghost equations coincide with equations for parameters of 
infinitesimal transformations of residual degeneracy. Therefore, according to pl.sp . we can immediately get M\9^ — 0, 
where 9^ are the Grassman fields. For the gauge pi.Sp with i? = 1 we get 

-d^^9^ + ^— = Q, (ILU) 
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V^99"%0' + ^^^ = O. (11.12) 

From (|II.12p it follows that those transversal components of vector 9" yield the equations dO'^/dt — 0, i.e. these 
degrees of freedom correspond to two trivial ghosts non-interacting with gravity. Note, that in pi. lip only the 
longitudinal component 9"^ is present. Equation (jII.12p has a status of equation connecting the longitudinal component 

^jj* and function 9^ . It means that the longitudinal field 0| is not dynamically independent. Equation for dynamically 
independent degrees of freedom 9^ is obtained as follows. First, the operator is applied to the equation pi.l2p . next 
the result is substituted into time-differentiated equation pi. lip . After the substitution, one gets Klein-Gordon-Fock 
equation, 

QlVlg-^ + d^^gr^d09 = = , (11.13) 

where 9^9°. Reconstruction of the ghost Lagrangian pi.lOp from dynamical equation pi.lSp does not require an 
additional explanation. 

Thus, the prove of wave properties of ghosts in the class of synchronous gauges is done by two methods with clear 
correlations between objects and operations used in these methods. 

We will return to the discussion of gauges and ghosts in the Section IlII El 



B. Einstein Operator Equations 



Let us take into account the fact that the calculation of gauged path integral should be mathematically equivalent to 
the solution of dynamical operator equations in the Heisenberg representation. It is also clear that operator equations 
of quantum theory should have a definite relationship with Einstein equations. In the classic theory, it is possible to 
use any form of representation of Einstein equations, e.g. 

(-5)" if Rim - U'^'g^'^Rirr. - ff^'^^A^) = , (a) 



(-5)" - - '5Na) = , (&) (11.14) 

(-5)" {R^k - \mkg''''Rim - 5''>^a) = , (c) 

where, for example, n = 0, 1/2, 1. Transition from one to another is reduced to the multiplication by metric tensor 
and its determinant, which are trivial operations in case when the metric is a C-number function. If the metric is an 
operator, then the analogous operations will, at least, change renormalization procedures of quantum non-polynomial 
theory. Thus, the question about the form of notation for Einstein's operator equations has first-hand relation to 
the calculation procedure. Now we show that in the quantum theory one should use operator equations pi. 14b ) with 
n = 1/2, supplemented by the energy-momentum pseudo-tensor of ghosts. 

In the path integral formalism, the renormalization procedures are defined by the dependence of Lagrangian of 
interactions and the measure of integration of the field operator in terms of which the polynomial expansion of non- 
polynomial theory is defined js^]. The introduction of such an operator, i.e. the parameterization of the metric, is, 
generally speaking, not simple. Nevertheless, it is possible to find a special parameterization for which the algorithms 
of renormalization procedures are defined only by Lagrangian of interactions. Obviously, in such a parameterization 
the measure of integration should be trivial. It reads: 

(IL15) 

where 'l'*^ is a dynamic variable. The metric is expressed via this variable. It is shown in [soj that the trivialization 
of measure pi.lSI) takes place for the exponential parameterization that reads 

V^a'" = V=^r'(exp *)f = V^f (^<5f + + i*r^^^ + . . .) , (11.16) 

where g**^ is the defined metric of an auxiliary basic space. In that class of our interest, the metric is defined by the 
interval 
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where ^ap is the metric of 3D space with a constant curvature. (For the flat Universe is the EucHd metric.) 

The exponential parameterization is singled out among all other parameterizations by the property that ^'f are the 
normal coordinates of gravitational fields [3i|. In that respect, the gauge conditions (jll.SP are identical to 0. The 

fact that the "gauged" coordinates arc the normal coordinates, leads to a simple and elegant ghost sector (|II.10[) . The 
status of ^j?*^, as normal coordinates, is of principal value for the mathematical correctness while separating the classic 
and quantum variables (see Section IlIPp . Besides, in the framework of perturbation theory the normal coordinates 
allow to organize a calculation procedure, which is based on a simple classification of nonlinearity of quantum gravity 
field. It is important that this procedure is mathematically non-contradictive at every order of perturbation theory 
over amplitude of quantum fields (see Section III El IIIFp . 

Operator Einstein equations that are mathematically equivalent to the path integral of a trivial measure are derived 
by the variation of gauged action by variables '^^ . The principal point is that the gauged action necessarily includes 
the ghost sector because there are no gauges that are able to completely eliminate the degeneracy. According to 
(|II.10p . in the class of synchronous gauges we get 



s^ - j d''x\ ^V-ar [R^k + i^d,6 ■dke]+ V-5A \ . (11.17) 



In accordance with definition pi.l6|) . the variation is done by the rule 
Thus, from ([IlTtI) it follows 



^ ^gg'^'Ra - x ^V~99'%"""'^ V^gr'fj^f"^ V^S^Aj = . (11.18) 

After subtraction of semi-contraction from pi.l8[) we obtain a mathematically equivalent equation 



(11.19) 

-gg'^^Ru - \s^V^gr'Rn.i - - {V^gg'%"""'^ + V^gS^A) = . 
In (|lLT8l) . pLTQ]) there is an object 

f^t"''^ = - ^ {d,e ■ dke + dj ■ d,e ~ m^g'^^diQ ■ d„^e) , (11.20) 

which has the status of the energy-momentum pseudo-tensor of ghosts. 

In accordance with the general properties of Einstein's theory, six spatial components of equations (jll.lSp are 
considered as quantum equations of motion: 

^gf'R^i = X (v^/'f (f° - l^f yZ-^-^fi^^''-*) _ yZ-^f a) . (11.21) 

(Everywhere in this work the Greek metric indexes stand for a, /3 = 1, 2, 3.) In the classic theory, equations 
of constraints £g = and £q = are the first integrals of equations of motion (|II.2ip . Therefore, in the quantum 
theory formulated in the Heisenberg representation four primary constraints from (jll.lQp , have the status of the initial 
conditions for the Heisenberg state vector. They read: 



-99''Roi ~ IV^r'Rrni - K + V^gA) ^ |vl.) = 



2 



(11.22) 







If conditions (|II.22p are valid from the start, then the internal properties of the theory must provide their validity 
at any subsequent moment of time. Four secondary relations, defined by the gauge non containing the higher order 
derivatives, also have the same status: 



{^(v^5^')-i?'}l*)=0 



(11.23) 
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The system of equations of quantum gravity is closed by the ghosts' equations of motion, obtained by the variation 
of action pi.l7[) over ghost variables: 



(11.24) 

Ghost fields 9 and 6 are not defined by Grassman scalars, therefore xj^sh-ost) ^ tensor. Nevertheless, all 

mathematical properties of equations (jll.24p and expressions (jll.20p coincide with the respected properties of equations 
and EMT of complex scalar fields. This fact is of great importance when concrete calculations are done (see Section 

nn). 



C. Factorization of the Path Integral 

Transition from the formally exact scheme (|II.2ip — pi.24p ) to the semi-quantum theory of gravity can be done 
after some additional hypotheses are included in the theory. The physical content of these hypotheses consists of the 
assertion of existence of classical spacetime with metric gik , connectivity F^; and curvature Rik ■ The first hypothesis 
is formulated at the level of operators. Assume that operator of metric g^^ is a functional of C-numher function 
g"^^ and the quantum operator ip^. The second hypothesis is related to the state vector. Each state vector that is 
involved in the scalar product (outjin), is represented in a factorized form \'^) = jC')!'!/;), where \4>) are the vectors of 
quantum states of gravitons; |<i>) are the vectors of quasi-classic states of macroscopic metric. In the framework of 
these hypotheses the transitional amplitude is reduced to the product of amplitudes: 

(out|in) = {^out\^tn){i'out\lpzn) ■ (11.25) 

Thus, the physical assumption about existence of classic spacetime formally (mathematically) means that the path 
integral must be calculated first by exact integration over quantum variables, and then by approximate integration over 
the classic metric. 

Mathematical definition of classic and quantum variables with subsequent integrations are possible only after the 
trivialization and factorization of integral measure are done. As already noted, trivial measure (|II.15|) takes place in 
exponential parameterization pi.l6p . The existence of |in) = \'^) vector allows the introduction of classic C-number 
variables as follows 

a>f = (vfl^flvl/) , V^g^'^ = V^/(exp<i>)f . 

Quantum graviton operators are defined as the difference V^f = ^\ — Factorized amplitude pi.25p is calculated 
via the factorized measure 

dpL = d^Lg X d^^ , 

f 1 (11.26) 

- n S ^-af" n dg'' \ , dt,^ = \{\{ d^l . 

X I i^k J X i^k 

Factorization of the measure allows the subsequent integration, first by d/i^, d/xg, then by approximate integration 
over dfig . In the operator formalism, such consecutive integrations correspond to the solution of self-consistent system 
of classic and quantum equations. Classical equations are obtained by averaging of operator equations pi.l9p . They 
read: 

m£^\^)^0. {11.27) 
Subtraction of pi.27p from pi.l9p gives the quantum dynamic equations 

4'=-(vI/|ff|*)=0. (n.28) 
Synchronous gauge pi.23p is converted to the gauge of classical metric and to conditions imposed on the state vector: 



(11.29) 



14 



Quantum equations (jll.24p of ghosts' dynamics are added to equations (jll.27p — (jll.29[) . 

Theory of gravitons in the macroscopic spacetime with self-consistent geometry is without doubt an approximate 
theory. Formally, the approximation is in the fact that the single mathematical object \/—gg^^ is replaced by two 
objects — classical metric and quantum field, having essentially different physical interpretations. That "coercion" of 
the theory can lead to a controversy, i.e. to the system of equations having no solutions, if an inaccurate mathematics 
of the adopted hypotheses is used. The scheme described above does not have such a controversy. The most important 
element of the scheme is the exponential parameterization pi.l6[) . which separates the classical and quantum variables, 
as can be seen from (jll.26p . After the background and quantum fluctuations are introduced, this parameterization 
looks as follows: 

= V^f (exp (<!> + ^)) ' = V^g'' (exp , (11.30) 

Note that the auxiliary basic space vanishes from the theory, and instead the macroscopic (physical) spacetime with 
self-consistent geometry takes its place. 

If the geometry of macroscopic spacetime satisfies symmetry constrains, the factorization of the measure (III.26|) 
becomes not a formal procedure but strictly mathematical in its nature. These restrictions must ensure the existence 
of an algorithm solving the equations of constraints in the framework of the perturbation theory (over the amplitude 
of quantum fields) . The theory of gravity is non-polynomial, so after the separation of single field into classical and 
quantum components, the use of the perturbation theory in the quantum sector becomes unavoidable. The classical 
sector remains non-perturbative. In the general case, when quantum field is defined in an arbitrary Riemann space, 
the equations of constraints is not explicitly solvable. The problem can be solved in the framework of perturbation 
theory if background gik and the free (linear) tensor field tp'^ belong to different irreducible representations of the 
symmetry group of the background spacetime. In that case at the level of linear field we obtain (|II.26|) . because the 
full measure is represented as a product of measure of integration over independent irreducible representations. At 
the next order, factorization is done over coordinates, because the classical background and the induced quantum 
fluctuations have essentially different spacetime dynamics. Note, to factorize the measure by symmetry criterion we 
do not need to go to the short-wave approximation. 

Background metric of isotropic cosmological models and classical spherically symmetric non-stationary gravitational 
field meet the constrains described above. These two cases are covering all important applications of semi-quantum 
theory of gravity which are quantum effects of vacuum polarization and creation of gravitons in the non-stationary 
Universe and in the neighborhood of black holes. 



D. Variational Principle for Classic and Quantum Equations 

Geometrical variables can be identically transformed to the form of functionals of classical and quantum variables. 
At the first step of transformation there is no need to fix the parameterization. Let us introduce the notations: 

v-5 V-5 (11.31) 

According to (|II.30p , formalism of the theory allows definition of quantum field ip^ as symmetric tensor in physical 
space, gkiWk — i^ik = i'ki- Objects, introduced in pi.3ip . have the same status. With any parameterization the 
following relationships take place: 

lim X''' lim Y,k ^ gzk ■ 

We should also remember that the mixed components of tensors , do not contain the background metric as 
functional parameters. For any parameterization, these tensors are only functionals of quantum fields ipf which are 
also defined in mixed indexes. For the exponential parameterization: 

= sf + + l^i^f + ... , = - v3f + l^pl^pt + ... , 

^ ^ (11.32) 



g = g-d = ge 
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where d = det One can seen from (|II.32[) . that the determinant of the full metric contains only the trace of the 

quantum field. 

Regardless of parameterization, the connectivity and curvature of the macroscopic space T\f., R\ijy^ are extracted 
from full connectivity and curvature as additive terms: 

ik ~ ^ ik I -'ifc' ^klm — ^klm + '"^ klm ■ 

Quantum contribution to the curvature tensor, 

klm km -,1 kl ;m ' nl km nm kl ' 

is expressed via the quantum contribution to the full connectivity: 

= I {-y^mX% - YkmX% + %Y,^X"^^X^-„,) + ir,„ (^SlX^-, + 5{X^-, - Y,kX"^'X^-^) . (11.33) 
The density of Ricci tensor in mixed indexes reads 



§g'''Ru = V^lx'''Ru + l: 



2 




V V V \ vkl yjm yns i^v , V"' 



(11.34) 



Symbol " ;" in pi.33|) . pi.34|) and in what follows stands for the covariant derivatives in background space. The density 
of gauged gravitational Lagrangian is represented in a form which is characteristic for the theory of quantum fields 
in the classical background spacetime: 



S — f d'^xJ—q ( Corav — VrfA — 



-g[Cgrav-VdA-—X'^0,, 



^nrav — X Rik ^" ^ 

2k 8m: 



vkl I V V V V 1 vi"'- 0V-, V*' vkm 

^ I J-jnJ^sm n J- jm-i sn I ^ .j,^ ;( ^J-zk^ -m^ d 



(11.35) 



When the expression for Cgrav was obtained from contraction of tensor (|II.34p . the full covariant divergence in the 
background space have been excluded. Formulas (jll.34[l . (jll.35p apply for at any parameterization. 

Let us discuss the variation method. In the exact quantum theory of gravity with the trivial measure (jII.15[) . the 
variation of the action over variables leads to the Einstein equations in mixed indexes (jII.18[) and pL19p . In the 
exact theory, the exponential parameterization is convenient, but, generally speaking, is not necessary. A principally 
different situation takes place in the approximate self-consistent theory of gravitons in the macroscopic spacetime. In 
that theory the number of variables doubles, and with this, the classical and quantum components of gravitational 
fields have to have the status of the dynamically independent variables due to the doubling of the number of equations. 
The variation should be done separately over each type of variables. The formalism of the path integration suggests a 
rigid criterion of dynamic independence: the full measure of integration, by definition, must be factorized with respect 
to the dynamically independent variables. Obviously, only the exponential parameterization pi.30p . leading to the 
factorized measure pi.26l) . meets the criterion. 

The variation of the action over the classic variables is done together with the operation of averaging over the 
quantum ensemble. In the result, equations for metric of the macroscopic spacetime are obtained: 

(*l^l*> = -2MV^gnk{^\G'l - lsfG\\^) = , (11.36) 

where — Gi / \/—g- Variation of the action over background variables, defined as "I>^' = (^'j^'^'I^I'), yields the 
equations: 

(*l^l*> = -2>.^/^(*|G,f I*) = . (11.37) 

Equations (III.36|) and pl.37p are mathematically identical. We should also mention that if the variations over the 
background metric are done with the fixed mixed components of the quantum field, these equations are valid for any 
parameterization. 
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Exponential parameterization (III.30p has a unique property: the variations over classic (before averaging) and 
quantum il)^ (without averaging) variables lead to the same equations. That fact is a direct consequence of the relations, 
showing that variations 5$f and Sipf are multiplied by the same operator multiplier: 



5\ 








= const , 




^gr = ^ 


/~gg''s^^f , 




= const . 



By a simple operation of subtraction, the identity allows the extraction of pure background terms from the equation 
of quantum field. The equations of graviton theory in the macroscopic space with self-consistent geometry are written 
as follows: 

(mE^l^) ^ i^^lG' - ^S'^GW^) = , (11.38) 



^ - ^SlG\ - {^\G'^ - Is'ydW^) = . (11.39) 

With the exponential parameterization, the formalism of the theory can be expressed in an elegant form. Let us go 
to the rules of differentiation of exponential matrix functions 

f„„x™t - i^u , x% = x^^it ,1 ■ (11.40) 

Taking into account pi.40p . we get the quantum contribution to the full connectivity pi.33p as follows 

I {-i^U - i'U + rfe„x'"vi;„) + i [sl^^J.,k + s'J,, - . (ii.4i) 



^zk 2 



Formulas (|II.35P could be rewritten as follows: 

Xi.-(exp^)l, VI =e*^', 



S= f (Cgrav - - ^Xie*e.^ , (11.42) 



1 1 



1>C OK \ Z 

As is seen from pi.42[) . for the exponential parameterization, the non-polynomial structures of quantum theory of 
gravity have been completely reduced to the factorized exponents'^. 



^ We are using the standard definitions. Matrix functions are defined by their expansion into power series as any operator functions: 

t/(y) = J^c„y" . 

n 

. The derivative of n— th degrees of matrix by the same matrix is defined as 

dV 

The derivative by numerical (non matrix) parameter z is 

dz dz 

If matrix function 0(V) and its derivative W = dO" /dV are elementary functions, then 

dz dz 

Formulas ||II.40| | — l|II.42| l are the consequence of these definitions. It worth to mention, that in matrix analysis in all intermediate 
formulas one should be careful with the index ordering. 
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The explicit form of the tensor, in the terms of which the self-consistent system of equations could be written is as 
follows 



1 ... 



j^lm (ik _lk \ _ xkrnjl _^ }_gk ( x^inU ^ X'^J," 1 



(11.43) 



1 r 

4 



Let us introduce the following notations: 



^(1) 



X 



ik Ak 



(2)-- -5 . 2 



(11.44) 



With use of (jll.44p . let us extract from pi.43p the terms not containing the quantum field, and the terms linear over 
the quantum field: 



Ef = R 



(11.45) 



i{grav) i{ghost) ' 



where 



X 



Ira f 7k 7k \ \rkm 71 \ ^ xf^ f vf^n 7^ i x^^m Tn 

;i) ^^V'^;m " V'm;i ) " ^(1) ^Arm + [X V„;m + ^(l)V'm;, 



(11.46) 



is the EMT of gravitons; 



i{ghost) 



X^' (9,10., + 1 



k \^ml/ 



5fX 



(11.47) 



is the EMT of ghosts. In the averaging of l|II.45p . it was taken into account that j^*) = by definition of the 

quantum field. Averaged equations for the classic fields pi.SSp take form of the standard Einstein equations containing 
averaged EMT of gravitons, renormalized by ghosts: 



mE'^m ^ Rf - -5fR - dfxA - ><{^ft\^) = . 
Quantum dynamic equations for gravitons (|II.39P could be rewritten as follows: 

As is seen in the equations (III.49p . in the theory of gravitons all nonlinear effects are in the difference between the 
EMT operator and its average value. System of equations pi.48p . pi.49[) is closed by the quantum dynamic equations 
for ghosts, which could be also written in 4D covariant form: 



(11.48) 



(11.49) 



{X'''^,k)■^ - 0, iX'''0■,k)■^ = 

Equations pi.SOp provide the realization of the conservative nature of the ghosts' EMT: 



(11.50) 
(11.51) 
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E. Differential Identities 

In the exact theory, which is deahng with the fuU metric, there is an identity: 

Dk {^g'^'Ru - \5':r'Rhn - S'y^K + i [g^' {eaO,, + e.J,i) - S'^g'^'l^e.i] | = , (11.52) 

where is the covariant derivative in the space with metric gik- This identity is satisfied by Bianchi identity and 
by the ghost equations of motion. In terms of covariant derivative in the background space, identity (jll.52p could be 
rewritten as follows: 

Elu - \ (In d),uE^ + f^iE\ - f^E^ ^ E^m " ^1 A' = . (11.53) 

For the exponential parameterization, taking into account (jll.4ip . the expression (111.53^ can be transformed to the 
following form 

A'fc + [Ef - \5^E\^ - . (11.54) 

Identity transformation E^ = {'^\E^\^) + L*^ and the subsequent averaging of (jll.54p yields: 

(*l4'l*>;fe + \WU (^f - \5^e£^ I*> - . (11.55) 

Here we have used explicitly the fact that (^'ji/'f l^*) = 0, (5'|Z/^|\E') = 0, by definition. Next, expression pi.48p is 
substituted into pi.SSp . Taking into account the Bianchi identity and the conservation of the ghost EMT, we obtain: 

As is seen from (jll.56p . quantum equations of motion pi.49|) provide the conservation of the averaged EMT of 
gravitons: 

Take notice, that tensors E^ and in pi.54p . pi.56p are multiplied by the linear forms of graviton field operators 
only. Such a structure of identities is only valid for the exponential parameterization. This fact is of key value for 
the computations in the framework of perturbation theory. The order n of the perturbation theory is defined by the 
highest degree of the field operator in the quantum dynamic equations for gravitons pi.49[) . The EMT of gravitons 
which is consistent with the quantum equation of order n contains averaged products of field operators of the order 
n + 1 (e.g., the quadratic EMT is consistent with the linear operator equation). We see that by defining the order of 
the perturbation theory, we have identity pi.56|) . in which all terms are of the same maximal order of the quantum 
field amplitude: 

^^t^i^)* = (^f'"' - m . (11.58) 

Such a structure of the identity automatically provides the conservation condition pi.57|) at any order of perturbation 
theory^. 



* In the framework of the perturbation theory, any parameterization, except the exponential one, creates mathematically contradictory 

models, in which the perturbative EMT of gravitons i^l'^if^g^^vy^'^ ^'^^ conserved. In our opinion, a discussion of artificial methods 

of solutions of this problem, appeared, for example, if linear parameterization gij, = gi^ + i/^it, is used, makes no sense. The algorithm we 
have suggested here is well defined because it is based on the exact procedure of separation between the classical and quantum variables 
in terms of normal coordinates. We believe there is no other mathematically non— contradictive scheme. 
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F. One— Loop Approximation 



In the framework of one-loop approximation, quantum fields interact only with the classic gravitational field. 
Accordingly, equations (jll.49|) are being converted into linear operator equations: 



1 



Of course, these equations are separated into the equations of constraints (initial conditions): 



and the equations of motion: 



LS]|*)=0, i?|*)=0, Ll\^)=Q 



Li ~ - . 



The equations for ghosts pi.SOp are also transformed into the linear operator equations: 

= , eil^o . 



(11.59) 



(11.60) 



(11.61) 



(11.62) 



In the one-loop approximation, the state vector is represented as a product of normalized state vectors of gravitons 
and ghosts: 



l*> = l*s>l*s/.) 
Equations for macroscopic metric (jll.48p take the form: 

1 



- -6^R = >C ((*JT;'^,_)|vI/,) + {^gh\ft^g,,,,t^\^gh) + S^A) 



(11.63) 



(11.64) 



The averaged EMTs of gravitons and ghosts in equations (|II.64I) are the quadratic forms of the quantum fields. 
Assuming that X''' = .g'*^, X*^. = -i/-'*^, X^^. = /2 in (|lL46|l . (|IL47)) . we obtain: 



(11.65) 



^Uo.*) = (^"*^^' + - ^'^"'^;0 ■ (11-66) 
Quantum equations (jll.59|) . (|II.62[) provide the conservation of tensors (jll.65p . pi.66p in the background space: 

(*9|T;V-)I*^^>; = 0' ^'^9h\T\,nost)\'^9^); k^O. (11.67) 

The ghost sector of the theory pi.59p — pi.67p corresponds to the gauge pi.29|l . Note, however, that all equations of 
the theory, except gauges, are formally general covariant in the background space. That provides a way of expanding 
the class of gauges for classic fields. Obviously, we can move from the initial 4-coordinates, corresponding to the 
classic sector of gauges pi.29[) . to any other coordinates, conserving quantum gauge condition 

Vi^l*) = . (11.68) 
It is not difficult to see, that in the classic sector any gauges of synchronous type are allowed: 

goo = N^{t), 50c = . (11.69) 

where N(t) is an arbitrary function of time. 
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An important technical detail is that in the perturbation theory the graviton field should be consistent with an 
additional identity. In one-loop approximation that identity is obtained from the covariant differentiation of equation 
(llL59|l : 

EE + xK5l) = . (11.70) 

The appearance of conditions (|II.70|) reflects the fact that we are dealing with an approximate theory. As it was 
already mentioned in Section Hi CI the partition of the metric into classic and quantum components, and, respectively, 
the factorization of the path integral, can be only done under the condition that additional constrains are applied to 
the geometry of background space. These constrains are manifested through the structure of the Ricci tensor of the 
background space which should provide the identity l|II.70p for the solutions of dynamic equations for gravitons. In 
the Heisenberg form of quantum theory the additional identity can be written as conditions on the state vector: 

4|*) ^ {Rf + ><M^) ^i,m = . (11.71) 

Status of all constrains for the state vectors are the same and are as follows. If dUMI, (HlMI, (jTLTT]) exist at the 
initial moment of time, the internal properties of the theory should provide their existence at any following instance 
of time. 

While one is conducting a concrete one-loop calculation, there is a problem of gauge invariance of the total EMT 
of gravitons and ghosts. As was mentioned by De Witt after the separation of the metric into background and 
graviton components, the transformations of the diffeomorphism group (jll.4[) can be represented as transformations 
of the internal gauge symmetry of graviton field. In the framework of one-loop approximation, these transformations 
are as follows: 

6i^^-5^rf:^+4 + Tr^ . (11.72) 

The problem of gauge non-invariance is twofold. First, the EMT of gravitons pi.65P is not invariant with respect to 
transformations in (|II.72p . Second, the ghost sector (the ghost EMT), inevitably presented in the theory, depends on 
the gauge. Concerning the first problem, it is known that the operation removing gauge non-invariant terms from the 
EMT of gravitons belongs to the operation of averaging over a quantum ensemble. In the general case of arbitrary 
background geometry and arbitrary graviton wavelengths we encounter a number of problems (when conducting this 
operation), which should be discussed separately. 

In the particular case of the theory of gravitons in a homogeneous and isotropic Universe, the averaging problem has 
a consistent mathematical solution. It was shown in Section FlII Al that removing the gauge non-invariant contributions 
from the EMT of gravitons from the quantum ensemble has been set gauge-invariantly. To address the second aspect 
of the problem, we should take into account that the theory of gravitons in the macroscopic space with the self- 
consistent geometry operates with macroscopic observables. Therefore, in this theory one-loop finiteness, as the 
general property of one-loop quantum gravity, should have a specific embodiment: by their mathematical definition, 
macroscopic observables must be the finite values. This requirement on the theory is realized in the class of allowable 
gauges of full metric consistent with the hypothesis of the existence of macroscopic space and macroscopic observables 
(see Section UlIEp . Gauge (|II.29p used above belongs to this class. 



III. SELF-CONSISTENT THEORY OF GRAVITONS IN THE ISOTROPIC UNIVERSE 

A. Elimination of 3— Vector and 3— Scalar Modes by Conditions Imposed on the State Vector 

We consider the quantum theory of gravitons in the spacetime with the following background metric 

ds^ = g.kdx'dx'' = N^{t)dt^ - a^{t){dx^ + dy^ + dz^) . (III.l) 

In this space the graviton field is expanded over the irreducible representations of the group of three-dimensional 
rotations, i .e. ov er 3-tensor ■i/'f^j), 3-vector -^f^^) = (V'o(«)' '^^(v)) ^nd 3-scalar V-f^^-, = (^o(s)' V'o(s)' V'f(s)) modes. 
Equations (jll.59|) are split into three independent systems of equations, so that each of such systems represents each 
mode separately. The state vector of gravitons is of multiplicative form that reads 

I*,) = |vl/,)|vl/,„)|vl/,) . 
The averaged EMT (jll.63[) is presented by an additive form that reads: 

(*9l^;Va.)l*9) ^ + (*.|t;^.)|*„) + (^.It;^,)!*.) . (111.2) 
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The averaged EMT contains no products of modes that belong to different irreducible representations. This is because 
the equality = is divided into three following three independent equalities 

{^s\i^1is)\'^s) = o, (vi/.|v3f(„)|f„)-o, {^t\^y^t) = o. (III.3) 

Equalities (jlll.SP are conditions that provide the consistency of properties of quantum ensemble of gravitons with 
the properties of homogeneity and isotropy of the background. In the homogeneous and isotropic space, the same 
equalities hold for Fourier images of the graviton field. Therefore, the satisfaction of these equalities is provided by 
the isotropy of graviton spectrum in the k-space and by the equivalence of different polarizations. 

3-tensor modes V'f(i) and their EMT jT-^^^ |^t), respectively, are gauge invariant objects. Gauge non-invariant 

modes V'fj^), "0^(5) a-^e eliminated by conditions that, imposed on the state vector, read 

V'f(„)l*.)=^f(.)l*.>=0. (III.4) 

Note that the conditions (jIII.4[) automatically follow from equations pi.59p and conditions (jll.66[) . As a result of this, 
a gauge non-invariant EMT of 3-scalar and 3-vector modes is eliminated from the macroscopic Einstein equations, 
and we get 

i^rlfH^^l^,) = , (*.|T;^-,)|*,) = . (IIL5) 

The important fact is that in the isotropic Universe, the separation of gauge invariant EMT of 3-tensor gravitons 
is accomplished without the use of short-wave approximation. In connection with this, note the following fact. In 
the theory, which formally operates with waves of arbitrary lengths, the problem of existence of a quantum ensemble 
of waves with wavelengths greater than the distance from horizon is open [l3j . In cosmology, the existence of such 
an ensemble is provided by the following experimental fact. In the real Universe (whose properties are controlled 
by observational data beginning from the instant of recombination), the characteristic scale of casually-connected 
regions is much greater (many orders of magnitude) than the formal horizon of events. The standard explanation of 
this fact is based on the hypothesis of early inflation. Taking into account these circumstances, we do not impose any 
additional restrictions on the quantum ensemble. 

The procedure described above is based on the existence of independent irreducible representations of graviton 
modes only. But in this procedure, gauge-non-invariant modes are eliminated by using of a gauge, i.e. they are 
eliminated by using of gauge-non-invariant procedures. The gauge-invariant procedure of getting the same results is 
presented below. 



1. Elimination of Scalar Modes 



We consider equations pi.59p with the pil.ip background using the conformal time: N = a, Ndt — > a{rj)dri . 
(Symbol "t" belongs now to the physical time, for which N ~ 1). The metric of the 3D flat isotropic Universe is 
conformally similar to Minkowski's metric, and it reads 



ds"^ = 0^(77) [dif - dx^ - dy^ - dz^) = a'^{r])gikdx''dx'' . 

Let us introduce the new variables that can be interpreted as quantum fluctuations of covariant metric 

1 



hf 



(III.6) 



(III.7) 



In terms of variables (|III.7p . the equations pi.59[) read (after calculations of covariant derivatives and Ricci tensor 
components in the pII.6P metric): 



l,k . %,k,l 



1 t.'t 



h^ i - K - 5, [h; 



a 

— n 

a 



-2\- 2^ ) ( UiH^hf Uiu'^h 



Sf 2 ^ 



1 



i~5':{2h^^-h,i) 



{nil 



nrnn'h]" --h] - -SfxAa-'h = , 



where 



n' =n^^ (1, 0, 0, 0) 
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In equation pil.8|) . operations with indexes are defined in the Minkowski space; comas mean derivatives of metric 
fluctuations over the coordinates of the Minkowski space; dashes mean derivatives of scale factor over the conformal 
time 77. In equation (jIII.8|l and further on in this Section, we do not impose any gauge on metric fluctuations. 

As was shown by Lifshitz in 1946 f?\, fluctuations of metric can be expanded in Fourier series over 3-scalar, 
transverse 3-vector and transverse 3-tensor plane waves. Projections of general equations (jIII.9[) onto scalar, vector 
and tensor basis functions lead to three independent systems of equations — for each type of mode separately. Fourier 
images of scalar fluctuations of metric and parameters of gauge transformations are defined as follows 



;il](k) = ^k, /i°(k) = (k)--zfc„xk, /if(k) = ^(Ak + Ak) - 



fc2 



A. 



(IIL9) 



a-2^o(k)=77°(k)=a;k, 



-a^V^ik) = r;"(k) = -zfcVk 



All operations with space indexes of vector-tensor basis in equations (|III.9|1 and further on are conducted with the 
Euclid metric. Gauge transformations pi.69|) for Fourier images of scalar fluctuations read 



(^k - 2 ^ 



6— cjk + 2fc z/k, 
a 



Xk 



Aik ^ Mk 

For brevity, the following notation is used below: 

^k = Ak + Ak, Mk = + 2fc^Xk, 



Xk - {i^i + Wk) 
Ak — 2fc^z/k • 



(III.IO) 



V'k 



<y5k + Mk 



There are two liner combinations of Fourier images of metric fluctuations that are invariant with respect to transfor- 
mations piL10|) nil, which is an important sequence of the theory. They read 



Jk = k'^ipk + - 
a 



-Nu 



a 
a 



Ah 



Jk ^ Jw , 



(III.U) 



The Fourier image of the equation of motion pi.6ip is expanded over the tensor basis. It reads 



In accordance with piI.12p , we obtain 

Li(k) 



ka}f_ 
fc2 



L2(k) - 



(III.12) 



<^k 



\k + \ (Mk - f'k) 



- >cKa^ Vk = , 



(III.13) 



i2(k) = -(iv^ 



ML 



k (fiu) 



k^Nu 



(III. 14) 



To eliminate gauge-non-invariant scalar fluctuations, it is necessary to prove the existence of such initial conditions 
that are independent of gauge and flx the only trivial solution of equations pll.isp . pII.14p . 

As it was mention above, the initial conditions must contain the equations of constrains. After Fourier transforma- 
tions, primary constrains can be written as follows: 



p ^ a' - a'2 
— A^k + -Mk - 3^(^k + ^ ( 3^ - }<Aa^ ] 



1^-,,) = 



(III. 15) 



-A^i^ - -^k 
3 a 



I*.,) - 



-A^4 - -(^k + 2 2 -] Xk 

3 a \ a a'^ 



I*,) = 



(III.16) 
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The condition (jllLlSp is L[5(k) 
equations Qol*) =0 and QqI^) 



0, and piI.16p conditions are obtained from I/°(k) = and ioC*-) ~ ^- Now, 
obtained from pi.71[) are also included in the number of initial conditions: 



2 2— 



— - 2— + >*Aa2 



, 



(III.17) 



t'2 



<y3k + 2 — Xk 

a 



2 h nAa'' Vk 

a 







(III.18) 



All constrains pil.isp — piI.lSp are contained in the original non-gauged equations of the theory. It means that 
their mathematical structure is independent of the choice of gauge. Thus, any quantum ensemble of gravitons in the 
isotropic Universe must satisfy to UII.15\) — IIIII.18\) . Of course, initial conditions cannot be fully defined by these 
constrains because the quantum ensemble is still not actually defined. What is actually defined at this point, are 
ensemble's properties that follow from the isotropy of background. The full determination of ensemble properties can 
be done by imposition of gauge. Such a procedure is gauge non-invariant, and this is the reason why it was disputable 
for many years. (For discussion of the problem of gauge invariant description of scalar fluctuations see, e.g. [33|.) 

To solve this problem, one needs to use the invariant and to impose the following gauge invariant conditions 

on the state vector 



a \ 



Ah 



, 



(III.19) 



Jkl*. 



(III.20) 



The relation piI.19[) can be immediately used as an initial condition because it does not contain higher derivatives. 
Higher derivatives also can be excluded from pil.20[) by non-gauged equation piI.14p via a simple algebraic procedure. 
As a result of these operations, one gets the last initial condition that reads 



2fc2(^k 



fc2 



-iVk 



a V 



Afk 







(III.21) 



Thus, we have initial conditions that are presented by the closed system of algebraic equations (IIII.ISP — piI.19[) 
and equation pil.2ip with respect to Fourier images of metric and their first derivatives over time. Because of its 
homogeneity, this system of equations has the only trivial solution that reads 



TVkl*,) = Ni\^s) = Mkl*.) = Xkl^s) = ^kl^-s) = (^kl^-s) = ^'v\^s) = 



(III.22) 



The substitution of pil.22p into the equations of motion piI.13p , piI.14p shows that higher derivatives are also zeroes 
at the initial instance of time, i.e. 



K\^s) - ML\^s) = 



(III.23) 



It follows from (|III.23p , that conditions (|III.22p defined at the initial instance of time are valid for any future instances 
of time. Thus, scalar fluctuations are excluded from the theory. It is important to emphasize that gauge was not used 
in the procedure that was described above. Scalar fluctuations are eliminated if gauge invariant conditions, which are 
imposed on the state vector, are added to the equations of constrains that are already contained in the theory itself. 



2. Elimination of Vector Modes 



Vector fluctuations and vector parameters of gauge transformations are presented by their Fourier images after 
their expansion over 3-transversal plain waves. They read: 



/i°(k) = 0, /i° (kA) = -/i?(kA) = -5«(kA)xkA , /if (kA) = i [A:„^^(kA) + A:^^„(kA)] iikA 
-a277<,(kA) r;"(kA) = -S'"(kA):/kA , A:„S'"(kA) = . 



(III.24) 
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In pil.24p and further, A is the index of polarization of vector modes. Gauge transformations read 

XkA XkA - t'kA ' "kA — > "kA - i^kA ■ 



(III.25) 



There exists the linear superposition of Fourier images which is invariant with respect to (|III.25|1 transformations. It 
reads 



-fkA = ui ^ - 



kA - XkA 



(III. 26) 



Primary constrains i^j'I't,) ~ 0, Lq|5'„) = and additional constrains Qal'^v) = generate the following initial 
conditions for the state vector 



fc'/kAl^-.) = , 



k^I, 



kA + 4 2— XkA 

a 



a / a 
a \ a 



II 



(III.27) 



XkAl*.) = 



The equation of motion = contains only the invariant in Eq. pil.26[) . It is integrated and reads 

f/ , a' f _ n f _ ^kA 

^kA 1 -"kA — U , -fkA — ■ 



(III. 28) 



According to equations (jIII.27[) and (jIII.28p . the following conditions are imposed on the state vector at the initial 
instant of time 

WkAl^t,) = <aI*«> = XkAl*^,) = . 

They are satisfied at any further instant of time. As it can be seen from the above consideration, to eliminate vector 
fluctuations it is sufficient to take into account only the constrains that exist in the equations of theory. 



B. Canonical Quantization of 3— Tensor Gravitons and Ghosts 

The parameters of gauge transformations do not contain terms of expansion over transverse 3-tensor plane waves. 
Therefore, Fourier images of tensor fluctuations are gauge-invariant by definition. We have 



(III.29) 



/i[!(k) = 0, hl{\,) = -/i^(k) = , hi{\,a) = -^f (ka) = -Q^Mi'^. , 

ko.Q${ka)=0, g^(ka) = 0, 
where a is the index of transverse polarizations. The operator equation for 3-tensor gravitons is 

- (ka)V'ka(i)e^''" , V^k. + SH^ka + ^^k. = , (III.30) 

kcT 

where H = d/a is Hubble function and dots mean derivatives with respect to the physical time t. 

The special property of the gauge used is the following. The differential equation for ghosts is obtained from the 
equation for gravitons by exchange of graviton operator with the ghost operator. It reads 

1.2 

0(i,x) =y 0k(i)e'*'" , ^k + 3i/^k + — ^k = . (III.31) 
k 

Other gauges that automatically provide finiteness of macroscopic quantities in the one-loop quantum gravity (see 
Section lill E[) have the same property. 

Macroscopic Einstein equations (|II.64|) read 

= }<{eg + A) , (III.32) 



2H + 3H^ ^xiA-pg) , 



(III.33) 
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where 

kcr k 

are the energy density and pressure of gravitons that are renormahzed by ghosts. Formulas (jIII.34ll were obtained 
after ehmination of 3~scalar and 3-vector modes from equations pi.65[) and (jll.66p . We also took into account the 
following definitions 

Also we have the following rules of averaging of bilinear forms that arc the consequence of homogeneity and isotropy 
of the background 

(^sl^k^^k'<T'|«'g) = (^'g|?Ak<,V'k^|*g)<5kk"5<,<,' , (^-g/i l^k^'k' | *g/i) = (^gh. I ^k^'k | ^-g/i) <5kk' ■ 

The self-consistent system of equations piI.BOp — (|III.33[) is a particular case of general equations of one-loop 
quantum gravity pi.59[) . pi.62[) . pi.64[) — pi.66[) . In turn, these general equations are the result of the transition to 
the one-loop approximation from exact equations pi.46|) — pi.SOp that were obtained by variation of gauged action 
over classic and quantum variables. To canonically quantize 3-tensor gravitons and ghosts, one needs to make sure 
that the variational procedure takes place for equations pII.30P pIL33P directly. To do so, in the action (111.42^ we 
keep only background terms and terms that are quadratic over 3-tensor fluctuations and ghosts. Then, we exclude 
the full derivative from the background sector and make the transition to Fourier images in the quantum sector. As 
a result of these operations, we obtain the following 



S = J dt (^--^^ - Aa^iV + Lgrav + L ghost 

Lgrav + Lghost = ^ E ( l^^ktr^kCT ^ Nak^^^^^\^„ ) ^ 4~ E ( Jj^^^^ ^ Nak'^0k9i^ 
^ k(T ^ ^ k ^ 



(III.35) 



In pil.35p . the background metric is taken to be in the form of pil.ip . and the N function is taken to be a variation 
variable (the choice of this function, e.g. TV = 1, to be made after variation of action). Here and further on, the 
normalized volume is supposed to be unity, so y = J '^^^ — 1- The terms which are linear over the graviton field are 
eliminated from pil.35p because of zero trace of 3-tensor fluctuations. Variations of action over N and a are done 
with the following averaging. These procedures lead to equations pil.32p . (IIII.33P and expressions pil.34p . Variation 
of action over quantum variables leads to the quantum equations of motion (IIII.30p and pil.3ip . 

In accordance with the standard procedure of canonical quantization of gravitons, one introduces generalized mo- 
menta 



dL ■ 

Then, commutation relations between operators that are defined at the same instant of time read 



^k. = ^ = — V-^, . (III.36) 





_ a' 




TTktJ , V'k'CT' 




V'ko- ' '^k'<T' 



= -iMkk'<^..' . (111.37) 

Formulas pil.36p and pil.37p are presented for the iV = 1 case. Note also that the derivative in pil.36p should be 
calculated taking into account the -i/jj^^ — ip^]^-„ condition. 

The ghost quantization contains three specific issues. First, there is the following technical detail that must be 
taken into account for the definition of generalized momenta of ghost fields. The argument in respect to which the 
differentiation is conducted needs to be considered as a left co-multiplier of quadratic form. Executing the appropriate 
requirement and taking into account Grassman's character of ghost fields, we obtain 
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Second, the quantization of Grassman's fields is carried out by setting the following anti-commutation relations 

^3 r . 



0k , ^k' 



-ih6] 



+ 



kk' 



(III.39) 



-ih6] 



kk' 



Third is the bosonization of ghost fields, which is carried out after quantization of (jIII.39|l . The possibility of the 
bosonization procedure is provided by Grassman algebra, which contains Grassman units defined by relations uu — 
—uu — 1. Therefore, conjunctive Grassman fields can be always presented in the following form 



k ' 



(III.40) 



where t?k is Fourier image of complex scalar field which is described by the usual algebra. The substitution of (jIII.40p 
in (jIII.39|) leads to the following standard Bose commutation relations 



a 



-ihSkk' , -r- 
4>f 



^k , K' 



-ihSkk' 



(III.41) 



The Hermit conjugation transforms one of them to the other. 



C. State Vector of the General Form 



To complete the self-consistent theory of gravitons in the isotropic Universe, one needs to present the algorithm of 
introduction of the graviton-ghost ensemble into the theory. Properties of this ensemble are defined by Heisenberg's 
state vector which is expanded over the basis that has a physical interpretation. Any possible basis is the system of 
eigenvectors of an appropriate time independent Hermit operator. The existence of such operators can be proved in 
a general form. Let us consider the following operator equation which is an analog of operator equations of gravitons 
and ghosts 

jjk + iHijk + -^Vk = . (HI.42) 

Coefficients of equation (|IH.42[) are continuous and differentiated functions of time along all cosmological scales except 
for the singularity. Thus, with the exception of the singular point, the general solution of equation (jHI.42[) definitely 
exists. Below we will show that the existence of a state vector follows only from the existence of general solution of 
equation (|HI.42p (see also [13]). 

Suppose gk, h^ are linear independent solutions to pif .42p . so that their superposition with arbitrary coefficients 
gives the general solution to (|HI.42[) . With no loss of generality, one can suppose that these solutions are normalized 
in some convenient way in each concrete case. From the theory of ordinary differential equations it is known that gk, 
hk functions are connected to each other by the following relation 

Cfc 

gkhk - hkCjk = , (HI. 43) 

where Ck is a normalization constant. The comparison of pil.42|) with (|HI.30|) and pil.3ip shows that solut ions of 
operator equations are presented by the same functions. For operators of graviton field we have 

^ka = ikafffc + Bi,„hk , (HI.44) 

where A\^„, i?ko- are operator constants of integration. Directly from these operator constants, one needs to build the 
operator which gives rise to the full set of basis vectors. 

It is important to keep in mind that commutation property of operator constants A\^u, Ska- and physical interpre- 
tation of basis state vectors are determined by the choice of linear independent solutions of equation (III.4<^ - The 
simplest basis is that of occupation numbers. The choice of linear independent solutions as self-conjugated complex 
functions corresponds to this basis. 

In accordance with pil.43p . if gk = fk, hk — f^ the normalization constant is pure imaginary. Let's take Ck = i, 
so we obtain 
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To build the graviton operator over this basis, one need to carry out the multiphcative renormahzation of operator 
constants taking into account that field is real. This yield 

As result of these operations, we get the graviton operator and its derivative that read 

.Ak. - V4;^(cka/fe+cV./fc) ' ka = ^^{4afk+^-^-Jk) ■ (111.46) 

Standard commutation relations for operators of graviton creation and annihilation are obtained by the substitution 
of pil.46p into pil.37p and taking into account (|III.45p . They read 

[Ckcr , cj,^,]_ = (5kk"5o-o-' , [Cko- , Ck'o-']_ = , [c^^ , Cj^/^/j _ = . (III. 47) 

In accordance with pil.47[) . the operator of occupation numbers fikcr = c^„cit.a exists that gives rise to basis vectors 
jrikcr) of Fock's space. Non-negative integer numbers riko- = 0, 1, 2, ... are eigenvalues of this operator. 

In accordance with pil.24[) . the observables are additive over modes with given ka. Therefore, the state vector is 
of multiplicative structure that reads 

ko- 

where l^'ko-) is state vector of kcr-subsystem of gravitons of momentum p = Kk and polarization a. In turn, in a 
general case, l^'ko-) is an arbitrary superposition of vectors that corresponds to different occupation numbers but the 
same kcr values. Suppose that Cn^^ is the amplitude of probability of finding the ka-subsystem of gravitons in the 
state with the occupation number riko-- If so, then the state vector of the general form is the product of normalized 
superpositions 

kcr rikcr "kCT 

After the bosonization in the ghost sector is done, one gets equations of motion and commutation relations that are 
similar to those for graviton. The same set of linear independent solutions fk , f^ that was introduced for operators 
of graviton field is used for operators of ghost fields. What is necessary to take into account here is originally complex 
character of ghost fields, which leads to -d^ 7^ i^-k- As a result, operators of ghost and anti-ghosts creation and 
annihilation appear in the theory. They read 

79k = v^(akA. + 6lk/fc) > i9+ = \/4^ (a+ A! + 6_kA) ■ (111.49) 
The substitution of pII.49P into pll.4ip leads to standard commutation relations 

[ttk , ak/]_ = <5kk' , [flk , ak']_ = [ftk ' ^k']_ = ' 



kk' , 6k , 6k 



bt^b+ =0, 



(III.50) 



flk , 6k' 




flk , 6j, 




Ok , ^k' 




a+ , 6+ 



Applying the reasoning which is similar to that described above, we conclude that in the ghost sector, the state vector 
of the general form is also given by product of normalized superpositions. It reads 



i^s'^) = n E -^"^ i"k> n E 

k Kk k fik 



(III.51) 



The set of amplitudes C„^^, An^, ^a^i which parameterizes Heisenberg's state vector actually determines the initial 
condition of quantum system of gravitons and ghosts. 
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Formulas pil.48p and pil.51[) can be also used in case when real functions are chosen as linear independent solutions 
of equation (|III.42[) . The justification for this is due to the fact that real linear independent solutions can be obtained 
from complex self-conjugated ones by the following linear transformation 

= ^ ifk + m , hk = ^ih~ m ■ (111.52) 

After transition to the basis of real functions in (jIII.46P and (jIII.49p , we get 



(III.53) 



where 



1 ^ ^ i 

1 i 
9k - ^ («k + b\) , pk = --^ (ak - b\) . 



(III. 54) 



Relations piL54p allow to work with real linear independent solutions and to use simultaneously state vectors pil.48p 
and piI.Sip for the representation of occupation numbers. Note that in the framework of the basis of real functions, 
operator constants are operators of generalized coordinates and momenta: 

P^^, Ql^'a' = -j^kk"5o-o-' , [Pk, <?kv]_ = -«<5kk' . (III. 55) 

To complete this Section, let us discuss two problems that are relevant to intrinsic mathematical properties of the 
theory. First of all, let us mention that "bosonization" of ghost fields is a necessary element of the theory because only 
this procedure provides the existence of state vector in the ghost sector. Mathematically, it is because the structure 
of the classic differential equation pil.42p and properties of its solution pil.45p are inconsistent with the Fermi-Dirac 
quantization. In terms of original ghost fields we have 

9^ = VI^{ak.fk + f^-kfk) , ^k = \/4^ (ak/fc + /3^k/fc) . (III.56) 

Substitution pil.56p into pil.39p and taking into account pil.45p leads to anti-commutation relations for operator 
constants that read 

[ak , ak']+ = -Skk' , [/3k , Pk']^ = Skk' . 

The [/3k , /3k'] _|_ = <5kk' relation can formally be considered as anti-commutation relation for operators giving rise the 
Fermi space of ghost states. There is no such a possibility for ak , cik operators because their anti-commutation is 
negative. If one considers these operators as complete mathematical objects that are not subject to any transforma- 
tions, then it is impossible to build an operator over them that gives rise to some space of states, and this is because 
of non-standard anti-commutation relation. The problem is solved by the fact of the existence of Grassman units 
which are necessary elements of Grassman algebra. At the operator constants level, the bosonization is reduced to 
the following transformation 

ak = uflk , ak = wa+ , ^k = ubu , ^k = ub^ . 

This leads to operators with pil.50p commutation properties. 

The choice of basis is the most significant problem in the interpretation of theory. In the theory of quantum fields of 
non-stationary Universe, the choice of linear independent basis fk, is ambiguous, in principle. This differentiates 
it from the theory of quantum fields in the Minkowski space. In the latter, the separation of field into negative 
and positive frequency components is Lorentz-invariant procedure [s^]. A natural physical postulate in accordance 
to which the definition of particle (quantum of field) in the Minkowski space must be relativistically invariant leads 
mathematically to fk = (2u;fc)~^/^e~*"''*. In the non-stationary Universe with the pil.6p metric, the similar postulate 
can be introduced only for conformally invariant fields and at the level of auxiliary Minkowski space. At the same 
time, the graviton field is conformally non-invariant. This can be seen from the following. Using the conformal 
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transformation = yk/o- and transition to the conformal time drj = dt/a, one can see that equation pil.42p is 
transformed to the equation for the oscillator with variable frequency that reads 



]yk 







(III.57) 



Effects of vacuum polarization and graviton creation in the self-consistent classic gravitational field correspond to 
parametric excitation of the oscillator (|III.57p . 

The approximate separation of field on negative and positive frequency components is possible only in the short 
wavelength limit. Regardless of the background dynamic, linearly independent solutions of equation pil.57p exist, 
and they have the following asymptotes 



fk 



-ikrj 



2k 



fk 



Jkrj 



2k 



k^ 



> 



(III. 



Effects of vacuum polarization and particle creation are negligible for the subsystem of shortwave gravitons. In this 
sector, quanta of gravitational field can be considered, with a good accuracy, as real gravitons that are situated at 
their mass shell. The conservation of the number of such real gravitons takes also place with a good accuracy. In the 
shortwave limit, choosing linear independent solutions of the pil.58[) form, occupation numbers rikg. are interpreted 
as numbers of real gravitons with energy ek = hk/a, momentum p — Kk/a and polarization a. The possibility of 
such an interpretation is the principle and the only argument in favor of choice of this basis. For the subsystem of 
shortwave gravitons, initial conditions are permissible not in the form of products of superpositions but in the form 
of products of state vectors with determined occupation numbers. In accordance with the usual understanding of the 
status of shortwave ghosts, their state can be chosen in the vacuum form. The gas of shortwave gravitons is described 
in more detail in Section flV Al 

In the k^ ~ W /'-'■I vicinity, there is no criterion allowing a choice of preferable basis. It is impossible to introduce 
the definition of real gravitons in this region because there is no mass shell here. This is the reason why we will use the 
term "virtual graviton of determined momentum" in discussions of excitations of long wavelengths. Under the term 
"virtual graviton" we mean a graviton whose momentum is defined but whose energy is undefined. Each set of linear 
independent solutions corresponds to the distribution of energy for the determined momentum. This distribution can 
be set up, for example by the expansion of basis function in the Fourier integral. Thus, the choice of basis is, at the 
same time, the definition of virtual graviton. One needs to mention that different sets of probability amplitudes C„^^ 
correspond to different definitions of the virtual graviton for the same initial physical state. Note also that limitations 
that are defined by asymptotes (|III.58|1 do not fix basis functions completely. 



D. One— Loop Finiteness 



The full system of equations of the theory consists of operator equations for gravitons and ghosts (|III.30p . pll.3ip . 
macroscopic Einstein equations (|III.32[) . pil.33|l and formulas pil.341) for the energy density and pressure of gravitons. 
The averaging of pil.34p is carried out over state vectors of general form pil.48p and piI.Sip . The one-loop finiteness 
is satisfied automatically in this theory. The finiteness is provided by the structure of ghost sector, and it is a result 
of the following two facts. First, in the space with metric pil.ip the ghost equation pil.3ip coincides with graviton 
equation pil.30p . Second, the number of internal degrees of freedom of the complex ghost field coincides with that of 
3-tensor gravitons. We will show this by direct calculations. 

Let us introduce the graviton spectral function which is renormalized by ghosts. It reads 



Zero and first moments of this function are the most important objects of the theory. They are 



(III.59) 



k \ cr 



(III.60) 
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The energy density and pressure of gravitons that are expressed via moments (IIIL60p can be obtained by transfor- 
mations identical to pil.34p with use of equations of motion pil.30p and pil.3ip . They read 

"^^ = l^^ + i^- -^^-^^ + 1^^- (111.61) 
D^Wo + 3HWa . 

In addition, the following relation between moments is derived from equations of motion 

D + 6HD + AWi + IQHWi = . (III.62) 
This relation ensures that the graviton energy-momentum tensor is conservative: 

eg + 2,H{eg+Pg) = Q 

As it was shown above, field operators can always be chosen from the basis of complex self-conjugated functions that 
are the same both for gravitons and ghosts. One needs to also mention that the interpretation of short wave gravitons 
as real gravitons determines the asymptotic of basis functions (see piLSSI) ). After the commutation of operators of 
creation and annihilation are done, graviton contributions to the moments of the spectral function Wn, f^ = 0, 1 can 
be presented in the following form 

Wr.(^grav) = E ^ E ' '/^k-^k. | = 8><;i^ ^Jt.fk + 

2„ ' ^ ' (111-63) 



k ° 



In the right-hand-side of piI.GSp , the first term is the functional which is independent of the structure of Heisenberg 
state vector. It reads 

wZa.) = 8-^E "^r.h = ^f k^-^'nhdk . (III.64) 

k 

The integral pil.64p describes the contribution of zero oscillations whose spectrum is deformed by macroscopic 
gravitational field. Asymptotic piI.SSp shows that this integral is diverges. In such a situation, the usual way is to 
use regularization and renormalization procedures. As a result of these operations, quantum corrections to Einstein 
equations appear. These corrections are the conformal anomalies and terms that came from Lagrangian B? ln(i?/Ag) 
where Ag is a scale parameter that comes from renormalization (see Appendix IXII A| ) . The theory that we present here 
does not use such operations. There is a contribution of ghost zero oscillations in the moments of spectral function. 
Its sign is opposite to pil.64p . It reads 

■.2n , " . (111.65) 



The observables pil.6ip are expressed via sums Wn{grav) + Wn{ghost)- In those sums, the exact graviton-ghost 
compensation takes place in the contribution from zero oscillations. 

The final expressions for the moments of spectral function are obtained by using the explicit form of state vectors 
IIII.48P and piI.Sip . They read 

}^2n 

= 8><;iE — (^kl/feP + U^f + U^fl) , (III.66) 



k 



where 



oo 



^k = E E \Cn.AWa- E E \^r.AW (111.67) 
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and 



(oo \ / oo 

nk=0 / \nk'=0 



(III.68) 



are spectral parameters. They are defined by initial conditions for the chosen normalized basis of linear independent 
solutions of equations pil.42p . (For sake of brevity, in (|III.68[) and below we use the following notation k' = — k, a' = 
— (7.) Note that the relation pil.66[) does not contain divergences. Divergences in the relation pil.661) may appear 
only because of non-physical initial conditions. The spectrum of real gravitons that slowly decreased for /c — > cx) is 
an example of such a non-physical initial conditions. 

The spectral function pII.59P depends of three arbitrary constants as it is averaged over the state vector of general 
form. It reads 

iyk-8>^ft(7Vk|/,|2 + [/*/*2 + c/k/2) . (III.69) 

In pil.69p , the basis of normalized linear independent solutions contains information on the dynamics of operators of 
graviton-ghost field; integration constants iVk, ?7^, t/k contain information on the initial ensemble of this field. Due 
to the background's homogeneity and isotropy the moduli of the amplitudes and average occupation numbers do not 
depend on the directions of wave vectors and polarizations: 

oo oo oo 

K(3)> = E I'^'ik.P'^ka ' K-(g/i)> = E l-^»kP'^k , (^-te/i)) = E l^"kP^k • (III.70) 

)ik<T = l nk = l nk = l 

Phase of amplitudes, in principle, may depend on the directions and polarizations. One must bear in mind that in the 
pure quantum ensembles, for which the averaging over the state vector is defined, phases of amplitudes are determined. 
If the phases are random, then the additional averaging should be conducted over them, which corresponds to the 
density matrix formalism for mixed ensembles. The question of phases of amplitudes is clearly linked to the question 
of the origin of quantum ensembles. In particular, it is natural to assume that the ensemble of long-wavelength 
gravitons arises in the process of restructuring graviton vacuum. This process is due to conformal non-invariance of 
the graviton field and can be described as particle creation. In this case, there is a correlation between the phases 
of states with the same occupation numbers, but mutually opposite momenta and polarizations: the sum of these 
phases is zero. 

If the typical occupation numbers in the ensemble are large, then squares of moduli of probability amplitudes are 
likely to be described by Poisson distributions. For this ensemble we get 



Cuk,, = \/^K(g)]exp(i(/3„kJ , Ak = y^-P["/c(gft)]exp(ix«k) > ^«k' = \/^[»^fe(s/i)] exp(iXftk/ ) : 



(III.71) 



^K] = ^^exp(-(nfe)) . 

The substitution of piI.TOp . pil.71|) to pil.67p . pil.681) leads to 

N^^ = Nu = 2(nA;(g)) - {nk(gh)) ~ {nk(gh)) ■ (III.72) 



(III.73) 



Ag) ^ 1 a(9'») 1 



where 



Cf ^e''^'= ^ E ( E^["M9)] exp(iv«k. - *<Pnk+i,J ) I E -P["fc(9)] exp(j</?„^,^, - i</?nk'+i,.' ) | ' 
\ "k / \ ftk' 



(III.74) 
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Limit equalities = Ck =1 are satisfied if the phase difference between states of neighboring occupation numbers 
does not depend on values of occupation numbers. It is also easy to see that piL72p and pil.73p apply, with somewhat 

different definitions, to any ensemble with and C,\?^^e^^'' parameters. 

We already mentioned above that different basis functions that correspond to different definitions of the virtual 
graviton can be used for the same initial physical state. Limitations due to the prescriptions on the asymptotic 
expression piLSSp allow to fix only asymptotic expansions of basis functions for k oo. These expansions can be 
used, however, only for description of shortwave modes (Section IIV A|) . Meanwhile, all non-trivial quantum gravity 
phenomena take place in spectral region where characteristic wavelengths are of the order of the horizon scale. The 
choice of basis functions to describe these waves is not unique, and the set of amplitudes of probability C„^^ depends 
significantly on this set. At the level of equations pil.72p . piL73p . the ambiguity in the definition of the virtual 

graviton reveals itself in the ambiguity of values of parameters (fife(g)) and Clf'^e^'^'' . Similar ambiguity exists in the 
ghost sector. Two conclusions follow from that. First, it is necessary to work with the state vector of general form, at 
least during the first stage of the study of the system that contains excitations of long wavelengths. Concretization 
of the amplitudes C„^^ is possible only after using of additional physical considerations that are different for each 
concrete case. Second, a theory would be extremely desirable which is invariant with respect to the choice of linear 
independent solutions of equation pn.42[) . and, correspondingly, is invariant with respect to the choice of amplitudes 
of probability C„^^ defining the structure of Heisenberg's state vector, respectively. In Section |V] we will show that 
such a formulation of the theory exists in the form of equation for the spectral function of gravitons renormalized 
by ghosts. The mathematically equivalent formulation of theory exists in the form of infinite BBGKY chain or 
hierarchy where joint description of gravitons and ghost is carried out in terms of moments of spectral function 
Wn, n = 0, 1, 2, N ^oo. 



E. Class of Legitimate Gauges 

The theory presented above is actually the result of transformations of equations which are set up by the original 
gauged path integral pi.l|) . Mostly, these transformations are mathematically identical. There are only three issues 
of the theory that are absent in the original integral. First, the hypothesis of the existing of the classical spacetime 
of definite but self-consistent geometry is introduced. Second, in the self-consistent system of classical and quantum 
equations, the transition to the one-loop approximation is made. Third, a gauge is chosen, which automatically 
provides the one-loop finiteness of self-consistent theory of gravitons in the isotropic Universe. 

The first and second issues were discussed in the process of incorporating these into the theory formalism. We 
will return to an additional discussion of these in Section |Xl . The gauge can be chosen arbitrarily but must be 
mathematically consistent. Moreover, the gauge is necessary to strictly define the path integral as a mathematical 
object. Nevertheless, the use of a gauge unavoidably raises a question about gauge dependency of the results obtained. 

Some of the results of the theory presented above are gauge invariant or can be obtained by use of any gauges. All 
procedures used to build the theory in Section [ill can be done by means of an arbitrary gauge. The general structure 
of theory and its key property (consistency of the system of classic and quantum equations in each order of the theory 
of perturbations over the amplitude of gravitational field) are gauge independent. In Section IlIIl the elimination of 
scalar and vector fluctuations, quantization of 3-tensor gravitons and the construction of graviton state vector of the 
general form are also gauge invariant operations. Only the ghost sector consisting of equations of motion and ghost 
EMT are gauge dependent. Actually, however, there is no arbitrariness in the choice of ghost sector. 

There is almost no doubt that the future theory unites the gravity with other physical interactions must be finite 
to all orders of theory of perturbations. From this point of view, the one-loop finiteness of self-consistent theory 
of gravitons is a manifestation of fundamental properties of quantum gravity and its generalizations. It is fair to 
say that the one-loop finiteness of quantum gravity can be considered as a prototype of the future comprehensive 
theory. From the structure of quantum gravity itself without fields of matter, it follows that in the self-consistent 
theory of gravitons (the simplest version of this theory) the one-loop finiteness can be only achieved by graviton- 
ghost compensation of diverged contributions to macroscopic observables. In this version of theory, there are no 
other mathematical algorithms that are able to provide the one-loop finiteness. The important fact also is the gauge 
invariance of graviton contribution to divergences. The mathematical structure of the ghost sector that is able to 
compensate these divergences follows directly from this fact. In the self-consistent theory of gravitons in the isotropic 
Universe, which is one-loop finite, the ghost equations of motion and ghost EMT must agree with that presented in 
this work. 

We propose the following statement as one of main constructive principles of the theory. The one-loop quantum 
gravity as a theory of gravitons in the macroscopic spacetime with the self-consistent geometry must be finite by 
definition and, hence, by construction. In the framework of this postulate, we present physical principles that choose 
the class of legitimate gauges that automatically provide one-loop finiteness of macroscopic quantities. 
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It directly follows from the definition of the original path integral that the gauge is imposed over the full 

metric. This means that the background and fluctuations are considered in one and the same reference frame if the 
transition to the self-consistent theory is made, which operates with both classical and quantum variables. The 
question concerning the class of legitimate reference frames, i.e. the class of legitimate gauges, inevitably arises 
in calculations. It is because the self-consistent background geometry is actually a non-arbitrary. It was already 
mentioned in Section III CI that the factorization of measure of path integral (outside of short wave approximation) 
is possible only if the background and fluctuations belong to different representations of group of symmetry of the 
background geometry. Reference frames where the high symmetry of background spacetime is displayed clearly stand 
out both physically and mathematically. It is exactly in these reference frames that the separation of non-physical 
degrees of freedom from physical ones can be done on the basis of simple classification of quantum fluctuations over 
irreducible representations of the symmetry group of the background geometry. The class of such reference frames is 
defined by the gauges which are covariant (form-invariant) with respect to transformation of the group of background 
symmetry. 

In the self-consistent theory of gravitons in the homogenous and isotropic Universe we deal with the symmetry 
of background spacetime with metric ()III.1|) . The form of the (jllLip interval is conserved under 3D rotations with 
parameters with arbitrary dependence on time and under arbitrary transformation of the time itself 

dx°'' = n°'p{t)dx^ , t' = f{t) , (III.75) 

where ^"^it) is the Euler matrix, with different angles of rotations at different moments of time, and f{t) is an 
arbitrary function. Below is the explicit form of two legitimate gauges that are form-invariant in respect to the 
transformations (jIII.75p 



35™ = B{t) , = , (III.76) 



There are of course also some non-legitimate gauges. For example, in the case of \/—g<f^ — B{t), d\/ —gg""^ / dx^ — 0, 
there is no inverse operator {M\)~^ . 

To build the theory in Sections |TT] and IIIH we used the pil.76p gauge with B = 1 that leads to a simple and elegant 
ghost sector providing one-loop finiteness automatically. (Recall that, in final equations of self-consistent theory, 
there is a possibility to make a transition from geometrodynamical time which is defined by condition pi.29[) to any 
other time coordinate (see (|II.69I) in Section III Fp ). The ghost sector corresponding to the (|III.77p gauge is of a much 
more complex structure. If however we go over to the one-loop approximation and take into account the explicit form 
of background metric (jllLip , then we get the following simple expression 

det\\Ml\ \ ^ det\\^kV^9"''^^\ \ x det||7VaA|| x det||iVaA|| x det||iVaA|| , (III.78) 

where A is the Laplace operator in the 3D Euclidian space. As it follows from pil.78p . one of the ghost fields is 
described by the Klein-Gordon-Fock equation, and three other ghost fields satisfy the Laplace equation Ax = 0. 
The singular sources of the fields x ^-re absent, therefore the unique constrained solution of the Laplace equation is 
actually a trivial case x = 0. Thus gauge pil.77p . as well as pil.76p . creates a unique non-trivial ghost field. In 
the framework of one-loop approximation, the EMT and equations of motion of this field coincide with pi.62p and 
pi.66p . and the conditions of the quantization for the isotropic space coincide with pil.39p . 
In quantum gravity, the harmonic gauge is often used in calculations. It reads 



With regard to the self-consistent theory of gravitons in the isotropic Universe, note that the group of gauge invariance 
pll.78p contains within it a subgroup of 3D rotations parameters of transformation which is independent on time. 
This subgroup is sufficient to calibrate the background 

as well as to separate gauge invariant 3-tensor gravitons from scalar and vector fiuctuations. In the framework of 
harmonic gauge, the ghost sector in the one-loop approximation reads 



det ||M\|| = det WdkV^g'^d^W x det p,y3^g'*=a,|| x det || 9,-^^5*^- 5, || x det \\dkV^g''d,\\ , 



a=0 



(III.80) 
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Each of the four ghost fields, appearing in pil.80[) . has the dynamic properties of complex scalar field with Grassmann 
algebra. One of these fields compensates graviton vacuum divergences, and the other three fields generate their own 
divergences of the same type. The absence of one-loop finiteness for macroscopic observables in this version of the 
theory can be linked to the fact the gauge pil.79[) does not have a property of form-invariace with respect to 03{t) 
transformation. 

Thus, we are ready to formulate the following proposition. For the theoretical and experimental study of the 
non-stationary isotropic Universe, it is preferable to use only such gauges (reference frames) that satisfy to the 
following conditions. First, they must be invariant with respect to transformation of the symmetry group { III. 75\ ) of 
the background spacetime UII.1\) . Second, they must automatically produce the ghost sector ill. 6^) and llll.66\j which 
provides the one-loop finiteness of the graviton theory in such space. 

This statement is of the clear interpretation. The first part of the statement proposes to eliminate fictitious 
inhomogenities and anisotropy that are due to the motion of the reference frame itself. The second part of the 
statement proposes to organize calculations in such a way that allows not to take into account the fictitious renorm- 
group evolution of observables. Results obtained in the framework of self-consistent theory of gravitons for the 
isotropic Universe are independent on choice of gauge chosen from this class. 

Using gauges that do not provide the one-loop finiteness leads to a fatal internal inconsistency of self-consistent 
theory of gravitons. This is the aforementioned fictitious renorm-group evolution of observables. Such evolution 
is a direct consequence of renormalization of divergences. To remove divergences, in the process of regularization, 
one needs to introduce counter terms with the tensor structure in the Lagrangian of theory which are absent in the 
original Einstein theory. After modification of Einstein's theory an additional term of the type of B? / f^ appears in 
the Lagrangian of the new theory. This term contains a new fundamental constant 1//^. The essence of contradiction 
is in the fact that after this modification, one can not save the original properties of the graviton field which actually 
created a new post-renormalization Lagrangian (see Appendix IXIIBp . To be free of internal inconsistency, non- 
renormalizable quantum gravity (supergravity) must be finite. The one-loop finiteness off the mass shell (supported 
by graviton-ghost compensation of didvirgences) must be considered as a prototype of properties for future theory. 

Thus, the postulate of one-loop finiteness off the mass shell must be accepted as a condition of internal non- 
inconsistency of the self-consistent theory of gravitons. We can conclude also that the one-loop finiteness (and, 
accordingly, the use of gauges of the type (jin.76|) . (jin.77[) ) automatically allows direct participation of ghosts in the 
formation of macroscopic observables. 



F. Instability of Trivial Vacuum and Evolution of Graviton— Ghost Ensembles 

Let us show that the state vector set up at an instant of time taken at random on the cosmological scale cannot be 
that of the trivial vacuum. The term "trivial vacuum" refers to the state of gravitons and ghosts with zero occupation 
numbers regardless of their wavelengths. 

The operator equation for gravitons pil.30p describes the quantum oscillator of variable frequency. After the 
conformal transformation, we get 

(IIL81) 

a 

Using the method of Ref. [35l |. we introduce two new operator functions c\^a{il) and c^k-o-l'?) instead of the original 
operator function (j)kai'>])i with the following additional condition imposed on these functions 

(IIL82) 

Note that (|IIL82|) can be regarded as one of definitions of graviton that is possible if ljI > 0. Substitution of piL82[) 
into pil.37p shows that new operator functions must satisfy to the following time-conserved commutation relations 

[ckAv) , c+^,iT^)]_ = Si^k'Saa' , [ckaiv) , ^k'.' = , [5+,(ry) , 5+,, (??)]_ - . (IIL83) 



35 



Relations (|III.83|) are consistent with general properties of solutions of operator equations. The first equation follows 
from pil.82p and connects derivatives of new operator functions. Substitution of (|III.82[) to (|III.81[) gives the second 
equation. So, the system of equations now reads 

(111.84) 

2ujk 

The commutators (jIII.83[) are the integrals of motion for the system (jIII.84p . This fact demonstrates the consistency 
of the canonic quantization procedure with the operator dynamics. 

In the ghost sector, conformal transformations of ghost fields are done together with the extraction of Grassman 
units 

" " (III.85) 

Analogous to pil.82|) . the operator structure and additional condition are presented by the following relations 



(III.86) 

From pil.86[) . we get the following system 

(III.87) 

a+' = 7fe6-ke-2' ^ "^-'^'^ , 6'_k = Ik&i^' ^ ■ 

Commutators are obtained from pil.50p by replacement of operator constants Ok, by operator functions 
ak(?7), ^k(T?)- They are the integrals of motion of the system pll.87p . 
The following operators are taken at the initial instant of time 

nkaiva) = c^aMckcrim) i "-k(?7o) = ftk (^o)ak(?/o) , «k(?7o) = ^k (?7o)fek(??o) , (III. 88) 

They give rise to the Fock basis. The Heisenberg state vector of the general form is constructed over this basis. 
Macroscopic quantities are expressed via normal distribution functions pil.89p and anomalous distribution functions 
PII.QOP that read 



^k'^W = ^E(*sl^.WSk.W + 5^k-.('?)5-k-.(^)l^.) . 



(III.89) 



(III.90) 



Inhomogeneous ordinary differential equation of the third order for pil.89p functions (taken from equations pil.84p 
and pil.87p ) can be derived, taking into account commutation relations 

-2^k-, - (7fc^^k) + 87^Fk + 4^— = 47^^ ~ 87' • (III.91) 



dif 7fe dry d-q 7^ dr] drj 7fe drj 
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Anomalous and normal distribution functions are connected by the following equations 



di] 



(III. 92) 



7fe (i^k + l)e2^^"'='^'' . 



drj 



dr] 



If gravitons and ghosts are defined by conditions pil.82p and pil.86p at an arbitrary instant of time, the equation 
piI.Qip (or the system of equations (1111.92^ ) describes the evolution of quantum ensemble. The normal distribution 
functions F^{ri), F^^{ri) are the averages of graviton and ghost occupation numbers as functions of time. If the 
initial state is defined as trivial vacuum, then F^^rjo) = 0, F^^rj^) = 0. The inhomogeneity of equations pil.91[) and 

pll.92p leads to instability of trivial vacuum: for 771 > 770 it clearly leads to F^{'qi) ^ 0, F^{rii) ^ 0. In addition, the 
equations pil.9ip and pil.92p show that the graviton-ghost ensemble of the general form inevitably appears in the 
process of evolution. We get the product of normalized superpositions at the instant of time rj = rjo for an arbitrary 
state vector that is built by operators piI.SSp at the instant of time 771 > 7?o. As a matter of fact, such a state vector 
can be obtained automatically. To construct it, one needs to use the following operators that are solutions of operator 
equations pil.82p and piI.SSP 



Any different definition of gravitons and ghost leads to the same conclusion. 

Thus, states with graviton and ghost zero occupation numbers are degenerative and unstable (for a given deter- 
mination). During the evolution of the Universe, such states can appear only incidentally at a single specific instant 
of time. The status of ghosts as a quantum field with non-zero occupation numbers results only from the internal 
properties of the theory. First, there are no ghostless gauges in the quantum gravity. Second, the one-loop finiteness 
determines dynamic properties of ghost field, specifically its conformal non-invariance and zero rest mass. Third, in 
cosmology (in distinction to the scattering problem), the trivial ghost vacuum is unstable in principle. 

The properties of quantum gravity listed above differentiate the gravity from the Yang-Mills gauge fields. To em- 
phasize these differences, let us mention that the Yang-Mills theory has the ghostless gauges. This fact deprives ghosts 
of status of physical objects. Furthermore, in the 5'-matrix theory, all observables are calculated for asymptotical 
states in which all objects of the theory are situated at their mass shells. The vacuum and particles in these states 
are unambiguously defined, and the vacuum is stable. It allows defining ghost states as vacuum states if ghost gauges 
are in use. Note that the elimination of ghosts from the physical sector of the theory is the procedure of restoration 
of gauge invariance of ^'-matrix. This statement is relevant to both the Yang-Mills theory and theory of graviton 
scattering. The theory of gravitons in the cosmological space of self-consistent geometry is different. Its specific 
feature is that there are no asymptotical states in the non-stationary Universe; vacuum is unstable and conformal 
non-invariant gravitons and ghosts of zero rest mass are situated off the mass shell. In the ultra-short wave limit 
only, ghosts can gain the status in cosmology similar to the status they have in S'-matrix theory. This is because 
all non-stationary effects are negligible in this limit in the Universe (Section IIV Al) . There is no physical or/and 
mathematical basis to exclude ghosts from the set of physical objects of the theory if we deal with long wavelength 
modes. Direct participation of long wave ghosts in the formation of macroscopic observables leads to creation of 
physical states that have no analogies in the classic theory of gravity (Sections IIV Bl and IVi p . We name these states 
as "vacuum graviton-ghost condensate" . 

Note, that the ghosts and ghost condensates that might play a possible role in the formation of properties of the 
earl y U niverse appear in some generalizations of the gravity theory and were discussed in many papers (see e.g. 
[36l ISTI. [38I \3^). In this paper, we emphasize the fact that the long wave ghosts gain the status of physical objects 
even in the framework of one-loop quantum gravity. 



Let us consider the gas of gravitons of wavelength that is much shorter than the distance to the cosmological horizon. 
We exclude the long waves from the model. Also, the calculation of observables is done approximately, so that non- 
adiabatic evolution of quantum ensemble is not taken into account. In the framework of these approximations, it 
is possible to save the pure vacuum status of ghosts because their role is just to provide the one-loop finiteness of 
macroscopic quantities. Long wave excitations we will consider in Section [IV Bl 



"-kcr(?7i) = c'^aMckaim) I »^k(»7i) = a+(?7i)ak(?7i) , nk(??i) = ('7i)^k(?/i) , 



(III.93) 



IV. APPROXIMATE SOLUTIONS 



A. 



Gas of Short Wave Gravitons 
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The calculation of observables for the gas of short wave gravitons can be done by general formulas (jIII.61|) , (|III.66|) 
— pil.68|) after the definition of basis functions and the state vector. For the short wave approximation, the full 
asymptotic expansion of basis functions exists that satisfies the normalization condition (|III.45P and asymptotes 
piI.SSp . Of course, to use the method of asymptotic expansions, basis functions must be taken in the following form 



fk 



1 



a\/2el 



k — 



1 



(IV.l) 



'70 



where 



efe = ek{p, p', p", ...) 



a" 

P = 

a 



is a real functional of scale factor and its derivatives. In the short wave approximation, this functional is expanded 
into the local asymptotic series, which satisfies to the following boundary condition^ 



Efc = f^kip, p', p", •■•) fc , P, p', p" 



. 



(IV.2) 



There are no arbitrary constants in this expansion if the (jIV.2p condition is satisfied. 

The following linear ordinary differential equation of the third order with respect to l/e^ functional follows from 
the equation (jIII.57|l for yk = afk, o-fk functions 



k^ + p 



— = 



(IV.3) 



The solution of equation pV.3[) satisfying to the asymptotic condition (jIV.2p reads 

1 1 °° 



^k 



(IV.4) 



Powers of Jfc operator from (IIV.4|) are defined as follows 



Jk-V= -; / — — 

A J UJk 



Vn 



(IV.5) 



• 1 = 1 



Jk ■ {Jk ■ 1) 



'^k 



{Jk ■ 1) . 



The integral is calculated explicitly for arbitrary s, so that • 1 is a local functional of p and its derivatives. It 
follows from pV.5[) that a small parameter of asymptotic expansion is of the order of ~ l/k"^. The pV.4[) solution is 
approximate because non-local effects are not included to the local asymptotical series. Calculation of these effects 
is beyond of limits of this method. 

The asymptotic expansion pV.4[) . PV.SP defines the 1/e/c functional, and hence, it defines basis functions pV.ip . 
The substitution of pV.ip to pil.66p produces asymptotic expansions of moments of spectral function that read 

A>ch fc^" f 



Ic+ c+ 



(IV.6) 



^ Note that the p^"^[r)) asymptotic exists for cosmological solutions of usual interest. For instance, p'"^(??o) = as = — oo for 
the inflation solution. For »7o = +oo it takes place for the FRW solution for the Universe filled with ordinary matter. 
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State vectors from (IIV.6|1 can be concretized from the general considerations. It was mentioned in Section llll CI that 
such terms as vacuum, zero oscillations and quantum wave excitations are well defined for the p^"-* (77) — > condition. 
Under the same condition, state vectors that are built on basis vectors of the Fock space are easily interpreted. First 
of all, this statement is relevant to gravitons. Eigenvalues ntcr and eigenvectors \n\^a) of nko- — c^^c^a operator 
describe real gravitons in asymptotic states. In the short wave approximation, the concept of real gravitons is valid 
for all other stages of the Universe evolution. Thus, in this particular case, the state vector of the general form can be 
reduced to the product of vectors corresponding to states with definite graviton numbers rikcr = 0, 1, 2, ... possessing 
definite momentum and polarization. It reads 

i*s)=ni"k.)- (IV.7) 

kcr 

In asymptotic states, short wave ghosts are only used to compensate non-physical vacuum divergences. In accordance 
with such an interpretation of the ghost status, we suppose that ghosts and anti-ghosts sit in vacuum states that read 



i*^"')=nniOk>iOk'> . (1V.8) 

k k' 

Averaging over the quantum state that is defined by pV.7p and pV.8|) vectors, we get 

*5|cLcV.I*S> - *5|C-k-.Ck.|*g> = l^^ikl^S'O = I ^-kOk I ^g/, > - 0, 

ko- 

To calculate macroscopic observables in this approximation, it is sufficient to keep only the first terms of expansion 
of moments of spectral function that contain no higher than second derivative of scale factor. In this approximation, 
moments of spectral function read 

ko kcr 

Taking into account pV.Op . we get energy density and pressure of high-frequency graviton gas from piI.Gip that read 

ko ko 

3a4 



Kpg 



Relations pV.9p and pV.lOp arc valid if c? jk^ ~ ^ , \tl\^^ ^ i.e. the square of ratio of graviton wavelength to 
horizon distance is much less than unity. In case of large occupation numbers, these results are of the quasi-classical 
character and can be obtained by the classical theory of gravitational waves [lO| . 

As can be seen from pV.lOp . the high-frequency graviton gas differs from the ideal gas with the equation of 
state p = e/3 by only so-called post-hydrodynamic corrections. In accordance with the approximation used, these 
corrections are of the order of }?H^ <C 1 in comparison with main terms. Thus, the following simple formula can be 
used 

inpg ~ % , Cgi = fcnko . (IV.U) 



yiEg 

a 



ko- 



B. Quantized Gravitons and Ghosts of Super-Long Wavelengths 



In the framework of this theory, it is possible to describe the ensemble of super-long gravitational waves (fc^ ^ 
|a"/a|) by an approximate analytical method. Such an ensemble corresponds to the Universe whose observable part 
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is in the chaotic bunch of gravitational waves of wavelengths greater than the horizon distance. The chaotic nature 
of the bunch is provided by non-zero wave vectors of these waves, so that observable properties of the Universe are 
formed by superposition of waves of different polarizations and orientations in the space. Such a wave system can 
produce an the isotropic spectrum and isotropic polarization ensemble consistent with the homogeneity and isotropy 
of the macroscopic space. 

Such an ensemble of super-long waves can be formed only if the size of causally-bounded region is much greater 
than the horizon distance, which is possible in the framework of the hypothesis of early inflation (or other scenarios 
(see, e.g. [40]). However, the problem of kinematical stability of an ensemble exists even in the framework of the 
hypothesis of early inflation. The case is due to the fact that the ensemble of long waves is destroyed during the 
post-inflation epoch if the Universe is expanded with a deceleration. When long waves come out of horizon, they are 
transformed to the short waves. Below we show that the kinematical self-stabilization of an ensemble is possible in 
the framework of self-consistent theory of long waves. 

Long waves under discussion correspond to virtual gravitons. To describe them approximately, one needs to use 
asymptotic expansions of basis functions over the small parameter ~ fc^. As well as in the case of short waves, the 
basis can be chosen in the representation of self-conjugated functions that are parameterized by the universal real 
functional efc(a). This preserves the definition (|IV.1[) and equation (jIV.3p . However, due to of our interest in the 
asymptotical expansion l/efc(a) over ^ fc^ parameter, it is necessary to rewrite equation pV.3p in the following form 

= _4fcV ( - V . (IV.12) 









'1 




a" 















Let us introduce the geometric-dynamic time dr — drj/o? and the following functional 



oo 

(") 



2e!a2 - ^'^ - X^/*^ 

n— 



Note that the r time coordinate corresponds to the original gauge \J—gg^^ — 1. Equation pV.12p and the spectral 
function (|HL69p now read 

= -4fcV|-(a2a), (IV.13) 



r 



(IV. 14) 



where Tk is a numerical parameter. Its value is unimportant because the constant's contribution to phases of basis 
functions is absorbed by phases of contributors that form vectors of the general form pil.48p and piI.Sip . Observables 
(IIII.6ip are expressed via moments of spectral function. The latter read 



D 



-y 



a" dr^ 

k 



8>ch \ ^ 

k 



dr"^ 



(IV.15) 



The iteration procedure over ~ fc^ parameter for equation pV.13p is constructed accordingly to the following rules 



^fe = 

dT^ 



40) - 



(0) 



= Pu + RkT + QkT^ , 



(IV. 16) 



n ^ 1 . 
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In particular, we get 

~2k^Pkd 



^=-2fc2p,a4 + ... . (IV.17) 



The virtual graviton is defined by integration constants Pfc, Qk^ Rk of the main term of asymptotic expansion. 
Because the efc functional of (|IV.1|) is real (and therefore the s}k^ functional of pV.16p is also real), we obtain 
following inequality 

APkQk-Ri>0, Pk>0, Qk>0. (IV.18) 

The dependence of constants Pj,, Q^, Rk and phase $fe on A; for /c ^ is defined by the finiteness condition for k^W\^ 
and <PW\!^/ dr^ , and taking into account the inequality (jIV.18|) we obtain 

Pk = 0{k-^), Rk - 0{k-^), Qk - 0(A:°), $fe = Oik") . 

The main terms of asymptotical expansions of moments PV.15[) , energy density and pressure of long wave gravitons 
can be obtained from (|IV.16|) for ^ and (|IV.17|) for ^ . They read 



Cg2 Cg3 _ Cg2 Cg3 



(IV. 19) 



where 

Cg2 ^^hJ2 k'Pk (iVk + C/k + f^k) , 

k 

Cgz = xhy\ Qk (iVk + U^ + t/k) . 



(IV.20) 



k 



For the first time, approximate solutions for the energy density and pressure in the (|IV.19[) form were obtained for 
classical long gravitational waves in [111, ■ In the theory of classical gravitational waves [ll|, [l^ , the constants 
of integration Cg2 and Cg^ must be positive. The crucial formal difference between classical and quantum long 
gravitational waves is in the fact that the last ones allow an arbitrary sign of Cg2 and Cg^ (negative as well as 
positive). The physics of this crucial difference will be discussed below (Section IIVC|) . 

In a particular case of (5-type graviton spectrum, which is localized at the region of very small conformal wave 
numbers, PV.19P can be considered as exact solutions. One needs to to go over from summation to integration 

oc 

Y^:-^j^,jd^k...^±,jedk.... 
^ 

After that, these solutions can be obtained by the following limits 

fci 

k^Pk — = const(k), Qk Q — const(k) , 
ai 

(IV.21) 

271-2 

Nu + U^ + Uu-^ -^J^o^{k - kq) , JVo = const{k) , kq ^ . 

In pV.21|) ki and ai are the constants of dimension of conformal wave number and scale factor, respectively. They 
provide the correct dimension to parameter lim k^Pk- 

fe^O 



C. Scenarios of Macroscopic Evolution 



In accordance with (|IV.19p , the system of long wave gravitons behaves as a medium consisting of two subsystems 
whose equations of state are pi = — ei/3 and p2 = £2- But, the internal structure of this substratum cannot be 



41 



determined by measurements that are conducted under the horizon of events. The substratum effect PV.19|) on 
evolution of the Universe, is seen by an observer as an energy density and pressure of the "empty" (non-structured) 
spacetime, i.e. vacuum. The question is: does the graviton vacuum have a quasi-classic nature, or has its quantum 
gravity origin been revealed in some cases? 

Let us review the situation. First, the superposition of quantum states in state vectors of the general form l|III.48[) 
and pil.51|l could be essentially non-classical. Second, the clearly non-classical ghost sector is inevitably presented 
in the theory. Its properties are determined by the condition of one-loop finiteness of macroscopic quantities (Section 
IIIIEp . The ghost sector is directly relevant to the (jIV.lQp solution. Let us consider pil.72p and (jIIL73p . assuming for 
the sake of simplicity {nk[gh)) — (nk{gh))- Parameters of solution (jlV.lQp are expressed via parameters of graviton- 
ghost ensemble as follows 

Cg2=2KhY^k^Pk [{Ukig-,) (1 + Ci'UoS ipk) ~ (Uk^gh)) {1 + Clf''^ COS Xk) , 

^ (IV.22) 

Cg3^2>ih^Qk {nk{g)){l + Ck ^ cosipk) - {nk{gh)){l + Ck''^ cosxk) ■ 

k 

It follows from pV.22p that Cg2 > 0, Cgs > if the graviton contribution dominates over ghosts in the quantum 
condensate. We will name such a condensate " quasi-classical" . Its energy density is positive, and it can be formed by 
usual super-long gravitational waves. If the ghost contribution dominates over gravitons in the quantum condensate, 
then Cg2 < 0, Cgs < 0. Such a condensate of negative energy density has no classical analogy. 

Summarizing the results of Sections IIV Al and IIV Bl we see that in cosmological applications of one-loop quantum 
gravity we deal with the multi-component system consisting of short wave graviton gas gl and two subsystems of 
graviton-ghost condensate g2, g3. Taking into account pV.lip and PV.19|) . we get the following equation for the 
scale factor 

3^ = % + % + % . (IV.23) 

a a* a° 

In the first scenario, the long wavelength condensate is of negative energy, which means that the contribution of ghost 
dominates over gravitons. The evolution of such a Universe is of oscillating type. The solution reads 



2 Cgl ^ V C-gi - 4Cg2C<,3 U\Cg2\ 



(IV. 24) 



"^'^ 2\Cg2\^ 2\Cg2\ 

There is no classic analogy to the solution pV.24p . It can be used for scenarios of evolution of the early quantum 
Universe. In the region of minimal values of the scale factor amin = «i, the 53 condensate bounces the Universe back 
from a singularity. The transition from the expansion to the contraction epoch at the region of maximal scale factor 
dmax = 0,2 IS provided by g2 condensate. Because of correlation of signs of Cg2 < and Cgs < 0, the non-singular 
Universe oscillates. Recent scenarios of oscillating Universes based on condensates of hypothetical ghost fields are 
under discussion in the current literature as an alternative to the idea of inflation (see, e.g. [lO])). Actually, we have 
shown that the same-type scenario is constructed with the standard building blocks of quantum gravity the well- 
known De Witt-Faddeev-Popov's ghosts located far from the mass shell. Thus, a very attractive idea is that one and 
the same mechanism of graviton-ghost condensate formations in the framework of one-loop quantum gravity based 
on the " standard" Einstein equations (without hypothetical fields and generalizations of Einstein's general relativity) 
could be responsible for both Dark Energy effect (see Section lIXp and cyclic evolution of the early Universe (instead 
of inflation). 

The second type of scenario applies if gravitons dominate over ghosts in the condensate of positive energy. The 
solution reads 



2^Cg2{Cg2ai + Cgia^ + Cg3) + 2C g2a^ + Cgl = (2^^;;^^+ Cgl) cxp [ J J • (IV.25) 



The g3 condensate forms the regime of evolution in the vicinity of singularity; meanwhile the asymptote of cosmological 
solution for 77 —> 00 is formed by g2 condensate. Short wave gravitons gl dominate during the intermediate epoch. 
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The ratio of graviton wavelength to horizon distance is constant during the following asymptotical regime 




This means that the long wave condensate g2 forms the self-consistent regime of evolution that provides its kinematic 
stability. 



V. BBGKY HIERARCHY (CHAIN) AND EXACT SOLUTIONS OF ONE-LOOP QUANTUM GRAVITY 

EQUATIONS 

A. Constructing the Chain 

Approximate methods used in Sections IIV Al and IIVBI provide an opportunity to describe only limit cases which 
are ultra shortwave gravitons and ghosts against the background of almost stable Fock vacuum and super-long wave 
modes, forming nearly stable graviton-ghost condensate. Now we are examining self-consistent theory of gravitons 
and ghosts with the wavelengths of the order of distance to the horizon: 



-.2 



H\ \H\ . (V.l) 



When describing modes (|V.1[) . one should keep in mind two factors. First, in the area of the spectrum (jV.ip . 
there are no reasonable approximations, which could be used to solve equations pil.30[) and pil.31[) . if the law of 
cosmological expansion a{t), H{t) is not known in advance. Second, the (jV.ip modes are quasi-resonant. Quantum 
gravity processes of vacuum polarization, spontaneous graviton creation by self-consistent field and graviton-ghost 
condensation are the most intensive in this region of spectrum. From (jV.ip it is also obvious that the threshold for 
quantum gravitational processes involving zero rest mass gravitons and ghosts is absent. These processes at the scale 
of horizon occur at any stage of evolution of the Universe, including, in the modern Universe. 

The theory that allows quantitatively describe quasi-resonant quantum gravitational effects is constructed in the 
following way. For the spectral function of gravitons and ghosts Wu, as defined in pil.59[) . a differential equation 

is derived. For this, the first equation (jIII.30p is multiplied by the i/'ko- (and then by the ip^^), conjugated equation 

piL3p is multiplied by the ipu^ (and then by the V'ko-); and the equations obtained are averaged and added. Similar 
action is carried out with equations for ghosts, after which the equations for ghosts are subtracted from the equations 
for gravitons. These operations yield: 

Wu - 2Fk + 3Wk + ^l^k = , (V.2) 
Fk = -6i/Fk - ^VFk , (V-3) 

where 



a 

Further, equation ljV.2p is differentiated. Expressions for Fk, F-^ via Wk are substituted into the results of differenti- 
ation. For the spectral function the third-order equation is produced 

M>k +9i?W'k + 3 (ij + 6i/") + ^ (w^k + 2i/Wk) = 0. (V.4) 

It is now necessary to draw attention to the fact that Wk(i) is Fourier image of the two-point function, taken at 
t = t': 

W{t, t'; X - x') = WHt, ^Wk{t', x') ~ 2e{t, x)0(i', x')!*) , 

1 f (^-'^ 
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An infinite set of Fourier images is mathematically equivalent to the infinite set of moments of the spectral function 



\,2n 



71 = 0, 1, 2, CXD 



(V.6) 



Therefore, from the equation for Fourier images (|V.4p . we can move to an infinite system of equations for the 
moments. For this, equation (jV.4[) is multiplied by (/c/a)^" followed by summation over wave numbers. The result is 
a Bogoliubov-Born-Green-Kirkwood-Yvon (BBGKY) chain. Each equation of this chain connects the neighboring 
moments: 



D + QHD + 4VKi + IQHWi = 



(V.7) 



B{1, 2) =Wi +lbHWi + 3 {22H'^ + 3ijj W'l + 2 {aQH"^ + 18HH + H^j Wi + m-2. + 2mW-2. = 0, 



(V.8) 



B(n, 71 + 1) =Wn +3(271 + Z)HWr, + 3 (471^ + 127i + 6) i?^ + (27i + 



(V.9) 
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2 (271^ + 971 + 9) i/^ + 6(71 + 2)llB + if ] VK„ + 4W^„+i + 8(71 + 2)i/W„+i = 0, 71 = 2, cx) . 



Equations (|V.7P — (IV. 9P have to be solved jointly with the following macroscopic Einstein equations 



3i/2 = 4^ + 7^1 
16 4 



(V.IO) 



Note that an infinite chain of equations (|V.7p — (|V.9[) contains information not only on the space-time dynamics 
of field operators, but also about the quantum ensemble, over which the averaging is done. The multitude of solutions 
of the equations of the chain includes all possible self-consistent solutions of the operator equations, averaged over 
all possible quantum ensembles. Theory of gravitons presented by BBGKY chain, conceptually and mathematically 
corresponds to the axiomatic quantum field theory in the Wightman formulation (see Chapter 8 in monograph [285). 
Here, as in Wightman, full information on the quantum field is contained in an infinite sequence of averaged correlation 
functions, definitions of which simply relate to the symmetry properties of manifold, on which this field determines. 

In BBGKY chain (jV.7p . (jV.8p and (jV.9[) . unified graviton-ghost objects appear which are moments of the spectral 
function, renormalized by ghosts. The ghosts are not explicitly labeled so that the chain is can be built formally in 
the model not containing ghost fields. Mathematical incorrectness of such a model is obvious only with a microscopic 
point of view because in the quantum theory all the moments of spectral function diverge the stronger, the more the 
moment number is. The system of equations (jV.7p — (jV.9p does not "know" , however, that without the involvement 
of ghosts (or something other renormalization procedure) it applies to the mathematically non-existent quantities. 
The three following mathematical facts are of principal importance. 

(i) In one-loop quantum gravity, the BBGKY chain can be formally introduced at an axiomatic level; 

(ii) The internal properties of equations (F. 7| ) — IIV.10\) provide the existence of finite solutions to this system; 

(iii) Infinite solutions, there are solutions which do not meet the "classic" condition of positiveness of moments 
(see Sections \ V B\ and W C\) . 

It follows from these facts that there should be an opportunity and the need to implement a renormalization 
procedure to the theory. This procedure should be able to redefine the moments of the spectral function to finite 
values, but that leaves them sign-undefined. As it can be seen from the theory which is presented in Sections [TTl and 
lllli in the one-loop quantum gravity such a procedure is contained within the theory under condition that the ghost 
sector automatically provides the one-loop finiteness. 

We found three exact self-consistent solutions of the system of equations consisting of the BBGKY chain (|V.7|) — 
(jV.9[) and macroscopic given below in Sections lVBI and IV CI The existence of exact solutions can be obtained through 
direct substitution into the original system of equations. The microscopic nature of these solutions, i.e. dynamics of 
operators and structure of state vector is described in Sections IVII IVIII 
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B. Graviton— Ghost Condensates of Constant Conformal Wavelength 



In Section IIVBI the exact solution was found for the graviton-ghost condensate, consisting of spatially uniform 
modes (see pV.19|) — (|IV.21I) '). This solution satisfies to the first two BBGKY equations (jVjp . (jV.Sp for an arbitrary 
law of evolution H{t) and under condition that Wn = for n ^ 2. (Recall that in this solution D and Wi must be 
understood as the result of limit transition 0; and equality to zero of higher moments follows from the spatial 

uniformity of modes.) Now we describe the exact self-consistent solutions for the system, in which in addition to 
spatially uniform modes, quasi-resonant modes with a wavelength equal to the distance to the horizon of events are 
taken into account. In terms of moments of the spectral function, the structure of solutions under discussion is 



D = D{g2) + D{g3) + D{g4) , 



Wi ^Wi{g2) + Wi{gA) 



D{92) 



Wi{g2) 



Wn = W„(ff4), n ^ 2 

8Cg2 



(V.ll) 



nf A\ 48C34(i) flo 



Here Cgs, Cg2, Cg4(„), ao are numerical parameters. Restrictions on their values follow from the condition of the 
existence of the exact self-consistent solution. 

The solution is found by using of the consistency of functions (jV.ll[) with the relations arising from the macroscopic 
Einstein's equations (we are discussing model with A = 0): 
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In (jV.12p as well as further, we use notation Cg4(i) = CgA- Functions D and Wi from (|V.11|) transform the equation 
(|V.7[) to an identity. The substitution of Wi and W2 into (|V.8[) . taking into account (IV.12[) . leads to the following 
expression 
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B{l,2)^H-^ 



4(Cg4(2) - Cg4) In - ■^Cg2Cgi + C^^ - 2Cg^2) 



= 



(V.13) 



The infinite chain (jV.Op . in contrast to the equation (|V.8|) . contains moments of spectral functions of quasi-resonant 
modes. Nevertheless, it does result, only including (jV.13|) as a particular case 
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B{n,n + 1) = H— 



4(Cg4(„+l) — Cg4Cg4(„)) In — -^Cg2Cgi(n) + C g^C gi(^n) ~ 2Cg4(„+l) 



= 



n = 2, 00. 

The following relations between parameters follow from (jV.13[) and (jV.14[) 



(V.14) 
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Thus, moments of the spectral function of quasi-resonant modes satisfy to the following recurrent relation 



C, 



Wn+i{g^) = ^Wn{g'i)^\^] W,{gA). 
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(V.16) 



Comparison of (|V.16p with (|V.6p shows that in the exact solution under discussion all quasi-resonant modes have 
the same wavelength A — a/ ^J\Cg4\ = a/ko. In other words, in the space of conformal wave numbers the spectrum of 
quasi-resonant wave modes is localized in the vicinity of the fixed value |k| — fcp. 
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Depending on the sign of Cg4, we get two exact solutions to the macroscopic observables of graviton-ghost media 
in the form of functional of scale factor, 
(i) Oscillating Universe. 

Suppose that Cg4 > 0. In accordance with (|V.15[) . in this case all Cg4(„) > 0. The positive sign of all moments 
Wn{g4:) > suggests that gravitons dominate over ghosts in the ensemble of quasi-resonant modes. We also see 
that the parameter of spatially uniform mode g2 is negative, i.e. Cg2 < 0. As was shown in Section IIV C) signs of 
parameters of 172 and modes are the same, so Cg3 < 0. From this it follows that ghosts are dominant in case of 
spatially uniform modes. The energy density and pressure of graviton-ghost substratum read 

^eg = -\^ + ^±ln^ , ^p^ = _^_%ln^. (V.17) 
a° a"' a a° a-^ ea 

The parameter Cg2 is not explicitly showed up in (|V.17p because it is expressed via Cg4 in accordance with (|V.15p . 
There is an oscillating solution to the Einstein equation — xEg if solutions for the turning points = , amax 
exist, i.e. 

5^!^>e, f^V^Mnf^V. (V.18) 



In the vicinity of turning points energy density is formed by contributions of ghosts and gravitons, which are com- 
parable in their absolute values, but have opposite signs. Far from turning points, graviton quasi-resonant modes 
dominate. Simplifying the situation, we can say that in the oscillating Universe spatially uniform modes have essen- 
tially quantum nature, and quasi-resonant modes allow semi-classical interpretation. 

In the absence of a spatially homogeneous subsystem g3, the infinite sequence of oscillations degenerates into one 
semi-oscillation. Indeed, with Cg^ = the scale factor, as a function of cosmological time, reads 

a(7y)=aoexp(^-^l^^ , C<,4 > . (V.19) 

In accordance with (jV.lOp . the Universe originates from a singularity, reaches the state of maximal scale factor 
^max = ao and then collapses again to singularity, 
(ii) Birth in Singularity and Accelerating Expansion. 

Accordingly to (jV.16p , moments of the spectral function of quasi-resonant modes form an alternating sequence if 
Cg4 < 0. It reads 

W„(g4) = -(-l)"^^%[^ln^ , n = 1, 00 . (V.20) 

It is clear that the result (jV.20p can not be obtained for the quasi-classical ensemble of gravitational waves. The 
microscopic nature of this solution is discussed in Section IVII It is appropriate here to emphasize one more time 
that the theory, which is formulated in the most common way in the BBGKY form, captures the existence of such a 
solution. 

It is not difficult to notice that the solution which we arc now discussing is in a sense, an alternative to the previous 
solution. With Cg4 < 0, parameters of spatially homogeneous modes are positive Cg2 > 0, Cg^ > 0. Thus, spatially 
uniform modes admit semi-classical interpretation, but quasi-resonant modes have essentially quantum nature. The 
energy density and pressure of graviton-ghost substratum are 

^93 , 3|Cg4| a Cg3 \Cg4:\ , ea /Trr,i\ 
>ceg = ^ + §- In — , >^Ps = -| 1- In — . V.21 

Specific properties of solutions to Einstein's equations = xSg depend on initial conditions and relations between 
the parameters of graviton-ghost substratum. First of all, let us mention a scenario that corresponds to a singular 
origin with the strong excitation of spatially uniform modes 

Cg.^O, H>0, &^<e. (V.22) 

4C.g3 

In the case (|V.22|) . the Universe is born in the singularity and fairly quickly reaches the area of large scale factor 
values, where it expands with the acceleration: 

„.|a,|"^*l„.'^l, a»a,.{^Y\ (V,23) 
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Branch of the same sohition, with H < describes the collapsing Universe with a singular end-state. 
Two other scenarios correspond to the weak excitation of graviton spatially uniform modes 

+ 0, > e . (V.24) 

In the case of (|V.24I) , the region of legitimate values of the scale factor is divided into two sub-regions ^ a ^ oi 
and 02 ^ a < CX3 separated by a barrier of finite width 02 > ai. In the sub-region of small values of the scale factor, 
the Universe is born in a singularity, reaches the state with a maximum value of a = oi, and then returns to the 
singularity. In the limit C^^g — > the possibility of such an evolution disappears because of a\ — > 0. In sub-region 
of the large scale factor, the evolution of the Universe starts at the infinite past from the state of zero curvature. 
At the stage of compression, the Universe reaches the state with a minimum value of a = 02, and then turns into 
an accelerated mode of expansion. With Cg3 = 0, this branch of cosmological solution is described by the following 
function of cosmological time 

a(77)=aoexp(^i^^^ , C<,4 < . (V.25) 
Note that degenerate solutions (|V.19p and (|V.25p differ only in the sign under of exponent. 



C. Self— Polarized Graviton— Ghost Condensate in De Sitter Space 

It is easy to find that the system of equations (|V.7|) — (jV.lOp has a simple stationary solution R — const, D — const, 
Wn — const. This solution describes the highly symmetrical graviton-ghost substratum that fills the De Sitter space. 
It reads 

= + , a = aoe"' , 

36 3 , ^ 

(V.26) 

1 

This solution exists both for the A = case and for A ^ 0. The first moment of the spectral function satisfies the 
inequality Wi > —12xA is the only independent parameter of the solution. The remaining moments are expressed 
through by recurrence relations: 

8 n(2n + 3)(n + 3) , 

D = --Wi , Wn+i = ^ ^, -H^Wn , n > 1 . V.27 

3 2(n + 2) 

From (jV.26p and (jV.27p it clearly follows that the solution has essentially vacuum and quantum nature. The first 
can be seen from the equation of state Pg = —Eg. The second can be seen from the fact that the signs of the moments 
Wn+i/Wn < alternate. Another sign of the quantum nature of the effect is contained in the properties of graviton 
spectrum. The first of recurrence relations allows estimating of wavelengths of gravitons and ghosts that play a 
dominant part in the formation of observables 



X^^^J^ = lJl = const. (V.28) 
k \ \W2\ H\ 10 ^ ^ 

As can be seen from (jV.28p . during the exponential expansion of the Universe typical values of k rapidly shift to 
the region of exponentially large conformal wave numbers. The physical wavelength and macroscopic observables are 
unchanged in time. Such a situation occurs if the following two conditions apply. 

(i) In the k— space of conformal wave numbers spectra of graviton vacuum fluctuations are flat; 

(ii) In the integration over the flat spectrum, divergent components of integrals excluded for reason to be discussed 
in Section fVl Bl Observables are formed by finite residuals of these integrals. 

In Section rVlB I we will show that these conditions are actually satisfied on the exact solution of operator equations 
of motion, with special choice of Heisenberg's state vector of graviton-ghost vacuum. Microscopic calculation also 
allows expressing the first moment of spectral function through the curvature of De Sitter space 



(V.29) 
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where Ng is a functional of parameters of state vector, which is of the order of the number of virtual gravitons and 
ghosts that are situated under the horizon of events. Their wavelengths are of the order of the distance to the horizon. 
It must be stressed that the number of gravitons and ghosts Ng is a macroscopic value. 

The order of magnitude of Ng is determined by graviton and ghost numbers in the condensate. Let us emphasize 
that numbers of gravitons and ghosts and hence, Ng parameters are macroscopic qualities. Further down in this 
section it is assumed that the gravitons dominate in the condensate and that the parameter Ng > 0. 

Note that the result (|VT29)) can be easily predicted from the general considerations, including considerations of 
dimension. Indeed, the general formula (jIII.66p shows that the moment Wi is of dimension \Wi\ = [/]~^ ([/] is of 
dimension of length). It also contains the square of the Planck length as a coefficient. Because Wi is a functional of 
the metric, desired dimension can be obtained only using metric's derivatives. It follows from this that Wi = C ■ xhH^ 
where C dimensionless constant that contains parameters of vacuum condensate. Given (|V.29[) . the solution in its 
final form is as follows: 

D = ^ ^„ ^ ^^;^(2n - l)!(2n + l){n + 2) x '^^H^r.+2^ ^^-^^ (V.30) 

TT^ 2^" TT^ 



^9 = -V, = ^J?^ , (V.31) 
The macroscopic Einstein's equation is transformed into the equation for the inflation exponent 

3^2 ^ 3^^4 ^ ^^^2) 

OTT 

Because Ng is a macroscopic parameter, the solution under discussion can be directly relevant to the asymptotic 
future of the Universe. In this case, the number of gravitons and ghosts under the horizon of events and A-term in the 
equation (jV.32p should be considered as parameters, whose values were formed during the earlier stages of cosmological 
evolution. According to Zel'dovich [4l| . A-term is the total energy density of equilibrium vacuum subsystems of non- 
gravitational origin. The problem of the A-term formation is so complex that little has changed since the excellent 
review of Weinberg We are limited only to showing the order of magnitude of A ^ 3 • 10~'*^?i~^ GeV* allowed 
by observational data. (See also Appendix IXII ) 

Some possibilities of co-existence of graviton condensate and A-term will be discussed for A ^ 0, iVg > 0. (For 
other possibilities see Section fVl Bl ) The curvature of the De Sitter space for the asymptotical state of the Universe 
is calculated by means of the solution to the equation (|V.32p . It reads 



The energy density of vacuum in this state contains contributions of subsystems formed by all physical interactions 
including the gravitational one 

eiZ^ = ^Ht+A. (V.34) 

The relative input of graviton-ghost condensate into asymptotic energy density of the vacuum depends on parameters 
of the Universe. If the following inequality 

x^hKNn 

< 1 , (V.35) 



Btt' 

applies because of a small number of gravitons and ghosts, then the quantum-gravitational term is small and one 
must use the following solution 



Hlc.'-><K{l + ^Ng] . (V.36) 



1 . / x^hK 

-xK 1 H 

3 V 247r2 

If the inequality (|V.35P is satisfied because of a small A-term then the asymptotic state is mostly formed by the 
graviton-ghost condensate 
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It can be seen from (jV.33l) for A > 0, the number of gravitons and ghosts that can appear in the Universe is Umited 
by maximum value 



^sima.) = - 10122 . (V.38) 

In this hmiting case (jV.38p . the equipartition of the vacuum energy takes place between graviton-ghost and non- 
gravitational vacuum subsystems 




D. The Problem of Quantum— Gravity Phase Transitions 

Three exact solutions of the equations of quantum gravity (with no matter fields and in the absence of A-term) are, 
in our view, impressive illustrations of physical content of the theory. (Of course, we can not exclude the existence of 
other exact solutions). Before the integral 

* / 

(V.40) 

is calculated, three solutions are given by three different functionals Keg{a) that are displayed in the right-hand-side 
of the macroscopic Einstein equation 

l^(J)| r^i) (I) 
6a° a'^ a 

jj2 ^ In ^ (//) (V.41) 

3a6 + a2 (II) ' 

"0 

H = -rrr • (III) 

XhNg ^ ' 

If each of solutions (|V.4ip is considered as independent of the others, then one can note that (IV. 411 1) and (jV.41I II) are 
one-parameter solutions, meanwhile (jV.41I IIIl does not contain any free parameter. After multiplicative transforma- 
tions of scale factor a aao and time t tao/\Cg4\^^'^ in equations (|V.41I I) and (jV.41I II'). and time transformation 
t ^ t{KhNg/8n^y^^ in the equation (fVTITI lII). we get 

3a° a"^ 
= 1 . (Ill) 

Formulas (jV.42I I). (jV.42I II) and (jV.42I III) are special solutions of nonlinear system of equations allocated by special 
initial conditions. The first step is to determine what relationship they have to a general solution of equations (|V.7|) 
— (jV.10|) . corresponding to fairly arbitrary initial conditions. (Recall that the initial conditions are set by definitions 
of virtual gravitons and ghosts and structure of the state vector). Immediately note that we do not have an answer 
to this question in the form of strictly proven mathematical theorems. The mathematical problem is that we are 
dealing with a non-linear system, the number of degrees of freedom of which is infinite. This fact is reflected both 
in the operator formalism (an infinite number of modes, interacting through self-consistent field) and the BBGKY 
formalism (an infinite number of equations for moments of the spectral function). 

In examining the problem, using numerical experiments, the infinite system of equations (|V.7p — (jV.lOp is trans- 
formed into a finite system by breaking the chain (|V.9p . In the framework of this method, three questions are raise, 
(i) The choice of approximation of higher moment Wn+i through the lower ones Wi, W2, Wn; (ii) dependence 
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of solution asymptotes of initial conditions; (iii) dependence of the solution of the number of equations N in the 
chain. The third of these issues is trivial enough, in the sense that the response to it is produced by a mere repetition 
of numerical experiments with the sequential increase in N. In all the experiments that we conducted, there was 
convergence of observables for > 10 — 12 to some final functions 

ait), H{t), D{t), Wi{t) . (V.43) 

(The experiments are mostly completed for N = 20, but some of them held up to = 50.) 

We assumed that the three exact solutions are independent attractors of nonlinear system of equations. Under 
this assumption, the mathematical classification of attractors corresponds to the physical classification of possible 
asymptotic regimes of the Universe evolution. Breaking the chain (|V.9p is governed by a choice of asymptotics, and 
this is our proposed response to the first of the above issues. To truncate the chain, recurrence relations from exact 
solutions are used 

W^AT+l = ^-^Wn , (/) 

Wn+i^-^^Wn . {II) (V.44) 

— <-) 

We found by means of numerical experiments that all the exact solutions are stable with respect to small perturba- 
tions. Indeed experiments themselves have been limited to small variations of initial conditions at the vicinity of values 
fixing the exact solutions. In all cases, small perturbations quickly died out. It is necessary, of course, to remember 
that the statements about the stability are made on the basis of numerical experiments using approximations (|V.44p . 

Next, we conducted experiments with the initial conditions that have nothing to do with the conditions relevant 
to any of the three exact solutions. Nevertheless, the exact solution was used in selecting the appropriate method 
for truncating the chain. In all cases we saw a clear line between the ways truncating the chain (|V.44I I). (|V.44I II). 
(|V.44I III) and asymptotics of numerical solutions (|V.41I I). (|V.41I II). (|V.41I III). In fact, fixing the asymptotics by 
way of truncating the chain does not depend on initial conditions. However, all stages of evolution depend on the 
initial conditions, including the initial stage (which is natural), the intermediate stage of evolution and the nature of 
transition processes before the system reaches its asymptotic state. 

The intermediate stage of evolution in all cases was related to one of asymptotic (|V.44I I') or (|V.44I II). In Section 
IVIII Al it will be shown that such character at intermediate stages of evolution stems from the general properties of the 
equations of the theory. The type of solution xeg ^ a~^f(a), where /(a) is a slow function of scale factor, is formed 
as a result of excitation of graviton-ghost modes primarily on quasi-resonant frequencies. The nontrivial fact is that 
the regime of transition to the asymptotics determined by approximations (jV.44I I) or (|V.44I II) depends on initial 
conditions. In numerical experiments, we have witnessed either a smooth transition or a transition, accompanied by 
non-linear oscillations of physical quantities (|V.43[) . In the case of approximation (|V.44I III). the transition to the 
asymptotic De Sitter space always proceeds in the regime of non-linear oscillations. 

According to the results of numerical experiments, we came to the following conclusions. 

(i) Three exact solutions describe three different stable (at least, meta-stable) phases of graviton-ghost vacuum. The 
physical assumption that the Universe must be in one of these phases in the process of evolution, is formalized by the 
choice of the way of truncating the BBGKY chain. 

(ii) Arbitrary enough initial conditions correspond to the non-equilibrium state of vacuum in phases or I or II. 
These conditions are divided into two classes: a consistent and not consistent with the equilibrium phase, given by the 
approximation. 

(iii) // the physical nature of initial non-equilibrium phase matches the equilibrium phase chosen as asymptotic, 
the solution, starting with the intermediate stage, describes smooth relaxation of the graviton-ghost vacuum to an 
equilibrium state. If initial and asymptotic states do not match, a phase transition is initiated in the system. The 
signs of such a transition are nonlinear oscillations of physical quantities. 

(iv) The self-polarized graviton-ghost condensate in the De Sitter space can emerge only as a result of quantum 
gravity phase transition. 

The rationale for introducing of the notion of phases of graviton-ghost vacuum is the fact that three exact solutions 
match spaces with different symmetry. The solution (jV.44I III) describes 4-space of constant curvature, with the 
highest possible symmetry. Solution (|V.44I II) (in the version of appropriate unlimited expansion) describes 4-space, 
the geometry of which tends asymptotically to the geometry of the Milln space. Finally, the solution (jV.441 1') (in 
the version corresponding to oscillations) describes 3-geometry, which is translation-invariant along the axis of time. 
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Taking into account considerations of symmetry (see above), the term "phase of graviton-ghost vacuum" that we 
introduced seems mathematically and physically justified. 

Representations of phase transitions are, of course, only heuristic nature. In the one-loop quantum gravity, multi- 
particle correlations in the system of gravitons and ghosts are not taken into account. For this reason, in this theory it 
is impossible to define the order parameter that plays the role of the master parameter when choosing a phase state. 
Phase transitions that were discussed above, were actually initiated by disparity between the choice of the asymptotic 
state and set of the initial conditions. Of course, such operations are meaningful only within the suggestion that the 
effect of non- equilibrium phase transition will be contained in future theory. 

Staying on the heuristic level, we can use the exact solutions (jV.411 1'). (|V.41I II). (|V.41I III) to demonstrate in 
principle the possibility of the existence of equilibrium phase transitions. Let us consider the exact solutions as the 
various branches of a general solution. A rough phase transition model is the passage from one branch to another while 
maintaining continuity of scale factor and its first derivative. As can be seen from (IV.41I I'). (jV.41I II). (|V.41I III). these 
conditions provide the equality of energies of graviton-ghost systems on both sides of the transition point. The second 
derivative of the scale factor and vacuum pressure are discontinued (have a jump) at the point of transition. The 
microscopic theory also makes it possible to see that at the point of transition the internal structure of graviton-ghost 
substratum is changed (see Sections I VII IVIip . 

Consider consistently simplified models of all of the phase transitions. The condition of the sewing together solutions 
(fOn n and (fOn m has the form: 
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Uc is the value of the scale factor value at the fitting point, common to the two branches. The condition of the existence 
of transition between phases / and // is reduced to the inequality shown in (|V.45[) . As we know, in phase / gravitons 
dominate in quasi-resonant modes, and ghosts dominate in spatially uniform modes. Following the transition, in 
phase II quasi-resonant modes are dominated by ghosts, but spatially uniform modes are dominated by gravitons. 
Further, the condition of sewing together of solutions (|V.41I I') and (|V.41I III) reads 
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(V.46) 



It must be borne in mind that the energy density in phase / xsg^"^ is limited above and below. Therefore, in phase 
/// the number of gravitons under the horizon of events must lie in a certain interval, whose borders are defined by 
parameters of phase /. The phase transition looks like a "freezing" of the distance to the horizon and of the value of 
the physical wavelength of quasi-resonant modes. 

Finally, the third possible transition is illustrated by sewing together of solutions ([011 11) and (fOTj lll): 
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(V.47) 



In the most general case, the solution (jV.41I II) describes the birth of the Universe from the singularity and its further 
expansion with the acceleration. In this scenario for any preset value of energy density there is a corresponding point 
on the evolutionary path. Therefore, the transition from phase // to phase /// can occur at any point by choosing 
the appropriate value Ng. 



VI. EXACT SOLUTIONS: DYNAMICS OF OPERATORS AND STRUCTURE OF STATE VECTORS 

In this section, we get the exact solutions for field operators and expressions for the state vectors that correspond 
to exact analytical solutions of BBGKY chain (jV.41I I) and (|V.41I IIIV Microscopic studies of exact solutions allow 
greater detail to identify their physical content. Solutions (jV.41I I) and (jV.41I III') are formed as a result of certain 
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spectrally dependent correlations between graviton and ghost contributions to the observables. These are full graviton- 
ghost compensation of contributions of zero oscillations (one-loop finiteness); full compensation of contributions 
in all parts of the spectrum, except the region of quasi-resonant (QR) and spatially homogeneous (SH) modes; 
incomplete compensation of contributions of QR and SH modes with non-zero occupation numbers; correlations 
between excitation levels and graviton-ghost contents of QR and SH modes, and, finally, some correlations of phases 
in quantum superpositions of graviton and ghost state vectors. 

The physical nature of solution (|V.41I II) turned out to be unexpected and nontrivial. In Section IVIII it will be 
shown that mathematically this solution describes instanton condensate, which physically corresponds to the system 
of correlated fluctuations arising during tunneling of graviton-ghost medium between states with fixed difference of 
graviton and ghost numbers. We explain also that self-polarized graviton-ghost condensate in the De Sitter space 
also allows instanton interpretation. 



A. Condensate of Constant Conformal Wavelength 



Let us consider the solution (|V.41I I) for C^g — 0, C4g — fcg: 
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a — ao exp 



The graviton wave equation with the (jVI.ip background reads 
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(VI.l) 
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The equation for the ghosts looks similar. Fundamental solutions of equation (IVI.2P are degenerate hypergeometric 
functions. It is unnecessary to consider those solutions for all possible values of the parameter fc^. First of all, it is 
obvious that the macroscopic observables can be formed only by simplest hypergeometric functions. Values fc^ that 
are — Q (spatially uniform modes) and k^ ~ /cq (quasi-resonant modes) stand out. For all other modes there 
is a precise graviton-ghost compensation. The reason why it is a mathematically possible follows from the general 
formulas (pTMl) . (1111721) . (1111731) ^. 

Let us start with quasi-resonant modes. Exact solutions of the equation (|VI.1|) and similar equation for ghosts for 
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where Qktr, Pka and q^, are operators whose properties are defined in pil.53|) . pil.54[) . (|TTL50|; 

a 

da an 



F(a) = a'o 
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Formally, all modes except with = and = k^, look like "frozen" degrees of freedom, which are excluded from consideration by 
the model postulate. By virtue of the principle of uncertainty, postulates of this type are outside the formalism of quantum field theory. 
We want to emphasize that in the finite one-loop quantum gravity there is no need to "freeze" degrees of freedom not participating in 
the formation of particular exact solutions. Instead of mathematically incorrect operation of "freezing" , the formalism of the theory 
offers mathematically consistent operations of graviton-ghost compensations. 
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Note that one of fundamental solutions to equation (|VI.2p is the Hermite polynomial Hi{i]), which corresponds 
to positive eigenvalue k'^/k'^ — 1. In the reproduction of solutions (|V.41I I) at the microscopic level, this fact is 
crucial. We will show that the choice of a state vector, satisfying the condition of coherence leads to the fact that 
only this solution takes part in the formation of the observables. The second solution, containing a function F(a), is 
a mathematical structure that does not correspond to the exact solution to the BBGKY chain. 

Averaging of bilinear forms of operators (IVI.3|) and (IVI.4|) over the state vector of the general form leads to the 
following spectral function 

PFk = - 2{^gH\S+M^gH) = [Ak + B^F\a) + C^F{a)] In ^ . (vi.5) 

The constants appearing in (jVI.Sp are expressed through averaged quadratic forms of operators of generahzed coor- 
dinates and momentums: 

a 
a 

Following the transition to the ladder operators in formulas ()IIL54| and calculations, carried out similar to (|III.66P 
— (|III.74|) . we get 

Ak = 2K(3))(1 + Ci'^ cos^k) - 2K(<,^))(1 + Ct'"^ cosxfe) , 

Bk = 2(nfe(3))(l - Cl'^ cos^fc) - 2(nfc(,„))(l - (jf^^ cosxk) , (VI-7) 

Ck-0 . 

For sake of simplicity, in (IVI.7P average numbers of ghosts and anti-ghosts are assumed to be the same; {n}^[gh)) = 

{TT'kigh))- 

Let us go back to the expression (jVI.5p . Obviously, the spectral function (|VI.5p creates moments (jV.16p only if 
Bk — — 0. The condition Ck = is satisfied automatically as a consequence of isotropy of macroscopic state, i.e. 
because of independence of average occupation numbers of the direction of vector k. Bk = imposes the conditions on 
amplitudes and phases of quantum superpositions of state vectors with different occupation numbers. It is necessary 
to draw attention to the fundamental fact: the solution under discussion does not exist, if phases of superpositions 
are random. Indeed, averaging the expression (IVI.7P over phases, we see that condition i?k = is satisfied only if 
("fc(g)) = ("'fc(g/i))- The last equality automatically leads to Ak = 0, i.e. which eliminates the nontrivial solution. 

Thus, the condition of the existence of the solution under discussion is the coherence of the quantum state. It is 
easy to notice (see pil.741) '). that equality i?k = 0, as a condition of coherence, is satisfied for zero phase difference of 
states with the neighboring occupation numbers of gravitons and ghosts: 

Ci'^ COS = Cl"'' cos Xfc = 1 ^ Ci'^ = Cl'"^ = 1, cos ^fc = COS xfc = 1 . (VI.8) 

Taking into account (jVI.8|) . we get the following final expression (|VI.9[) for the spectral function of quasi-resonant 
gravitons and ghosts 

Wk^Wk = ^ - In ^ . (VI.9) 

In calculating moments, summation over wave numbers is replaced by integration. Account is taken of that the 
spectrum as the delta-form with respect to the modulus of k = |k|. Also a new parameter Ng is introduced where 
Ng is the difference of numbers of gravitons and ghosts in the unit volume of = J d^x = 1 in the 3-space, which 
is conformally similar to the 3-space of expanding Universe. Index "g" in designation of Ng parameter indicates the 
dominance of gravitons in quasi-resonant modes. In accordance with this definition, the following replacement is 
performed 

27r2 

("-fc(s)> - {nk(gh)) -TT^a^^^ ~ ^o) ' (VI.IO) 
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Results of calculating of moments are equated to the relevant expressions of (jV.lip and (jV.16p . which were obtained 
by exact solution of the BBGKY chain: 

oo 

(VI.ll) 

, , 1 / „ a' \ 128xhN„ko an 48fc^ an 



a 



In accordance with (jVLlip . there is a relation between parameters fco, aq and Ng that appear in the microscopic 
solution 



Recall that in the solution under discussion, the Universe was born in singularity, expands to a state with a maximum 
scale factor amax = clq, and then is again compressed to the singularity. In this scenario, value oq can be defined as 
the size of the Universe, accessible for observation in the end stage of expansion. As can be seen from (jV1.12p . if oq 
is a macroscopic value, the difference in numbers gravitons and ghosts Ng ^ 1 is also a macroscopic value. 

Contributions of SH modes to the expressions for the moments are shown in (jV.lip . and the relation between the 
parameters C2g and C^g is shown in (|V.15[) . As a part of the microscopic approach, the construction of exact solutions 
for these modes is performed by the method of transaction to the limit, described at the end of Section flVBI The 
parameter of spatially homogeneous condensate is introduced similarly to (|VI.10[) : 

{noigh))il + Co""^ cos(^o) - ("o(<,))(l + Co"^ cosxo) ^ -j^NghS{k - kq) , Ko^O. (VI.13) 

The index "(7/1" in Ngh > indicates the dominance of ghosts over the gravitons in the spatially homogeneous 
condensate. The moments are: 

afa'' afa^ 

Definitions of parameters ki and ai are given in pV.2ip . The energy density and pressure of the system of QR and 
SH modes are given by (|VI.11|) and (|VI.14|) : 



S^chkoNg flQ 2>ch f koNg kiNnh\ S^chkoNn ao 

^£g = — 2 ? — ^ — r — 2^ T~ = — r^T^ 1^ ~ 

aia-' a a"^ \ a?, ai I aia'^ a 



2xh 1 


^koNg 


kiNgh 




\ al 


al 


2>ch , 


fkoNg 


kiNgh 


3^ ' 




al 



8>chkoNg J ao 2>ch f k^Ng kiNgh\ SxhkoNg ^ ao 
3ana^ ea 3a^ \ a,n a? J 30^0^ ea 



(VI. 15) 



In formulas (|V1.15p . the terms in brackets are eliminated by the condition (|V.15p . which is rewritten in terms of 
macroscopic parameters 

koNn kiN„h 

J)g ^ J_gh_ VI.16 

The solution (|VI.15p . (|VI.16p describes a quantum coherent condensate of quasi-resonant modes with graviton domi- 
nance, parameters of which are consistent with that of spatially homogeneous condensate with the ghost dominance. 

B. Condensate of Constant Physical Wavelength 

The De Sitter solution for plane isotropic Universe reads 

a^aoe"* ^ --^ , H ^ const . (VI. 17) 

For the background (jVl.lTp . the gravitons and ghost equations and their solutions read 



- + fc'^k. = , = y ^ [ck./(x) + e\_J*{x)\ , (VI.18) 
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K - -K + fc'^k = , ,9k = [ak/(x) + b\r{x)] , (VI.19) 

r] a y k 

where 

fix) = (l - ^) , ^^kr^- 

Ladder operators in (jVI.18|) . (jVI.lQp . have the standard property of pil.47p . (|III.50|) . which aUow their use of in 
constructing build basic vectors for the Fock space from which the general state vectors are constructed. 

The self-consistent dynamics of gravitons and ghosts in the De Sitter space are not trivial in the sense that the 
averaged bilinear forms of operators (jVLlSp . (jVI.19p which are explicitly and essentially depending on time, must 
lead to time-independent macroscopic observables. It must be emphasized, that the existence of such, at first glance 
unlikely solution, is guaranteed by the existence of the solution for the BBGKY chain. The key to the solution lies in 
the structure of the state vectors of gravitons and ghosts. 

Substitution of operator functions (jVLlSp . (|VI.19|) into the general expression for the moments (|V.6|) yields: 

oo 

^" = ^^'"^'/ «'^^'"+'{c^k(.,a.e)l/(a^)P + Ukicr)[n^)? + C/k(a„„) [/(a^)]' | , (VI.20) 


where 

= Uki^ave) - ^(*s|c+.Ck,|*9> - 1 4 «k | *g/0 - (^g/J &k | ^g/.) ; (VI.21) 



1 " (VI.22) 

Uk = C/k(arm) - - ^ (*g | C_k-<.Ck<. | - ('J'g/. | fe-kOk | 'J'g/O = C/k(cr) ■ 

a 

Here C/k(u)at)e) is the spectral parameter of quantum waves, which become real gravitons if ^77 1; C/k(cr)j ^^k(arm) are 
the spectral parameters of quantum fluctuations that emerge in the processes of graviton (and ghost) creation from 
the vacuum and graviton (and ghost) annihilation to the vacuum. 

Obviously, at the first stage of calculations we assume that the averaging in (jVI.2ip , (|VI.22|1 is conducted over the 
state v ectors of the general form (|III.48p . (|III.5ip . This allows us to go to formulas (jIII.67p . (|III.68P or (|III.72p ^ 
piL74p . Then it is necessary to take into account that the moments Wn must not depend on time, and that they also 
should be free of divergences. When analyzing the conditions for these demands, the specific form of the expression 
()VI.20p plays an important part. The measure of integration and the dependence of field operators on the wave 
number and time can be represented in the terms of the variable x = krj. A separate (additional) dependence on the 
wave number can be connected with the structure of spectral parameters. After substitution of the variable k = x/rj 
in the equation (|VI.2ip . it is seen that the first term in (jVI.20p is time- independent only if C/k(toat)e) is independent 
on the wave number. This means that the graviton and ghost spectra must be fiat. However, with the flat spectrum 
there is danger of divergences: if /7k(M,ai;e) = const (k) 7^ 0, then the first integral in (jVI.20p does not exist, because 
|/(x)p 1 with x ^ 00. 

The divergences can be avoided only with exact compensation of contributions from gravitons and ghosts to the 
spectral parameter C/k(tDot)e)- Let us point out, that in that case we are not talking about zero oscillations but 
about the contributions from the states with non-zero occupation numbers. The compensation condition leading to 

C^k(i«ai;e) = is: 

|C„.J = = . (VI.23) 

The result (|VI.23p has a simple physical interpretation. The quantum waves of gravitons and ghosts with the equation 
of state which differs from p — —e can not be carriers of energy in the De Sitter space with the self-consistent geometry. 
The total energy of quantized waves is equal to zero due to exactly the same number of gravitons and ghosts in all 
regions of the spectrum: 



(«kcri) + (?lkcr2> = ('^k) + (?^k) 



(VL24) 
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With equal polarizations of gravitons and the equality of numbers of ghosts and anti-ghosts, it follows from (jVI.24p 
that ('n.k(g)) = (nkigh))- Exact equality of the average number of gravitons and ghost is a characteristic feature of 
the De Sitter space with the self-consistent geometry. Let us mention that for the solution discussed in the previous 
section lVl Al that equality is absent in principle. It means that different solutions have different microscopic structures 
of the graviton-ghost condensate. 

Based on the reasoning analogous to the one described above, spectrum parameters ?7k(cr )i C4[;(ann) also must not 
depend on the wave vector k. However, the corresponding integrals in the second and third terms of (|VI.20p are 
not divergent. The absence of divergences is due to the fact that with a; ^ oo the integration is taken over the fast 
oscillating functions ~ e^^". To calculate these integrals, they should be additionally defined as follows: 



oc 

^limJd..--e-(C--)^ = T(-ir^^, 



oo 

2i hm / dxx^'^e-^'^-^'^ = . 



(VI.25) 



At every instant of time, the procedure of re-definitions of integrals (jVI.25P selects the contributions from virtual 
gravitons and ghosts with a characteristic wavelength (|V.28p and eliminate the contributions of all other graviton- 
ghost modes. This redefining procedure provides the existence of recursive relations (jV.27p in the exact solution of 
the BBGK Y chain . 

Thus, in (jVI.20[) we have a flat spectrum of gravitons and ghosts, f/k(iijat)e) = Oj t^k(cr) = ^k(ann) ^ ^ ^ const[k). 
The expression for the spectral parameter takes the form: 



U 



\ n / \ n / \ n / 



(VI.26) 



\Cn\ = \An\ = \Bn\ = , 



where Vn is a normalized statistical distribution. The average value of the number of gravitons and ghosts, having 
the wavelength in the vicinity of characteristic values (jV.28p . are calculated by the formula 



ngh) = [n) 



J2nV{n) . (VI.27) 



n=0 



Using the Poisson distribution in (|VI.26p . (|VI.27p . the values of integrals (|VI.25P and the formulas (|III.73|1 . (|III.74|) . 
we get the moments 

D = -}^^H^ , Wr. = ^4£^(2" - 1)!(2« + + 2) X ^ ^ ^ ^ (Vl.28) 

where 



Ng = {n) (Cg cos (p - Cgh cos x) ■ (VI.29) 

Zero moment Wg, which has an infrared logarithmic singularity, is not contained in the expressions for the macroscopic 
observables, and for that reason, is not calculated. In the equation for Wq, the functions are differentiated in the 
integrand and the derivatives are combined in accordance with the definition D = Wq + 377 Wq- At the last step the 
integrals that are calculated, already posses no singularities. 

Averaging of the parameter (jVI.29[) over the phases yields Ng — 0. Therefore the solution under discussion does not 
exist if the superposition of the phases are random. The coherence of the quantum ensemble, i.e. the correlation of 
phases in the quantum superposition of the basic vectors, corresponding to the different occupation numbers, points 
to the fact that the medium is in the graviton-ghost condensate state. The gravitons are dominant in the condensate 
if A'g > 0, and the ghosts are dominant if iVg < 0. 

The duality of the condensate and the indeterminate sign of the A-term create different evolutional scenarios. Of 
course, all these scenarios are present in the expression (|V.33p . which is obtained as a solution of the macroscopic 
Einstein equation (|V.32[) . In addition to the scenarios described in the Section IV C[ we will show the possibihty 



56 



of strong renormalization of energy of non-gravitational vacuum subsystems by the energy of the graviton-ghost 
condensate^. 

We have in mind a situation, in which the modulus of A-term exceeds the density of vacuum energy in the asymptotic 
state of the Universe by many orders of magnitude: 

1^1 -""^"-AA»1, (VI.30) 



where A/" is a huge macroscopic number. From (IV.33|) it follows that the effect of strong renormalization takes place if 



Let us mention that the strong renormalization of the positive A-term is provided by a condensate in which the ghosts 
are dominant, and for the negative A-term — by a condensate for which the gravitons are dominant. 

For clarity and for the evaluations let us introduce the Plank scale Mpi = {SttH/ >cy^'^ = 1.22 • 10^^ GeV, the scale 
of A-term Ma = (fi-'^|A|)-'^/'', and the scale of the density of Dark Energy in the asymptotical state of the Universe, 

Mde — (fi-^elTcV^''- We discuss the case when Mjje < -^^A- 

If non-gravitational contributions to A-term are self-compensating, then a realistic estimate of the Ma -scale can be 
based on the Zeldovich remark [41]. According to (4l| . non-gravitational A-term is formed by gravitational exchange 
interaction of quantum fluctuations on the energy scale of hadrons. In terms of contemporary understanding of 
hadron's vacuum, the focus should be on non-perturbative fluctuations of quark and gluon fields, forming a quark- 
gluon condensate (see Appendix IXip . In this case, A-term is expressed only through the minimum and maximum 
scales of particle physics which are the QCD scale Mqcd — 215 MeV and Planck scale Mpi = 1.22 • 10^^ GeV: 

/If 6 

;j3|A| = Ml = -ff^ ~ 10-42 GeV^ . (VI.32) 

In terms of these scales, it is turns out that a large number oi Af ^ A-lj^/M^^^ ^ 10^, which is defined in (jVI.3011 . 
can be obtained by the huge number of jiVgl^/^, for the same number of orders of magnitude greater than the ratio 
{Mpi/M{Cf' . Indeed, choosing A/", we find the value of |A^g|, which determines the ratio of vacuum energy density to 
the true cosmological constant in the asymptotic state: 



The vacuum energy density of asymptotical state is calculated as follows 



si^)^h-^Mj.,Ml^^^ . (VI.34) 

Thus, the macroscopic effect of quantum gravity — the condensation of gravitons and ghosts into the state with a 
certain wavelength of the order of the horizon scale — plays a significant role in the formation of the asymptotic values 
of energy density of cosmological vacuum. The current theory explains how the strong renormalization of the vacuum 
energy occurs, but, unfortunately, it does not explain why this happens and why the quantitative characteristics of 
the phenomenon are those that are observed in the modern Universe. Of the general considerations one can suggest 
that the coherent graviton-ghost condensate occurs in the quantum-gravitational phase transition (see Section [VPp . 
and the answers to questions should be sought in the light of the circumstances. 



Mechanisms that are able to drive the cosmological constant to zero have been discussed for decades (see [42l |43|| for a review) . Any 
particular scenarios were considered in [44l l45l 143. 147l l4 8ll . Renormalizations of cosmological constant by gravitons in the framework of 
one-loop quantum gravity were also considered in (20l l2lL |22| . l25j where the effect of reconstruction of zero oscillations of gravitational 
field in the self— consistent De Sitter space, i.e. effect of conformal anomalies was discussed. Conformal anomalies that arise due to 
regularization and renormalization procedures do not apply to this work (see also Appendix IXTTCj 
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VII. GRAVITONS AND GHOSTS AS INSTANTONS 



A. Self— Consistent Theory of Gravitons in Imaginary Time 



1. Invariance of Equations of the Theory with Respect to Wick Rotation of Time Axis 



As has been repeatedly pointed out, the complete system of equations of the theory consists of the BBGKY chain 
(IV. 7[) — (jV.9[) and macroscopic Einstein's equations (|V.10|) . On the basis of common mathematical considerations, it 
can be expected that solutions to these equations covers every possible self-consistent states of quantum subsystem 
of gravitons and ghosts and the classical subsystem of macroscopic geometry as well. In examining the model that 
operates with the pure gravity (no matter fields and A-term), one can identify the following unique property of the 
theory. Equations of the theory f K ?| ) — iV.l 0\) are invariant with respect to the Wick time axis rotation, conducted 
jointly with the multiplicative transformation of moments of the spectral function: 

t^iT, -iH, D -V, Wn {-irWn ■ (VII.l) 

Rules of transformation of time derivatives are obtained from (j VII. II) 

H -> -n , H in" , D^iv' , 

(VIL2) 

Wn ^ -*(-l)"W; , Wn -> -(-1)"W;- , Wn^ ^i-^r^n ■ 



In (|VII.2p and further on we use the notation JF' = dJ-fdr. The statement about the invariance of the theory can 
proved by direct calculations. As a matter of fact, transformations of quantities that appear in (IV. 7p — (|V.10|) by 
the use of the rules (jVII.ip and (|VII.2p lead to the BBGKY chain with imaginary time 

V' + mv + 4>v; + imwi = o , (vii.3) 



W;" +3(2n + 3)7^W;' +3 



[An^ + 12n + 6) Ti^ + {2n + l)rL 



(VII.4) 



'2n 



2 (2n2 + 9n + 9) + 6(n + 2)nii: + U" W„ + 4W;+i + 8(n + 2)HW„+i =0, n = 1, oo , 



and to macroscopic Einstein's equations with imaginary time 



1 1 
16 6 

16 4 



(VIL5) 



It is easy to see that for A = equations (|V.7p — (IV.lOp identically coincide with (|VII.3p — (|VII.5p after some trivial 
renaming. 

The invariance of the theory with respect to the Wick rotation of the time axis leads to the nontrivial consequence. 
Having only self-consistent solution of the BBGKY chain and macroscopic Einstein's equations, we can not say 
whether this solution is in real or imaginary time. Nevertheless, having a concrete solution of BBGKY chain, we can 
view the status of time during further study. To do so, it is necessary to explore the opportunity to obtain the same 
solution at the level of operator functions and state vectors. If this opportunity exists, the appropriate self-consistent 
solution of BBGKY chain and macroscopic Einstein's equations is recognized as existing in real time. In the previous 
Section lVll we showed that two exact solutions (|V.41I I) and (|V.41I III) really exist at the level of operators and vectors, 
and thus have a physical interpretation of standard notions of quantum theory. 

The problem is: What a physical reality reflects the existence of solutions to the equations |y. 7| ) — W.l Q|) ( or that 
the same thing, i f VII. 3\) — if VT/.5)) ). not reproducible in real time at the level of operators and vectors? The existence 
of the problem is explicitly demonstrated by the example of exact solutions (IV.41I II). Assume that this solution for 
Csg ~ 0, C4g ~ — fcg < 0, C g2 — 3/^0 /4 cxlsts lu rcal time: 



H' 



Oq exp 



(VII.6) 
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The wave equation for gravitons with the (jVII.6|) background reads 

i'L + klvi>L + fc'^k. = . (VIL7) 

The equation for the ghosts looks similar. Equation (|VII.7p differs from (|VI.2p just in the sign of coefficient before 
the first derivative. However, this difference is crucial: if k'^ /kg > it is impossible to allocate the finite Hermit Hi{rj) 
polynomial from degenerate hypergeometric functions that correspond to solutions of the equation (|VH.7|) . We have 
been left with the infinite series only. These series and integrals over spectrum of products of these series can not be 
made consistent with the simple mathematical structure of the exact solutions (jV. 411 11). For this reason the solution 
(|V. 411 11). as the functional of scale factor is not relevant to solving operator equations in real time. 



2. Imaginary Time Formalism 

As is known, the imaginary time formalism is used in non-relativistic Quantum Mechanics (QM) (examples see, e.g., 
in book [i^), in the instanton theory of Quantum Chromodynamics (QCD) [HO, IMl 113, IHj lEj] and in the axiomatic 
quantum field th eory (AQFT) (See Chapter 9 in the monograph [28|). The instanton physics in Quantum Cosmology 
was discussed in [53, [56|] . 

In QM and QCD the imaginary time formalism is a tool for the study of tunnelling, uniting classic independent 
states that are degenerate in energy, in a single quantum state. In AQFT, the Schwinger functions are defined in the 
four-dimensional Euclidian space — Euclid analogues of Wightman functions defined over the Minkowski space. It 
is believed that using properties of Euclid-Schwinger functions after their analytical continuation to the Minkowski 
space, one can reconstruct the properties of Wightman functions, and thereby restore the physical meaning of the 
appropriate model of quantum field theory. 

All prerequisites for the use of the formalism of imaginary time in the QM and QCD on the one hand, and in AQFT, 
on the other hand, are united in the self-consistent theory of gravitons. Immediately, however, the specifics of the 
graviton theory under discussion should be noted. Macroscopic space-time in self-consistent theory of gravitons, unlike 
the space-time in the QM, QCD and AQFT, is a classical dynamic subsystem, which actually evolved in real time. 
If in QCD and AQFT Wick's turn is used to examine the significant properties of quantum system expressed in the 
probabilities of quantum processes, then in relation to the deterministic evolution of classical macroscopic subsystem 
this turn makes no sense. Therefore, after solving equations of the theory in imaginary time, we are obliged to apply 
(to the solution obtained) the operation of analytic continuation of the space for the positive signature to the space 
of negative signature. It is clear from the outset that the operation is not reduced to the opposite Wick turn, but is 
an independent postulate of the theory. 

Before discussing the physical content of the theory, let us define its formal mathematical scheme. The theory is 
formulated in the space with metric 

ds^ = -dr^ - a^{T){dx^ + dy^ + dz^) . (VH.S) 

Note that in our theory, that is suppose to do with cosmological applications (as opposed to QCD and AQFT), one of 
the coordinates is singled out simply because the scale factor depends on it. This means that in the classical sector of 
the theory time r, despite the fact that it is imaginary, is singled out in comparison with the 3-spatial coordinates. In 
the quantum sector the r coordinate also has a special status. Operators of graviton and ghost fields with nontrivial 
commutation properties are defined over the space (jVll.Sp . Symmetry properties of space (jVII.Sp allow us to define 
the Fourier images of the operators by coordinates x, y, z, and to formulate the canonical commutation relations in 
terms of derivatives of operators with respect to the imaginary time r: 

= -iMkk"5aa' ■ (vn.9) 



a 



kcr 



dr 



k'cr' 



a" 



-7^ , 1?k 

dr 



-ihS] 



kk' 



a" 

4:>C 



dT 



-ihS] 



kk' 



(VH.IO) 



Note that (|V11.9p . (jVII.10|) are introduced by the newly independent postulate of the theory, and not derived from 
standard commutation relations pil.37p . pil.4ip by conversion oi t ^ it. (Such a conversion would lead to the 
disappearance of the imaginary unit from the right hand sides of the commutation relations.) Thus, the imaginary time 
formalism can not be regarded simply as another way to describe the graviton and ghost fields, i.e. as a mathematically 
equivalent way for real time description. In this formalism the new specific class of quantum phenomena is studied. 
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The system of self-consistent equations is produced by variations of action, as defined in 4-space with a positive 
signature: 



dr < 3 



a d^a a dN da a f da 




(VII. 11) 



Note that the full derivative with respect to the imaginary time is not excluded from Lagrangian. In (jVILlip the 
integrand contains the density of invariant ^/gR- The Lagrange multiplier N after the completion of the variation 
procedure is assumed to be equal to unity. The system of equations corresponding to the action (jVII.lip can also be 
obtained from the system of equations in real time by conversion of i — > it. Quantum equations of motion for field 
operators in the imaginary time read 

■ o + 37Y— 5 ^-V^kcr = , (VII. 12) 

dr-' dr 

+ m—^ - ^i?k - , (VII.13) 



dr^ dr 
where Ti = a /a. 

Equations (|VII.12|1 . (|VII.13j) differ from (|III.30|) . (|III.3ip by only replacement of fc^ ^ -A:^. At the level of analytic 
properties of solutions of the equations this difference, of course, is crucial. However, formal transformations, not 
dependent on the properties of analytic solutions to equations (|III.30p . (jIII.3ip and (|VII.12p . (|VII.13p . look quite 
similar. Therefore, all operations to construct the equation for the spectral function in imaginary time (analogue to 
equation (|V.4p ) and the subsequent construction of BBGKY chain coincide with that described in Section fV Al with 
the replacement of fc^ — > —k'^. Replacing — )■ —k^ changes the definition of moments only parametrically: instead 
of (fVTBl) we get 

W"^Ef4^) (E(*9l^k.V^k.|*s)-2(M'g;,|^+4|*s'0 ) , n = 0, 1, 2, ...,oo, 

k V / V . / (VII. 14) 

dr^ dr 

Further actions lead obviously to the BBGKY chain (|VII.3p . (|VII.4p and to the macroscopic Einstein equations (jVII.Sp . 
To solve equations (|VII.12p and (|VII.13p . we will be using only the real linear-independent basis 

■0kCT = V4><ft {Q]^agk + Pkahkj , "i^k = V^>ih (qkgk + Pkhk) 



9khk - hkQk - 3 



^ (VII. 15) 



As will be seen below, one of the basic solutions satisfies the known definition of instanton: an instanton is a solution 
to the classical equation, which is localized in the imaginary time and corresponds to the finite action in the 4-space 
with a positive signature. We will call the operator functions (|VII.15p the quantum instanton fields of gravitons 
and ghosts. Operator constants of integration Qko-, Pka and qk, Pk satisfy commutation relations (jIII.55p . Ladder 
operators are imposed by the the equations pil.54p and then used in the procedure for constructing the state vectors 
over the basis of occupation numbers. State vectors of the general form in graviton and ghost sectors are already 
familiar structure pil.48p and pil.5ip . Only the interpretation of occupation numbers is changed: now it is number 
of instantons Uka, nk, nk of graviton, ghost and anti-ghost types, respectively. 
Direct calculation of the moments of the spectral function leads to the expression: 

W„=4xft(-l)"^f-2j (Akgl + Bkhl) , (VIL16) 

k 

where 

Ak = ^(*g|QLQk.|«'g) - 2(vl/,,|g+9k|*3h) = 

(VII. 17) 

= 2(nfe(g))(l + Cfe^^ cos ipk) - 2{nk{gh))il + Ci^''^ cosxfc) , 
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(VII. 18) 

= 2(nfe(g))(l - Cfe^^ cos ipk) - 2(nfc(gft.))(l - C^^''' cosxfc) • 

The term containing products of basis functions guhk is eliminated from (|VII.16p by the condition of homogeneity of 
3-space. In (jVII.17[) and (jVII.18[) average values of numbers of instantons of ghost and anti-ghost types are assumed 
to be equal: ('^fc(g/i)) = {nk(gh))- One needs to pay attention to the multiplier (—1)" in (|VII.16P : the alternating sign 
of moments is a common symptom of instanton nature of the spectral function. 

Instanton equations of motion (|VII.12p . (|VII.13p arc of the hyperbolic type. This fact determines the form of 
asymptotics of basis function for fc|^| 3> 1 where ^ = J dr/a is conformal imaginary time. One of basis functions 
is localized in the imaginary time and the other is increasing without limit with the increasing of modulus of the 
imaginary time 

9k ~ ^ , hk ^ ^ , m » 1 . (VII.19) 
av Zk ay Ik 



In this situation, it is necessary to differentiate between stable and unstable instanton configurations. We call a 
configuration stable, if moments of the spectral function are formed by localized basis functions only. Without 
limiting generality, we assigned hk to the class of increasing functions. It is easy to see that the condition of stability 
Bk = that eliminates contributions of hk from (|VII.16|) is reduced to the condition of quantum coherence of instanton 
condensate: 

Ci^)cos^fe=C["'^cosxfe = l ^ Ci'^=d'"' = l, cosv^fc=cosxfe = l . (VII.20) 
Expressions for the moments are simplified and read 



/ l,2\ » 

4^ft(-irE^k(^^j gi- (VII.21) 



Exact solutions, with the stable instanton configurations, are described in the following Sections IVII Bl and IVII CI In 
principle, for a limited imaginary time interval, there might be unstable configurations, but in the present work such 
configurations are not discussed. (The example of the unstable instanton configuration see in [st}.) 

Note that moments (IVIL21I) can be obtained within the classical theory, limited, as generally accepted, to the 
solutions localized in imaginary time. In doing so, Ak acts as a constant of integration of classical equation. 

The above approach is the quantum theory of instantons in imaginary time. Here are present all the elements of 
quantum theory: operator nature of instanton field; quantization on the canonical commutation relations; basic vectors 
in the representation of instanton occupation numbers; state vectors of physical states in the form of superposition of 
basic vectors. With the quantum approach, a significant feature of instantons is displayed, which clearly is not visible 
in the classical theory. It is the nature of instanton stable configurations as coherent quantum condensates. 

Construction of the formalism of the theory is completed by developing a procedure to transfer the results of the 
study of instantons to real time. It is clear that this procedure is required to match the theory with the experimental 
data, i.e. to explain the past and predict the future of the Universe. As already noted, the procedure of transition to 
real time is not an inverse Wick rotation. This is particularly evident in the quantum theory: in (|VII.9p . (jVILlOp the 
reverse Wick turn leads to the commutation relations for non-Hermitian operators, which can not be used to describe 
the graviton field. 

The procedure for the transition to real time has the status of an independent theory postulates. We will formulate 
this postulate as follows. 

(i) Results of solutions of quantum equations of motion (|VII.12p , (jVILlSp , together with the macroscopic Einstein's 
equations (|VII.5p after calculating of the moments (that is, after averaging over the instanton state vector) should be 
represented in the functional form 

v = v{a,n,n',...) , Wn^Wnia,n,n' ,...) . (vii.22) 

(ii) It is postulated that functional dependence of the moments of the spectral function on functions describing the 
macroscopic geometry must be identical in the real and imaginary time. Thus, at the level of the moments of the 
spectral function, the transition to the real time is reduced to a change of notation 

V{a,n,n,...) D{a,H,H,...) , 

(VII.23) 

Wn{a,n,n\...) ^ Wn{a,H,H,...) . 



61 



(iii) Moments -D(a, ff, i?, ...) and Wi{a, H, H, ...) obtained by operations (|VII.23p . are substituted to right hand side 
of macroscopic Einstein equations that are considered now as equations in real time. FormaUy this means that the 
transition to the real time in the left hand side of equations (jVII.SP is reduced to changing of the following notations 

n H, ^ . (VII.24) 

Thus, the acceptance of postulates (jVII.22|) — (jVII.24ll is equivalent to the suggestion that in real time the self- 
consistent evolution of classic geometry and quantum instanton system is described by the following equations 

H = -^Dia, H, H, ...) - Vi(a, H, H, ...) , 

16 6 , ^ 

(VII.25) 

= ^D{a, H, H, ...) + ]wi{a, H, H, ...) , 
16 4 

under the condition that the form of functionals in right hand sides of (jVII.25|l is established by microscopic calculations 
in imaginary time. It is obvious also that in the framework of these postulates any solution of equations consisting of 
BBGKY chain and macroscopic Einstein equations (obtained without use of microscopic theory) can be considered as 
the solution in real time. 



3. Physics of Imaginary Time 



Mathematical and physical motivation to look for the formalism of imaginary time comes from the fact that there 
are degenerate states separated by the classical impenetrable barrier. In non-relativistic quantum mechanics the 
barriers are considered, that have been formed by classical force fields and for that reason they have the obvious 
interpretation. It is well known, that the calculation of quantum tunnelling across the classical impenetrable barrier 
can be carried out in the following order: (1) the solution of classical equation of motion inside the barrier area is 
obtained with imaginary time; (2) from the solution obtained for the tunnelling particle, one calculates the action 
S for the imaginary time; (3) the tunnelling probability, coinciding with the result of the solution for Schrodinger 
equation in the quasi-classical approximation, is equal w = e~^ . Obviously, the sequence described bears a formal 
character and cannot be interpreted operationally. Nevertheless, a strong argument toward the use of the formalism 
of imaginary time in the quantum mechanics is the agreement between the calculations and experimental data for the 
tunnelling micro-particles. 

A new class of phenomena arises in the cases when tunnelling processes form a macroscopic quantum state. The 
Josephson effect is a characteristic example: fluctuations of the electromagnetic field arise when a superconductive 
condensate is tunnelling across the classically impenetrable non-conducting barrier. Here, the tunnelling can he 
formally described as a process developing in imaginary time, but the fluctuations arise and exist in the real space- 
time. Experimental data show that regardless of the description, the tunnelling process forms a physical subsystem 
in the real space-time, with perfectly real energy-momentum. 

In Quantum Chromodynamics (QCD) physically similar phenomena are studied by similar methods. The vacuum 
degeneration is an internal property of QCD: different classical vacuums of gluon field are not topologically equivalent. 
In the framework of the classic dynamics any transitions between different vacuums are impossible. In that sense 
the topological non-equivalence plays role of the classical impenetrable barrier. There is an heuristic hypothesis 
in quantum theory — that the probability of tunnelling transition between different vacuums can be calculated as 
w — e^^ , where S is the action of the classical instanton. The instanton is defined as a solution of gluon-dynamic 
equations localized in the Euclidian space-time connecting configurations with different topologies. As in the case 
of Josephson Effect, it is assumed that the tunnelling processes between topologically non-equivalent vacuums are 
accompanied by generation of non-perturbative fluctuations of gluon and quark fields in real space-time. Let us notice 
that in QCD the instanton solutions, analytically continued into real space-time, are used to evaluate the amplitude of 
fluctuations. The fluctuations in real space-time are considered as a quark-gluon condensate (QGC). The existence 
of QGC with different topological structure in "off-adrons" and "in-adrons" vacuums, is confirmed by comparison of 
theoretical predictions with experimental data. One of remarkable facts is that the carrier of approximately the half 
of nucleon mass is in fact the energy of the reconstructed QGC. 

Now let us go back to the self-consistence theory of gravitons. In that theory, due to its one-loop finiteness, all 
observables are formed by the difference between graviton and ghost contributions. That fact is obvious both from 
the general expressions for the observables (see (jIII.72p . pil.73[) ). and from the exact and approximate solutions 
(described in the previous sections) as well. The same final differences of contributions may correspond to the totally 
different graviton and ghost contributions themselves. All quantum states are degenerated with respect to mutually 
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consistent transformations of gravitons and ghosts occupation numbers, but providing unchanged values of observable 
quantities. Thus the multitude of state vectors of the general form, averaging over which leads to the same values 
of spectral function, is a direct consequence of the internal mathematical structure of the self-consistent theory of 
gravitons, satisfying the one-loop finiteness condition. 

In that situation, it is very natural to introduce a hypothesis about the tunnelling of the graviton-ghost system 
between quantum states corresponding to the same values of macroscopic observables. By the analogy with the 
effects described above, one may suggest that 1) the tunnelling processes unite degenerate quantum states into a 
single quantum state; 2) tunnelling is accompanied by creation of specific quantum fluctuations of graviton and ghost 
fields in real space-time. With regard to the mathematical method used to describe these phenomena, today we may 
use only those methods that have been tested in adjacent brunches of quantum theory. It is easy to see that this 
program has been realized in Sections IVII A 11 IVII A 21 We solve the equations of the theory for imaginary time, but 
the amplitude of the arising fluctuations we evaluate by the analytical continuation (|VII.22[) — (jVII.24[) , analoguos 
to the ones used in QCD. The specific of our theory lie in the fact that at the final step of calculations we use the 
classical Einstein equations (IVII.25P describing the evolution of the macroscopic space in real time. The possibility 
of using these equations is determined by the action (jVII.ll|) . which, when calculated by means of the instanton 
solutions and averaged over the state vector of instantons, is identically equal zero. As a matter of fact, after using 
instanton equations (|VII.12|) . (|VII.13p and averaging, the action (jVII.lip is reduced to the form: 



= ^J dra^ 



3{H +H^]+ —V 



1 

16' 



(VII.26) 



The integrand in (IVII.26|) is equal zero in the Einstein equations with imaginary time (|VII.5[) . The fact that w — 
exp(— (^'IS'l^E')) = 1 means that the macroscopic evolution of the Universe is determined. That feature allows the use 
of equations (jVII.25|) . after the moments are analytically continued into the real time. 



B. Instanton Condensate in the De Sitter Space 



Among exact solutions of the one-loop quantum gravity, a special status is given to De Sitter space if the space 
curvature of this space is self-consistent with the quantum state of gravitons and ghosts. In Section FVI Bl it was shown 
that in the self-consistent solution, gravitons and ghosts can be interpreted as quantum wave fields in real space- 
time. Nevertheless, it should be mentioned, that the alternating sign of the moments (|VI.28p points to a possibility of 
instanton interpretation of that solution. Methods described in Sections IVII A 11 IVII A 21 when applied to De Sitter 
space, show that such interpretation is really possible. 

We will work with the imaginary conformal time £, — J dr/a. The cosmological solution is: 

a = aoe"^ = — ^ , -c5o < ^ < . (VII.27) 

At the level of the BBGKY chain, due to the fact that the theory is invariant with respect to the Wick rotation, the 
calculations performed to get the solutions coincide with the those described in Section IV CI At the microscopic level 
we use the exact solutions (|VII.12I) . (|VII.13p with the background (|VII.27|) : 



where 



V^k. = - \/ ^ [Qkaff (x) + Pk./i(x)] , = - [qugix) + Pk/i(x)] , (VII.28) 

a y k a \ k 



g{x) = J^l - , h{x) = + ^) a; = fc^ < . 

The expressions for the moments of the spectral function are reduced to the form: 



Wn = (-l)"+l^H2n+2 J ^^^2n+2 (^^^2 _^ ^^^^2) (yjj 29) 

— oo 

Equations for A^, are given in (jVILlZp . (IVII.lSp . From (jVII.29p it is obvious that the self-consistent values 
Wri = const can be obtained only for a flat spectrum of instantons. However, with the flat specter and Bk 7^ 0, the 
second term in (|VII.29p creates a meaningless infinity. Therefore Bk = 0, and that, in turn, leads to the condition 
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(jVIL20[) . i.e. to quantum coherence of the instanton condensate. The quantitative characteristics of the condensate 
are formed by instantons only, locahzed in imaginary time. 
It is easy to calculate of the converging integrals in (|VII.29P : 



2 



J dxx^"+^ (^1 - e^" = 22^(2«-l)''(2« +!)(" + 2) • (VII.30) 

— oo 

After analytical continuation into the real space-time, following the rules (jVII.22[) — (|VII.24p . we obtain the final 
result: 

^ ^ _ l2^hN,r.st , I^„-^^^(2n- 1)1(271 + l)(n + 2) x^^^^^^ij2"+2, n^l, (VII.31) 



where 



N.nst - (rig) - (ugh) ■ (VII.32) 



The comparison of the two models of graviton-ghost condensate in the De Sitter space reveals some interesting 
features. In both cases we deal with the effect of quantum coherence. Expressions (|VII.3ip differ from (|VI.28p only in 
the formal substitution Ng Ninst ■ However the conditions leading to the quantum coherence are different in these 
models. According to (|VI.29p . in the condensate of virtual gravitons and ghosts, the average value of graviton and 
ghost occupational numbers are the same, and the non-zero effect appears due to the fact that the phase correlation 
in the quantum superposition in the graviton's and ghost's sectors are formed differently. As it follows from (|VII.20p . 
(|VII.32p . in the instanton condensate the phases in the graviton and ghost sectors correlate similarly, but the non- 
zero effect appears due to the difference of average occupation numbers for graviton's and ghost's instantons. The 
absence of the macroscopic structure of the condensates does not allow the detection of the differences by macroscopic 
measurements. In both cases the graviton-ghost vacuum possess equal energy-momentum characteristics. 

The question about the actual nature of the De Sitter space is lies in the formal mathematical domain. In these 
circumstances one should pay attention to the following facts. While describing the condensate of virtual gravitons 
and ghosts, we were forced to introduce an additional definition of the mathematically non-existent integrals (jVI.25p . 
i.e. to introduce into the theory some operations that were not present from the beginning. It is the additional 
operations that have provided a very specific property of the solution — the alternating signs in the sequence of the 
moments of the spectral function. By contrast, the theory of the instanton condensate has a completely different formal 
mathematics. The theory is motivated by the concrete property of the graviton-ghost system which is degeneration 
of quantum states, and the construction of the theory is constructed by the introduction of mathematically non- 
contradictory postulates. The moments of the spectral function's with alternating signs is an internal property of the 
graviton-ghost instanton theory. When we considered the instanton condensate in the De Sitter space, no additional 
mathematical redefinitions were necessary (compare the formulas (IVI.25I) and (|VII.30p ) . We have the impression that 
the instanton version of the De Sitter space is more mathematically comprehensive. Therefore, one may suggest that 
the key role in the formation of the De Sitter space (the asymptotic state of the Universe) belongs to the instanton 
condensate, appearing in the tunnelling processes between degenerated states of the graviton-ghost vacuum. 



C. Instanton Condensate of Constant Conformal Wavelength 

The exact solution (jV.41I II) has a pure instanton nature. Now we will obtain that solution with the value C^g ~ 0. 
One can rewrite the formulas (|VII.6p . (|VII.7p for the imaginary time: 

7^2^^1n^, a = aoe^p(^) . (VII.33) 
a-' ao V 4 y 



^ + ^0^^^ - fc^^- = 0, ^ + kli'-^ - = . (VII.34) 

As we already know, the spatially homogeneous modes participate in the formation of the solution for the equation 
(|V.41I II). As follows from (|VII.34p . when fc^ 0, the description of the spatially homogeneous modes in imaginary 
time does not differ from their description in real time. The contribution from modes g2 is present in (jVII.33p . with 
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the relations Cg4 = —k^ < 0, Cg2 = 3fcg/4 taken into account. These relations are necessary to provide the existence 
of the self-consistent solution. In what follows we are considering the quasi-resonant modes only. 

For k"^ — kg, the signs of the last terms in the equations (|VII.34p provide the existence of instanton solutions we 
are looking for: 



kcr 



i?k = 




hW^- 



16>ch 



ko 



-pkF(a)]ln^/^-, 
ao 



(VII.35) 



(VII. 36) 



where 



a 

F{a) = alj 



da 



a 



"■0 ao 



2a2 1ni/2. 



ao 



Calculations which follow contain the same mathematical operations we have already described several times in the 
previous sections. After we remove contributors to the spectral function which contains F{a), we obtain the condition 
for the coherence of the condensate. Some details of the calculations is related to the alternating signs of the moments, 
i.e. with the multiplier (—1)", characteristic for the instanton theory. Particularly, in the expression for Wi{g4:), there 
is a general sign "minus" . But, according to the Einstein equations in imaginary time Wi{g4:) > 0. The positive sign 
of the first moment is provided by the dominant contribution of ghost instantons over the contribution of graviton 
instantons. With that taken into account, we obtain the final equations for the moments of quasi-resonant modes, 
obtained after the analytic continuation into the real space-time: 



W„(g4) = (-1)^ 



In^ 

i/4„ 



-1 



a^" ao 



ao 



1/4, 



(VII.37) 



ao 



Here the following definition has been used: 

(nkigh)) - {nk(g)) 



2^2 



3fcoa§ 



The graviton instantons are dominant for the spatially homogeneous modes; 



(VII.38) 



The parameter of the spatially homogeneous condensate is defined as follows: 



("0(g)) (1 + Co^^ COSV?o) - {n(){gh)){^ + Co"' cosxo) 



iah) , 



90 j 



90 







From expressions (|VII.37|) . (|VII.38p . one gets energy density and pressure for the system of quasi-resonant and 
spatially homogeneous instantons: 
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(VII.39) 
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In formulas (|VII.39p . the terms in brackets are eliminated by the condition (|V.15|1 . which is rewritten in terms of 
macroscopic parameters 

^oiV^ = . (VII.40) 

Solutions (|VII.39|) , (|VII.40[) describe a quantum coherent condensate of quasi-resonant instantons with the ghost dom- 
inance. The parameters of the condensate are in accordance with parameters of a spatially homogeneous condensate 
with graviton dominance. 



VIII. GRAVITONS IN THE PRESENCE OF MATTER 



A. Nonlinear Representation of the BBGKY Chain and Integral Identities 



The full system of equations of self-consistent theory of gravitons in the isotropic Universe consists of the BBGKY 
chain (jV.7p — (jV.Op and macroscopic Einstein equations. In equations (|V.7p — (jV.9[) . the Hubble function H and 
its derivatives if, H are coefficients multiplied by the moments of the spectral function. In such a form the chain 
conserves its form even if besides of gravitons, other physical fields are also sources of the macroscopic gravitational 
field. We are interesting in the evolution of the flat isotropic Universe at a stage when the contributions of gravitons and 
non-relativistic particles, baryons and neutralinos, are quantitatively significant. (The latter are presumably carriers 
of the mass of Dark Matter.) We assume also that non-gravitational physical interactions created the equilibrium 
vacuum subsystems with full energy (an effective A-term) of the order of A ~ 3 • 10^^'' hr^ GeV*. The macroscopic 
Einstein equations containing all sources mentioned above read 

Rl^\R-^e,o. ^^ = ^^ + ^1 + 1 (a + ^) , (VIII.1) 



{£tot + Ptot) 



H 
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(VIII.2) 



Equation (jVIII.2|l should be differentiated with respect to time, and then D from (|V.7P should be substituted into 
the result of differentiation. These operations produce one more equation 



H = H 



-D + Wi 



3kM 



2a3 J 12 



Wi 



(VIII.3) 



The BBGKY chain (jV.7p — (|V.9p takes into account the interaction of gravitons with the self-consistent classical 
gravitational field which is represented by the Hubble function and its derivatives. According to Einstein equations 
(IVIII.ip — (|VIIL3p , a self-consistent gravitational field is created by gravitons and other components of cosmological 
medium, i.e. by the matter and non-gravitational vacuum subsystems. Therefore, the self-consistent gravitational 
field is a way of describing of significantly non-linear properties of the system that are the result of gravitational 
interaction of elements of the system. After excluding higher derivatives of the metric from the BBGKY chain (jV.Sp 
and (jV.gp . the true non-linear character of the theory emerges. Substitution of (jVIII.ip — (jVIII.Sp into (jV.Sp and 
(IV. 9p gives the non-linear representation of BBGKY chain: 



D + 6HD + AWi + IGHWi = , 

Wn +3(271 + S)HWn + 



16 



(4n^ + 671 + S)D +(n + lYWi + (8ri^ + ISn + 9)—^ + 2(2?!^ + 6n + 3)}<A 



(VIII.4) 



71 f 1 . 

-3 2^^+^ 



—D + {2v? + 371 + 3)T4^i + 871^ + IStj^ + 9)^5- + 4 271^ + 9n + 9)><A 
2 a-^ 

+4W'„+i + 8(71 + 2)HWn+i = 0, 71 = 1, 00 . 
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In the general case, the system of equations (|VIII.2[) and (jVIIL4[) (to which the definition a/a = H is added) should 
be solved numerically with initial conditions determined by the scale factor, moments of the spectral function and 
their derivatives 

a(0); D{0), W^„(0), Vk„(0), W^„(0), n = 1, oo . (VIIL5) 
The initial condition for the Hubble function should be calculated via the equation of the constraint (|VIII.1[) 



mO) = + jli.(0) + 1^,(0) -f i>.A + ^ . (VIIL6) 



Any solution of equations (|VIII.2[) and (|VIII.4p . which corresponds to initial conditions (jVIII.Sp . (|VIII.6[) . satisfies 
the identity which is local in time 

^^W-ii^W + ^^iW4><A+3^, (VIII.7) 
From the original equations of the chain (jV.7|) — (IV. 9|) one can obtain the chain of integral identities 

D{t) = -2Wi{t) + ^ i^D{0) + 2Wi{0) -2 J dtih'^ -^{h^Wi) 

{t t ti 

+ C„(,) / |i + C„(3) / f / (Villi 


t tl t2 

dti f dt2 f dt^ d /.2n+4w ^ I. 1 



oP' J a? J dt^ 





where 



a = ^ , C„(i) = WrM. C„(2) = Wn{Q) + 2ni/(0)W„(0), 
Cn{^) = + (An + 3)i/(0)W'„(0) + 2n[H{0) + {2n + 3)i/2(o)] W„(0) . 



(VIII.9) 



According to (jVIILSp . in the region attached to the point where initial conditions are defined, the solution can be 
represented in following form 

_lD«T^i^l/i(a) , Wn-^fnia), (VIII.IO) 

where /„(a) are slow-changing functions of the scale factor. Comparison of (jVIILSp with (jV.lip shows that the 
system tends to form a graviton-ghost condensate with constant conformal wavelength. 



B. Equation of State of Cosmological Medium Consisting of Dark Energy and Non-Relativistic matter. 

AGCDM Model 



In the framework of the theoretical model under discussion, the cosmological medium consists of three subsystems, 
each of these is described by its own energy density and pressure. They are the non-relativistic matter Emat = 
Mjc? , Pmat = 0; graviton-ghost coherent quantum condensate eg{a}, •pg{a\; and, possibly, the non-gravitational 
contribution of the A-term to the equilibrium vacuum energy. The original equations of the theory (jVIILip , (jVIII.2p 
can be written in the form (after simple and obvious transformations) 



ZH^ = >tetot{a\ = X ^£g{a} + A + 
-2H - SH^ = xptot{a} = x{Pg{a} - A) 



(VIII. 11) 



67 



where 

(VIII.12) 

are junctionals of scale factor. Explicit forms of these functionals are determined by the solution of the BBGKY 
chain in the nonlinear representation (|VIII.4p . The energy density and pressure of cosmological medium as a whole 
£tot{a}i Ptot{a} which are defined in (jVIILlip are also functionals of the scale factor. It is generally accepted that the 
combination of two formulae etot{i}, Ptot{a\ is called the equation of state of cosmological medium in parametrical 
form. 

The exact particular solutions of the equations of the theory are described in detail in Sections IV Bl IV CI IVI Al 
IVIBI I VII BIVII CI We will bear in mind these results during the discussion of the dependence of functionals (jVIII.12p 
on the scale factor. The most important is the result (jVIII.lOl) which follows from the integral identities (jVIILSp . 
which means that independently of concrete dynamics of the scale factor, the self -organization of the graviton-ghost 
medium leads to the formation of a condensate whose equation of state is of the form 

(VIII.13) 

C 

where Cg is a constant and /i(a) is a slow varying function of the scale factor. Approximate expressions (jVIILlSp can 
be obtained by a simple substitution of approximate relations ()VIII.10|) to the exact expressions (|VIII.12|) . According 
to (jVIILSp . the value and sign of the constant Cg and the form of the function /i(a) depend of initial conditions 
(IVIII.9[) . Let us mention also that the formation of the graviton-ghost medium of the approximate equation of state 
(IVIII.13I) is confirmed by numerical experiments that we conducted. 

The comparison of the characteristic dependence Eg (a) ^ from (|VIII.13p with the energy density of non- 

relativistic matter £mat(a) = Mjc? leads to two conclusions. First, the non-relativistic matter dominates over the 
condensate and the non-gravitational A-term during the epoch of sufhciently small scale factor, 

^»A. (VIII.14) 

Second, the epoch of the dominance of non-relativistic matter must be replaced by the epoch of condensate domination 
during which the inequality Cgf\{a)ld} ^ Mjct' is satisfied. This change of epochs is inevitable because in comparison 
to £mat{a), the energy density of graviton-ghost condensate £g{a) is a more slowly varying function of the scale factor, 
which increases with time. Further evolution of the cosmological medium depends on the relation between the energy 
density of the condensate and the value of the non-gravitational A-term. 

Let us assume that during the epoch of the dominance of the condensate over matter, the condensate (in some time 
interval) also dominates additionally over the A-term. This means that the following two inequalities are satisfied 
simultaneously 

%/i(a) » ^ , » A (VIII.15) 

In turns, this means that the cosmological model, which represents a medium consisting only of the condensate of 
the equation of state such a.s pg ~ eg/3, has good accuracy. This sort of equation of state follows from (|VIII.8p . 
(jVIILlOp and can be clearly seen from (|VIII.13p . Such a model (graviton-ghost condensate of constant conformal 
wavelength) was studied in Sections IV Bl IVI AlVII CI The exact solutions of the system of equations for the scale 
factor and graviton-ghost field obtained in these Sections are attractors. This mathematical status of these solutions 
is confirmed by numerical experiments. From this it follows that during the epoch of the Universe evolution when the 
inequalities (jVIII.14p are satisfied, the graviton-ghost condensate relaxes to the state in which its energy density and 
pressure are described by expressions that correspond to the attractor^ (see (|V.17I) . (|V.2ip . (IVI.lSp . (|VI.16p . (|VII.39p . 



In model l|VIII.16|l . the additive contribution of the spatial homogeneous condensate of the equation of state pg = £9 ~ is not 

taken into account. Because of the rapidly decreasing energy density with increasing scale factor, such a condensate can play some role 
only in the vicinity of cosmological singularity, while the epochs close to the contemporary Universe are of interest in this work. This 
is the reason why we also do not take into account the contribution of radiation whose equation of state is = Sr/S ~ l/a*. 
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(IVIL40I1 ): 



The sign of the Cg parameter in (jVIII.16[) depends on the microstructure of the graviton-ghost condensate. In the 
case of the condensate of virtual particles with graviton dominance Cg > 0, while in the case of instanton condensate 
with the ghost dominance Cg < 0. The type of condensate corresponding to the observable Dark Energy effect is 
determined by comparison of the theory with observational data. It will be shown in the next Section that Cg > 0. 

To discuss the future of the cosmological medium, it is necessary to bear in mind the fact of the existence of the 
De Sitter space. The geometry of this space is consistent with the vacuum states of all physical fields; it is stable 
and its symmetry is the highest among all possible symmetries. These are the reasons why the De Sitter space of the 
self-consistent geometry can be considered as an asymptotical state of the Universe. The results of comparison of the 
theory with observational data (see Section HX]) as well as internal properties of the theory itself are in favor of this 
assumption. The De Sitter solution of the BBGKY chain obtained in Section fV CI allows the interpretation in terms 
of virtual particles (Section I VI Bp as well as in terms of instantons (Section IVII Bp . The vacuum energy density in the 
asymptotic state e^^} = Aqo acquires the status of a fundamental cosmological parameter. The formulae to calculate 
eiri are given in (fOTjl . (|V33)l . (ESI. 

The self-polarized condensate in the De Sitter space is formed in the process of quantum-gravity phase transition 
(Section IV Dp . The energy density of graviton-ghost vacuum in the asymptotical state (with no A-term of non- 
gravitational nature) is 

247r2 

Aoo = ^TTT • (vni.i7) 

UNg 

If A ^ 0, it is necessary to make some assumptions on the absolute value and sign of A-term. For reasonable 
(relatively small) absolute values of the A-term, the theory proposes several more or less natural scenarios to form 
the asymptotical value e'^uac = Aqo consistent with existing observational data. In Section IV CI the possibility of 
equipartition of the vacuum energy between graviton-ghost and non-gravitational vacuum subsystems was shown. In 
this case, in accordance with (jV.38p . (jV.39p . the total vacuum energy density in the asymptotical state is 

(VIII.18) 

'^^^ g{max) 

If the absolute value of the non-gravitational A-term is significantly greater than the value which is acceptable 
from phenomenological considerations, then the scenario described in Section IVI Bl makes some sense. The theory 
predicts that the strong rcnormalization of the non-gravitational A-term of any sign by a graviton-ghost condensate 
(decreasing its original absolute value by several orders of magnitude) is possible. In this case, accordingly to (jVI.Sip . 
the asymptotical value of the total vacuum energy density is positive and can be estimated by the formula 



The results (|VIII.17p . (|VIII.18p . (|VIII.19p illustrate mathematically the existence of physical states of the condensate 
of constant energy density. From this fact it follows that under the following condition (when the inequality is satisfied) 

M 

eg{a} + A - Aoo > (VIII.20) 

the asymptotic behavior of the solution of the original equations (|VIII.lip describes the De Sitter space of constant 
4-curvature i?oo = — 4><Aoo. 

Thus, the internal mathematical properties of equations (jVIII.llI) allow the following classification of epochs of 
cosmological evolution: 

1. The epoch (jVIII.14p . during which the non-relativistic matter dominates over the graviton-ghost condensate 
and non-gravitational A-term. During this epoch the equation of state of the cosmological medium is etot{o) — 
Mja?' , ptot{a) — 0. Note that the existence of this epoch follows directly from WMAP data [11].) 

2. The epoch (jVIII.lSp . during which the graviton-ghost condensate of constant conformal wavelength dominates 
over the non-relativistic matter and non-gravitational A-term. During this epoch the equation of state of the 
cosmological medium is £tot(a) — Eg (a) , Ptotio) — Pgio-), where functions eg{a), Pg{a) are defined by expressions 
(IVIII.16I) . 
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3. The epoch (jVIIL20l) . during which the equiHbrium vacuum consisting of non-gravitational A-term and graviton- 
ghost condensate of constant physical wavelength dominates over the non-relativistic matter. During this epoch the 
equation of state of the cosmological medium is etotia) — ^oo , Ptot{o?) ~ — Aqo- 

It follows from inequalities (|VIII.14p . (|VIII.15p . (jVIII.20[) . that each epoch is clearly identified with the relevant 
region of the cosmological time scale. In these regions equations of state of the cosmological medium are very similar 
to equations of state obtained from exact solutions taking into account only the dominant component of the medium. 
An orderly sequence of exact solutions that take into account only the dominant terms, represents the most important 
properties of the exact solution of complete equations (jVIII.ll[) , augmented by BBGKY chain (|VIII.4p . Significant 
differences between the ordered sequence of exact solutions and solution of complete equations can be expected only 
in those parts of the cosmological scale where change of epochs takes place. 

The above presentation, based on the internal properties of the theory, clearly leads to the interpolation formulae 
for the energy density and pressure of cosmological medium: 

e*ot(a) = A^ + %ln^ + ^ , ^ _A«, - ^ In ^ . (VIII.21) 

The Universe filled by the medium with the equation of state (jVIII.2ip . described by Einstein's equations 



3— = >fetot(a) = >^ Aoo + ^ln— + — , (VIII.22) 

(X \ CL Qi CI / 



6- - ->c {etot{a) + 3ptot{a)) = k (ik^ _ % - ^ ) , (VIII.23) 



a 



In the terms of the quantities contained in the interpolation formulae (|VIII.2ip . the classification of epochs is performed 
according to the following inequalities and definition: 

^»A^ + %ln^ (VIII.24) 
a"^ a'^ a 

is the definition of the epoch of the matter dominance; 

^ In — > Aoo + ^ . VIII.25 

a a"^ 

is the definition of the epoch of the dominance of the condensate of constant conformal wavelength; 

A^»%ln^ + ^ (VIII.26) 

is the definition of the epoch of the dominance of the equilibrium vacuum. 

It is not difficult to notice the following fact. If one takes into account only the dominant component of the 
cosmological medium for each epoch, the streamlined set of approximate solutions obtained from (jVIII.22p . (jVIII.23p . 
practically coincides with the ordered collection of approximate solutions of the equations of one-loop quantum gravity, 
i.e. equations (|VIII.11|) . augmented by BBGKY chain (|VIII.4|) . 

The exact solution of equations (|VIII.22[) , (|VIII.23|) interpolates the exact solution of equations l|VIII.ll[) , (|VIII.4p . 
Differences between the exact solutions may be visible only in narrow transitional areas at the intersections of epochs. 
The most noticeable differences occur at the transition from the era of dominance of the condensate of a constant 
conformal wavelength to the era of dominance of the equilibrium vacuum. As was mentioned in Section FVDi numerical 
experiments show that the transition from pre-asymptotic stage to the asymptotic stage is followed by nonlinear 
fluctuations of physical observables. The latter forces us to assume that the asymptotic state of the graviton-ghost 
condensate is formed as a result of quantum-gravity phase transition on the scale of the whole Universe. It is worth 
mentioning that comparison of the model (|VIII.22p . (|VIII.23I) with the observed data shows that at the stage (|VIII.25p 
(or what is the same, at the stage (jVIII.lsp ) the Universe expanded with deceleration. Thus, the currently observed 
acceleration of expansion may be a consequence of a phase transition into an asymptotic state, i.e. a state of 4-space 
of constant curvature, with the highest symmetry among all possible symmetries. 

The issue of experimental verification of the phase transition is very complicated, and we will not discuss it in this 
paper. It is worth mentioning only that, in order to clarify the issue, it is necessary to conduct detailed computer 
experiments. However, comparison of the results of computer experiments with real observational data can be done 
only on the condition that the observational data to meet very rigid standards of completeness and accuracy. 
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At the level of accuracy of the modern cosmological observations, transition processes can be neglected. Thus, it 
gives interpolation formulae (|VIII.2ip the status of a model which follows only from the first principles of the one- 
loop quantum gravity. These principles are three: it (i) the inevitable appearance of the sector of nontrivial ghosts 
interacting with the macroscopic gravity; (ii) conformal non-invariance and zero rest mass of gravitons and ghosts; 
(iii) one-loop finiteness off the mass shell. 

The difference between the model (|VIII.2ip and the known ACDM model is in the additional term which is relevant 
to the variable component of the Dark Energy. We propose the name "AGCDM model" for the Graviton-Cold- 
Dark-Matter model (jVIII.21[) . We assume that the equations (|VIII.2ip are able to represent the energy density and 
pressure of the cosmological medium after the separation of matter and radiation, including the epoch of the formation 
of galaxies and the contemporary Universe, including its future. 

Currently, the priority task is the experimental verification of the variable component of Dark Energy and accordance 
with the properties of this component predicted by the AGCDM model. 



IX. GRAVITON-GHOST CONDENSATES IN THE ROLE OF DARK ENERGY. COMPARISON OF 

THEORY WITH OBSERVATIONAL DATA 

A. Approaches to the Dark Energy Problem 

Acceleration of the Universe and, consequently, the existence of Dark Energy, causing the phenomenon, was dis- 
covered experimentally in the works 0| . Today, we can identify three aspects of the Dark Energy problem. The 
first is obtaining reliable observational data on the density of Dark Energy £0^(0) as a function of the scale factor, 
or which is the same thing, as a function of the redshift z=:l/a— 1. It is well known that the simplest Dark Energy 
model is the cosmological constant. Observational aspect of the problem of Dark Energy lies in the answer to the 
question is it true that the density Eqe = A = const or does Dark Energy contain a variable component? From the 
point of view of available data, some preferred models are those with variable component, but, of course, further 
detailed study is needed (see, e.g., the SNAP project [H^l). The second aspect is that the physical nature of Dark 
Energy is a problem of fundamental physics. The statement of the problem will be determined to a substantial extent 
by existing observational data (that we expect is to be provided with the necessary completeness and accuracy) . The 
third aspect is the astrophysical aspect. It will become particularly relevant if the variable component of Dark Energy 
will be reliably identified during the era of large-scale structure formation in the Universe. 

It is common to consider that the problem of the physical nature of Dark Energy is a window into new physics. 
These expectations suggest that the theoretical model of Dark Energy is not limited to the cosmological constant. A 
cosmological constant could hardly be called new physics. Indeed, the most common principles of geometrized theory 
of gravitation say only that in a weak gravitational field the Lagrangian of the theory should be presented in a series 
of powers of invariants of the curvature. From this perspective, A-term is a zero term of expansion, and the scalar 
curvature, generating the left hand side of Einstein's equations, is the first term of the expansion of the Lagrangian 
of the theory of gravity. In other words, the presence of A-term in Einstein's equations is not a conceptual problem. 
On the contrary, the theory of gravity, taking into account the A-term looks more natural than the theory without it. 

The problem is the numerical value of A-term (it is anticipated that the value of A-term can be calculated in 
principle) . In this regard we wish to point out that the problem of calculating of A-term conceptually does not differ 
from the problem of calculating other fundamental constants of physics, e.g., electron charge, the Fermi constant of 
weak interaction, scale of Quantum Chromodynamics and others. Variations of these constants, as is known, also 
lead to a change in the properties of the Universe. The A-term problem is singled out rather emotionally. With the 
variation of the cosmological constant catastrophic changes in the Universe will follow, as can be seen by a simple 
analysis of a simple equation. We will continue the discussion of A-term in the Appendix IXII 

Possible deviation of the Dark Energy density of the cosmological constant is being actively studied by theorists 
(for a catalogue of models see review [goj) and is the subject of planned experiments (see SNAP project in ^^). A 
quintessence (evolving scalar field) was widely discussed in [6l|, [bI, [gI, [H, [63 • A free scalar field of a small mass as a 
model of Dark Energy was discussed in [6^ . Models with violation of the weak energy condition ("phantom energy" 
) were proposed in (btI . [68j . Models with Chaplygin gas and its modifications were considered in [63 | (see review [63| 
and references therein). The massless scalar field and Casimir effect were proposed in (toI . [71|. Some brane world 
models were proposed in [tI, [tI, [tI. [7p . The discussion of these and other ideas and models, including the "new 
physics" proposals can be found in [60l . [tgI . Wf\ and other reviews. All of these ideas and approaches have no 
bearing on this paper. Current prospects for the dynamical theory of the graviton vacuum with applications to the 
Dark Energy problem can be found in the review [3] ■ Technical approaches to the vacuum energy and the problem of 
cosmological constant in terms of quantum field theory in curved spacetime can be found in review [43|. The state of 
the art of observational data and theoretical prospects for the Dark Energy problem can be found infTSl. Itq} . 
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The problem of choosing a model for dealing with experimental data becomes very relevant. In essence, the choice 
of a model is the choice of the physical interpretation of the nature of Dark Energy. The problem is complicated by 
a lack of completeness and accuracy of data available on SNIa. The intervals of redshift < z < 1.8, available from 
observations of SNIa, are too narrow to choose between alternative models. In that situation, priority is given to 
models based on common and experimentally motivated principles of fundamental physics. Of course, A-term belongs 
to the first-principle models. In our view, a model based on the effect of condensation of quasi-resonant graviton and 
ghost modes in the self-consistent field of the isotropic Universe has the same status. The emergence of graviton-ghost 
condensates in the process of cosmological evolution is a consequence of only the most general properties of one-loop 
quantum gravity. These are conformal non-invariance and zero rest mass of gravitons; inevitability of appearance of 
the ghost sector; and the one-loop finiteness off the mass shell. Interpolation formulae of AGCDM model (|VIII.2ip . 
describing this effect, are based on exact solutions and are justified by integral identities (jVIILSp . 

With Cg^ AGCDM (|VIII.2ip turns into ACDM model. It is therefore obvious that the model of Dark Energy, 
based on one-loop quantum gravity, is consistent with existing observational data as well as the A-term model. In 
this situation, one must first determine whether the AGCDM model has statistically significant advantages over the 
ACDM model. 



B. Observational Data 



The set of formulas used for observational data processing follows below. The normalization of the scale factor in 
equations (jVIII.22p . (|VIII.23p can always be chosen in such a way that a{tu) = 1 in the contemporary Universe. To 
normalize other physical quantities, we use the Hubble constant in the contemporary epoch Hq — lOOh km/(s-Mpc), 
where h = O.lStom- 

In the dimensionless form, equations of AGCDM model read 



where 



^ H _ a ^ Q 



The values of parameters 



otlQ 'J^O "^^0 



S'A/ — '-'•^^-0.04' " ^ '-'•''-'-0.03 



are known from WMAP data [58[ . Below are expressions for the density and pressure of Dark Energy: 

enE{z) = eg{z) = f^A + f}g(l + zf ln[ao(l + z)] , 

PnE{z) = Pg{z) = -f^A - Ulgil + zf (ln[ao(l + z)] - 1) . 



(IX.1) 



(IX.2) 



In the flat Universe, four parameters of AGCDM model are connected by the following condition 

+ VLg In aa + n^, = l. (IX.3) 

The information on the normalized density of Dark Energy is contained in the dependence of conformal distance 
to supernovae SNIa of their redshifts. It reads 



R{z) 



dz 

W) ' (IX.4) 



H{z) = ^Jn^,{\ + zf + eaB(z) , 
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Figure 1: Distant modulii of supernovae SNIa. Observational data — 280 
counts points from [s oj. Fitt in g curv es: 1 — ACDM model (iDTg]) . (|IX.12|) : 
2 — AGCDM model (|IX.14|) . HIX.15|) . 



Observational data on SNIa are usually presented in terms of distant moduli, which reads 

tis{z) ^ M + 5\g[{l + z)R{z)], 
M = -Slgft +42.3841 + (A/b + 19.31) , 



(IX.5) 



where is the absolute bolometric brightness of a standard candle. The full set, containing data on all known 
supernovae (292 points), is given in [801 . In [8l|, a "Gold Data Set" was built from the data of [s^, [s^]^- In works 
iS HH, m, it is anticipated that = -19.31. In the equation (|IX3|) . M = 43.09^°;°2 for h, Ms values given 
above. Actually, however, questions about the precise values of normalized Hubble constant and the brightness of a 
standard candle remain open. For this reason, values /i^ (z) quoted in [8l|, [s^l are treated as data in arbitrary 
units. The value of M is determined in the process of processing these data. 

A set of observational data on supernovae SNIa is presented in Figure 1 (280 points selected from 292 points 80]). 
Solid lines in Figure 1 show a fit corresponding to the ACDM model (curve 1) and AGCDM model (curve 2). 

In the verification procedure of theoretical models an important role plays the redshift zq, at which the deceleration 
of the expansion of the Universe is changed to an accelerated expansion. The value below is recovered from the Hubble 
diagram 



From the condition Q{zq) 
reads 



1 + ^0 = 1-46 ± 0.13 . (IX.6) 
for a given value of 1 + zg, one more condition is obtained for the AGCDM model which 



2r!A - ng{l + zo)2 - r!,,(l + zo)^ = 



(IX.7) 



Note the following circumstance. If one fixes and 1 + parameters by their average values and substitutes these 
values into pX.3p and pX.7p . then only one independent parameter remains in the AGCDM model. The value of 
this parameter will be found by fitting of the Hubble diagram as a whole. It is convenient to select ^Ig as such an 



New versions of SNIa data see in |83l Is^l 
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independent fit parameter characterizing the variable component of Dark Energy. We hope that in the future, with 
increasing accuracy of measurements of iljj and 1 + zq the implementation of this program becomes possible. 

Currently, we can use intervals of acceptable values of parameters ^Im and 1 + zq to assess the status of the model 
of A~term. Assuming that fig = in (jlX.Sp and pX.7[) . we get 

0.73 < Qa < 0.8, (1 + zo)acdm = (^^) ^ =1.85^^;!^ (IX.8) 

As we can see, la interval predictions of ACDM model (jIX.8p does not overlap the la interval of measured values 
PX.6|) . On its own, this fact does not reject the A-term as a carrier of Dark Energy. For this, it would be necessary 
to detect the contradiction at the level of more than three standard deviations. In addition, the contradiction at the 
level of la is smoothed out if one takes into account systematic errors. Discussion of this range of issues can be found 
in the Weinberg book [78|. Nevertheless, when comparing (|IX.6|) and pX.8p doubts about the adequacy of ACDM 
model emerge. In this situation, consideration of AGCDM model strongly indicated. 



C. Results of the Data Fit and its Discussion 



Below, the best-fit parameters of multi-parameter models were calculated using the standard LSM statistical 
approach based on the calculation of the Fisher matrix [t^, [s^. To get the sense of the error margins for the 
parameters as they come from the errors in the current observational data, we show the square root of diagonal 
elements of the covariance matrix (1-sigma). Of course, for any further use of the calculated parameters the full 
covariance matrix must be used to take into account the correlation between parameters. 

We evaluate the AGCDM model in the comparison with the ACDM model. Therefore, first, let us introduce the 
results of the fit of the Hubble diagram using formulas of the ACDM model. The value of pX.6[) is included in the 
database of the fitting data, and the interval of allowed values of the parameter is defined by WMAP data. The 
results of the fit are shown on Figure 1, curve 1. The values of fitting parameters are: 



0.27 ±0.03 



0.73 ±0.03, (1 



1.75 ±0.07 



(IX.9) 



X = 43.37 ±0.02 , x^/dof^ 0.97 



We will use also the value of the statistical sum and the formula to calculate the number of degrees of freedom with 
their statistical fluctuations taken into account: 



1 



Zofit — Zoexp 



iV+ 1 



± 



xVdof 



2S 



N + 1 



(IX. 10) 



N + 1 



In (lIX.lOp . + 1 — n is the number of degrees of freedom; N is the number of experimental points of distant moduli 
fJ-exp{zi); term "1" corresponds to the experimental value 1 + zq; n is the number of independent fitting parameters 
(which are fl^ and Ai for ACDM model) . The interval of values of the statistical sum, characterizing the fit ACDM 
model (in the framework of appropriate ACDM's formulae) is calculated as follows: 



S±AS^xVdof 



N + l-n 



± 



N + 1 



(IX.U) 



{S ± AS)acdm = 136 ± 11 . (IX.12) 

The question whether the existing set of experimental data could pinpoint a single preferable model is reduced to 
the comparison of the intervals of the statistical sums. Let us assume that there is a set of models Xa, a — 1 ,2 , M, 
and the results of the fits by the formulas of each model are characterized by the respective intervals of the statistical 
sums Sa ± ASa- The models Xa and Xb are statistically equivalent if the intervals of the statistical sums Sa ± ASa 
and Sb ± ASb overlap. If the intervals do not overlap, the statistically preferable model is a model that corresponds 
to the minimal /dof. Certainly, in addition to the statistical criteria, it is necessary to take into account how the 
models correspond to general principles, how hypothetical the concepts are, etc. Assuming that the model of the 
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A-term is the simplest one among all of the first-principle models, the interval of values of the statistical sum (jIX.12p 
could be treated as the reference interval. 

The fit of data in Figure 1 using equations of the AGCDM model shows that the results practically do not depend 
on the parameter fl^ chosen within the interval 0.2 < Hm < 0.27 which is in accordance with WMAP data. As is 
shown below, that observation is connected with a specific prediction of the AGCDM model: in the area of the Hubble 
diagram for supernovae SNIa, 0.007 < z < 1.755 , the density of the Dark Energy is on the par with the density of 
non-relativistic matter. In that situation, Qm cannot have the status of fitting parameter. We have conducted three 
fitting procedures with three fixed values of = 0.20; 0.24; 0.27. In all three cases 

X^/dof ^0.91, Al = 43.33 ± 0.03 , 1 + zq 1.34 ± 0.04 . (IX.13) 

The values of other parameters are as follows: 

n„=0.20: f^A = 2.25 ±0.53, rjg = 2.23 ± 0.68 , oo = 0.52 ± 0.02 ; 

f7„=0.24: ^ 2.23 ±0.54 , = 2.15 ± 0.68 , oq 0.51 ± 0.02 ; (IX.14) 

O„=0.27: f^A = 2.21 ±0.55 , rjg = 2.09 ± 0.68 , ao 0.49 ± 0.02 . 

The interval of values of the statistical sum is calculated by equations (jlX.lOp and (jlX.lip with n — 3 (the fitting 
parameters are fig, ao, M): 

{S ± AS) AGCDM ^127 ±11 . (IX.15) 

As follows from (jIX.12p and (jIX.15|l , the intervals for the statistical sums for ACDM and AGCDM models are largely 
overlapping. Therefore, from the point of view of available observational data, these models should be considered 
statistically equivalent, with small advantage of AGCDM model over ACDM in accordance with the x^-criterion. 

ACDM and AGCDM models start diverging from each other with increasing completeness and accuracy of ob- 
servational data. We can even make them statistically distinguishable from each other on the basis of available 
observational data by data smoothing (using a simple operation of noise reduction). The first smoothing operation 
is conducted by averaging of distant moduli with the redshift conserved. In the second step, at the interval up to 
z = 1.4, the averaging is done by a moving window, the size of three neighboring points. Thus, the errors in the 
averaged data are reduced by a factor of v^. To maintain the statistical weight of the transition point, its error is 
also reduced by a factor of \/3: 1 + zq = 1.460 ± 0.075. 211 points obtained by that operation, iJexp{zi) are shown in 
Figure 2. The fit of smoothed points by the formulas of ACDM model provides the following results: 

f^M = 0.27 ±0.02 , rJA = 0.73 ± 0.02 , (1 + zq) = 1.75 ± 0.05 , 

(IX. 16) 

43.34 ±0.01 , xV'^o/=1.15 



{S±AS)i^j,M = ^21±12. (IX.17) 

With the fit of the same points by using the AGCDM model, again we observe the insensitivity of the x^-criterion to 
the variations of the parameter within the experimentally allowed interval: 

xVrfo/ = 0.93, = 43.32 ±0.02 , 1 + zq = 1.34 ± 0.03 . (IX.18) 

The values of other parameters change slightly with variation of the parameter Qm- 

r2„=0.20: f^A = 2.19 ±0.35, f7g 2.16 ± 0.45 , oq = 0.53 ± 0.02 ; 

r2„=0.24: Oa = 2.16 ±0.36 , f7g = 2.08 ± 0.45 , ao = 0.51 ± 0.01 ; (IX.19) 

f^M = 0.27: f^A = 2.15 ± 0.36 , f7g = 2.02 ± 0.45 , ao = 0.50 ± 0.01 . 



The statistical sum of the smoothed data has the following value and deviation: 

(5±A5)f^CDM-97.6±9.5 . 



(IX.20) 
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Figure 2: Distant modulii of supernovae SNIa. Smoothed observational data 
— 211 points. Fitting c urves: 1 — ACDM model pX.16p . pX.17|l : 2 — 
AGCDM model (|IX.19|) . (|IX.20p . 
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Figure 3: Density and pressure of Dark Energy and density of non-relativistic 
matter in the region of observations of SNIa and in the future Universe. 

Com paring the results (jIX.9j) and (llX.iep . belonging to ACDM model with the results (|IX.12p . (IIX.13P and (jlX.lSp . 
(jlX.igp . belonging to the AGCDM model, it is easy to see that in the framework of the given model the values of 
parameters obtained from the fits of the original and smoothed data, are statistically equivalent. This fact is not 
something amazing because the smoothing procedure and the fitting are mathematically related. The use of smoothed 
data is motivated by the necessity to compare and make a preliminary choice between two different models using the 
criterion of statistical likehhood. 



76 




Figure 4: Acceleration of cosmological expansion in the region of observations 
of SNIa and in the future Universe. 



When the two models are compared, two facts appear. First, the application of the x^-criterion to the smoothed 
data of the AGCDM model shows that the model gains a more significant advantage. Second, comparison of (jlX.lZp 
and (|IX.20[) shows that on the smoothed data the intervals of statistical sums are not overlapping, with the lesser value 
of the statistical sum belonging to the AGCDM model. Therefore, by statistical criteria obtained from the smoothed 
Hubble diagram for supernovae SNIa, the AGCDM model has an advantage. 

In addition to the Hubble diagram for supernovae SNIa, the information about Dark Energy is contained in the 
Hubble diagram for radio-galaxies [s^ and gamma-ray bursts [s^] . It is also contained in the cosmological parameters 
extracted from the CMB data and correlation functions characterizing the large scale structure of the Universe. Results 
of processing of the full data set will be reported in a separate paper. Because of big statistical errors of appropriate 
data, we believe that the use of Hubble diagrams for radio-galaxies and gamma-ray bursts is to no purpose in this 
work where we discuss statistical criteria to choose between ACDM and AGCDM models. As to the third type of 
data (following from CMB and correlation functions), we would like to mention the following. Information about the 
density of Dark Energy on the cosmological scale from the instant of the last scattering Zis up to the present time is 
contained in the shift parameter [s^ 



n 




= 1.716 ±0.062 



(IX.21) 



The (jIX.21[) parameter can be used under condition that the Dark Energy model, used for the interpretation of SNIa 
data in the < z < 2 interval, can be extrapolated up to z ~ 1000. In our view, such an extrapolation (three 
orders of magnitude for redshifts and nine orders of magnitude for curvature and energy density) looks improbable. 
In any case, the interpolation AGCDM model is not allowed to be extrapolated in such a way as can be seen from 
the algorithm of its construction (see Section IVIII Bp . This is the reason why we do not discuss constraints following 
from (|IX.2ip . The parameter A extracted from the acoustic peak data is of a different status. Acoustic oscillations 
in the photon-baryon plasma prior to recombination give rise to a pick in the correlation function of galaxies. This 
effect was recently been measured in a sample of luminous red galaxies and leads to the value [88] 



A = 



Zlrc 



dz 



1/3 



= 0.469 



Vo.98/ 



-0.35 



±0.017 , 



(IX.22) 



where Z^hg = 0.35 is the redshift at which the acoustic scale has been measured and n 
of the primordial power spectrum. 



0.961 is the spectral index 
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The AGCDM model with the pX.lOp parameters, which were determined from the SNIa data, predicts the following 
interval for the allowed values of the A parameter 0.413 < A < 0.479 if the parameter r2„ is taken from the allowed 
interval 0.20 < Hm < 0.27. The AGCDM model, which is precisely consistent with the value A = 0.472, is as follows 

X^/dof:^O.m, 7W = 43.32 ± 0.02 , 1 + zq = 1.34 ± 0.03 , 

(IX.23) 

= 0.26 ±0.01 : r^A = 2.15 ±0.27 , f^g = 2.04 ± 0.40 , ao = 0.50 ± 0.02 . 

The procedure of noise reduction in the Hubble diagram and the obtained results bear, of course, an illustrative 
character. In fact, we have demonstrated only the near term possibilities of the experiment: i.e. possibility to get 
the variable component of the Dark Energy and compare the results of the measurements with the prediction of the 
AGCDM model derived from the exact results of one-loop quantum gravity. 

The density and pressure of the Dark Energy and the energy density of non-relativistic matter are plotted in Figure 
3. The graphs are calculated by the formulas of AGCDM model (jIX.2[) with parameters corresponding to = 0.240 
in (jlX.igp . The acceleration of the Universe, calculated by the formula pX.l[) with the same parameters, is shown 
in Figure 4. Two features on the graphs in Figure 3 are of special interest: (i) in the entire area of observation of 
SNIa the density of non-relativistic matter and Dark Energy are compatible in their values; (ii) the density of Dark 
Energy first decreases with decreasing z, and then increases, starting approximately from z ^ 0.2. Let us mention the 
astrophysical aspect of these predictions of the AGCDM model: (i) at the epoch of creation of large scale structures 
in the Universe, Dark Energy has played a quantitatively important role in the global cosmological dynamics; (ii) 
the area of relatively small values of redshifts represents significant interest for the observations, because in that 
area the reconstruction of the graviton-ghost condensate is happening: pre-asymptotic state of the condensate with 
constant conformal wavelength is transforming into the asymptotic state of the condensate with the constant physical 
wavelength. With the increasing accuracy of experiments, it is possible that in this area, the nonlinear fiuctuations 
of Dark Energy could be discovered (see Sectioi fVDj) . 

X. CONCLUSION 

From the formal mathematical point of view, the above theory is identical to transformations of equations, deter- 
mined by the original gauged path integral (jll.l[) . leading to exact solutions for the model of self-consistent theory 
of gravitons in the isotropic Universe. To assess the validity of the theory, it is useful to discuss again but briefiy the 
three issues of the theory that are missing in the original path integral. 

(i) The hypothesis of the existence of classic spacetime with deterministic, hut self-consistent geometry is introduced 
into the theory. It is not necessary to discuss in detail this hypothesis because it simply refiects the obvious experimen- 
tal fact (region of Planck curvature and energy density is not a subject of study in the theory under discussion). Note, 
however, that the introduction of this hypothesis into the formalism of the theory leads to a rigorous mathematical 
consequence: the strict definition of the operation of separation of classical and quantum variables uniquely captures 
the exponential parameterization of the metric. 

(ii) The transfer to the one-loop approximation is conducted in the self-consistent classical and quantum system of 
equations. Formally, this approximation is of a technical nature because the equations of the theory are simplified only 
in order to obtain specific approximate solutions. After classical and quantum variables are identified, the procedure 
of transition to the one-loop approximation is of a standard and known character p7j . In reality, of course, the 
situation in the theory is much more c omp lex and paradoxical. On the one hand, the quantum theory of gravity 
is a non-renormalized theory (see, e.g. [83]). Specific quantitative studies of effects off one-loop approximation are 
simply impossible. On the other hand, the quantum theory of gravity without fields of matter is finite in the one-loop 
approximation 5] . The latter means that the results obtained in the framework of one-loop quantum gravity pose 
limits to its applicability that is mathematically clear and physically significant. The existence of a range of validity 
for the one-loop quantum gravity without fields of matter is a consequence of two facts. First, there are supergravity 
theories with fields of matter which are finite beyond the limits of one-loop approximation. Second, the quantum 
graviton field is the only physical field with a unique combination of such properties as conformal non-invariance 
and zero rest mass. For this field only there is no threshold for the vacuum polarization and particle creation in the 
isotropic Universe. Therefore, in the stages of evolution of the Universe, where i/^, \H\ <C m"^ {m is mass of any of 
the elementary particles), quantum gravitational effects can occur only in the subsystem of gravitons. It is also clear 
that in any future theory that unifies gravity with other physical interactions, equations of theory of gravitons in 
one-loop approximation will not be different from those we discuss in this work. Therefore the self-consistent theory 
of gravitons has the right to lay claim be a reliable description of the most significant quantum gravity phenomena in 
the isotropic Universe. 
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(iii) Dynamic properties of ghost fields are captured by the condition of one-loop finiteness of the theory off mass 
shell of gravitons and ghosts. The class of legitimate gauges picked out by this condition includes gauges that are form- 
invariant with respect to transformation of the symmetry group of the background geometry. This point is the most 
nontrivial part of the theory because it is essentially an additional mathematical condition on the theory ensuring its 
internal consistency. The condition of one-loop finiteness off the mass shell largely determines the mathematical and 
physical content of the theory. Given that the main results of this work are exact solutions and exact transformations, 
the evaluation of he proposed approach is reduced to a discussion of this point of the theory. Let us enumerate once 
more logical and mathematical reasons, forcing us to include the condition of one-loop finiteness off the mass shell 
into the structure of the theory. 

a) Future theory that will unify quantum gravity with the theory of other physical interactions may not belong to 
renormalizability theories. If such a theory exists, it may only be a finite theory. One-loop finiteness of quantum 
gravity with no fields of matter that is fixed on the mass shell can be seen as the prototype of properties of the 
future theory. 

b) Because of their conformal non-invariance and zero rest mass, gravitons and ghosts fundamentally can not be 
located exactly on the mass shell in the real Universe. Therefore, the problem of one-loop finiteness off the mass shell 
is contained in the internal structure of the theory. 

c) In formal schemes, which do not meet the one-loop finiteness, divergences arise in terms of macroscopic physical 
quantities. To eliminate these divergences, one needs to modify the Lagrangian of the gravity theory, entering 
quadratic invariants. This, in turn, leads to abandonment of the original definition of the graviton field that generates 
these divergences. The logical inconsistency of such a formal scheme is obvious. (The mathematical proof of this 
claim is contained at Appendix IXII Bl ) 

d) In the self-consistent theory of gravitons, one-loop finiteness off the mass shell can be achieved only through 
mutual compensation of divergent graviton and ghost contributions in macroscopic quantities. The existence of gauges, 
automatically providing such a compensation, is an intrinsic property of the theory. 

From our perspective, the properties of the theory identified in points a), b), c) and d), clearly dictate the need to 
use only the formulation of self-consistent theory of gravitons, in which the condition of one-loop finiteness off the 
mass shell (the condition of internal consistency of the theory) is performed automatically. We also want to emphasize 
that, as it seems to us, the scheme of the theory given below has no alternative both logically and mathematically. 

Gauged path integral factorization of classic and quantum variables, which ensures the existence of a self- 
consistent system of equations =^ transition to the one-loop approximation, taking into account the fundamental 
impossibility of removing the contributions of ghost fields to observables =4» choice of the ghost sector, satisfying the 
condition of one-loop finiteness off the mass shell — appears to us logically and mathematically as the only choice. 

As part of the theories preserving macroscopic spacetime being clearly one of its components, we see two topics 
for further discussions. The first of these is the replication of the results of this work by mathematically equivalent 
formalisms of one-loop quantum gravity. Here we can note that, for example, in the formalism of the extended 
phase space with BRST symmetry, our results are reproduced, even though the mathematical formalism is more 
cumbersome. The second topic is the reproduction of our results in more general theories than the one-loop quantum 
gravity without fields of matter. Here is meant a step beyond the limits of one-loop approximation as well as a 
description of quantum processes involving gravitons, while taking into account the existence of other quantum fields 
of spin J ^ 3/2. In the framework of discussion on this topic, we can make only one assertion: in the one-loop 
N = 1 supergravity containing graviton field and one gravitino field, the results of our work are fully retained. This 
is achieved by two internal properties of = 1 supergravity: (i) The sector of gravitons and graviton ghosts in this 
theory is exactly the same as in the one-loop quantum gravity without fields of matter; (ii) The physical degrees of 
freedom of gravitino with chiral h = ±3/2 in the isotropic Universe are dynamically separated from the non-physical 
degrees of freedom and are conformally invariant; (iii) The gauge of gravitino field can be chosen in such a way that 
the gravitino ghosts automatically provide one-loop finiteness of A^ = 1 supergravity. As for multi-loop calculations 
in the N = 1 supergravity and more advanced theoretical models, we have not explored the issue. 

Of course, a rather serious problem of the physical nature of ghosts remains. The present work makes use in practice 
only of formal properties of quantum gravity of Faddeev-Popov~De Witt, which point to the impossibility in principle 
of removing contributions of ghosts to observable quantities off the mass shell. A deeper analysis undoubtedly will 
address the foundations of quantum theory. In particular, one should point out the fact that the formalism of the path 
integral of Faddeev-Popov-De Witt is mathematically equivalent to the assumption that observable quantities can 
be expressed through derivatives of operator-valued functions defined on the classical spacetime of a given topology. 
On the other hand, finiteness of physical quantities is ensured in the axiomatic quantum field theory by invoking 
limited field operators smoothed over certain small areas of spacetime. Extrapolation of this idea to quantum theory 
of gravity immediately brings up the question on the role of spacetime foam [IS] (fiuctuations of topology on the 
microscopic level) in the formation of smoothed operators, and consequently, observable quantities. To make this 
problem more concrete, a question can be posed on collective processes in a system of topological fiuctuations that 
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form the foam. It is not excluded that the non removable Faddeev-Popov ghosts in ensuring the one-loop finiteness 
of quantum gravity are at the same time a phenomenological description of processes of this kind. 

Study of equations of self-consistent theory of gravitons, automatically satisfying the condition of one-loop finite- 
ness, leads to the discovery of a new class of physical phenomena which are macroscopic effects of quantum gravity. 
Like the other two macroscopic quantum phenomena of superconductivity and superfluidity, macroscopic effects of 
quantum gravity occur on the macroscopic scale of the system as a whole, in this c;asc, on the horizon scale of the 
Universe. Interpretation of these effects is made in terms of gravitons-ghost condensates arising from the interference 
of quantum coherent states. Each of coherent states is a state of gravitons (or ghosts) with a certain wavelength of the 
order of the distance to the horizon and a certain occupation number. The vector of the physical state is a coherent 
superposition of vectors with different occupation numbers. 

A key part in the formalism of self-consistent theory of gravitons is played by the BBGKY chain for the spectral 
function of gravitons, renormalized by ghosts. It is important that equations of the chain may be introduced at an 
axiomatic level without specifying explicitly field operators and state vectors. It is only necessary to assume the 
preservation of the structure of the chain equations in the process of elimination of divergences of the moments of 
the spectral function. Three exact solutions of one-loop quantum gravity are found in the framework of BBGKY 
formalism. The invariance of the theory with respect to the Wick rotation is also shown. This means that the solutions 
of the chain equations, in principle, cover two types of condensates: condensates of virtual gravitons and ghosts and 
condensates of instanton fluctuations. 

All exact solutions, originally found in the BBGKY formalism, are reproduced at the level of exact solutions for 
field operators and state vectors. It was foimd that exact solutions correspond to various condensates with different 
graviton-ghost microstructure. Each exact solution we found is compared to a phase state of graviton-ghost medium; 
quantum-gravity phase transitions are introduced. 

We siispect that the manifold of exact solutions of one-loop quantum gravity is not exhausted by three solutions 
described in this paper. Search for new exact solutions and development of algorithms for that search, respectively, is a 
promising research topic within the proposed theory. Of great interest will also be approximate solutions, particularly 
those that describe non-equilibrium and unstable graviton-ghost and instanton configurations. 

Self-consistent theory of gravitons allows an easy generalization that takes into account participation of non- 
relativistic matter in the formation of common self-consistent gravitational field. From the equations of this theory it 
follows that the era of dominance of non-relativistic matter should be replaced by an era of dominance of graviton- 
ghost condensate. This result is of direct relevance to the physics of Dark Energy. The pre- asymptotic condition of 
Dark Energy is interpreted as a condensate of virtual gravitons and ghosts of a constant conformal wavelength. As an 
asymptotic condition, the theory predicts self-polarized graviton-ghost condensate of constant physical wavelength 
in the De Sitter space. 

The view of the nature of Dark Energy is formulated in the form AGCDM model which interpolates exact solutions 
of the one loop quantum gravity. The proposed theory is consistent with existing observational data on Dark Energy 
extracted from the Hubble diagram for supernovae SNIa. According to the criteria for statistical reliability, the 
AGCDM model has certain advantages over the ACDM model. A graviton-ghost condensate lays claim to being a 
variable component of Dark Energy. Result of observational data processing suggests that during the era of large- 
scale structure formation in the universe, a graviton-ghost condensate played a measurably significant part in shaping 
global cosmological dynamics. 

Further applications of our theory to the physics of Dark Energy will be to conduct numerical experiments with the 
BBGKY chain, taking into account the non relativistic matter. The results of these experiments may be able to explain 
the finer details of the Hubble diagram for supernovae SNIa, which will be identified with increasing completeness 
and accuracy of observation data. Future observations, in our view, must focus on the variable component of Dark 
Energy. For reliable identification of this component it is necessary, first move as far as possible into the area of large 
redshift, and second to explore in detail the area of small red shift, in which, quite possibly, there is a transition from 
the pre-asymptotic state of graviton-ghost condensate to its asymptotic state. 
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XI. APPENDIX I. COSMOLOGICAL CONSTANT AS ENERGY DENSITY OF EQUILIBRIUM 

VACUUM 



The cosmological constant in Einstein equations does not contradict the general principles of geometrized gravity 
theory. Moreover, the gravity theory which includes the A-term looks more natural than a theory without A-term. 
According to Zeldovitch [l3| , A-term is interpreted as the energy density of the vacuum. Today, that definition should 
be "updated" a little bit: the subject is the energy density of the equilibrium vacuum subsystems of non-gravitational 
origin in 4-D space-time: elac = A. The calculation of eiac is the task of the future "Theory of Everything" , 
based, possibly, on the superstring theory. The existing experimentally verified theory (the Standard Model of quark 
and leptons interactions) allows making some general conclusions. The most important of these is the fact that for 
every fundamental physical interaction there is an associated vacuum subsystem with non-zero energy density. They 
are well-known Higgs condensate in the theory of electro-week interaction and quark-gluon condensate in Quantum 
Chromodynamics . 

The Higgs condensate is forming as a result of spontaneous breaking of electro-week symmetry U{1) x SUl{2) 
down to the electromagnetic symmetry Uemi^)- The energy density of Higgs condensate can be estimated as (in that 
section h = 1): 

= - i^'^i + + 6A4 - UM^) , (XI.l) 

where AIh, Mz, A/„,, Mt are masses of Higgs boson, intermediate vector bosons and t-quark; — 47r/29 is the 
gauge constant of SUl{2) group on VF-boson mass-shell. The first term in (|XI.ip is the energy density of spatially 
homogeneous (vacuum) component of the scalar Higgs field (0|i?|0) — the second term is the change of the 

energy of zero-fluctuations of quantum fields, which have obtained nonzero rest mass from the interaction with the 
vacuum field v = 246 GeV. Masses of all particles M ^ v, therefore 

£ew(«) = -jA«4c.-(120GcV)\ (XI.2) 

where A is a constant or (if high radiation corrections are taken into account) is a very slow function of vacuum field. 
The numerical value used in (jXI.2p corresponds to the mass of Higgs boson AI„ = 2Af„. ~ 160 GeV. There is no 
doubt in the existence of the electro-week vacuum subsystem at the scale Xew ~ 120 GeV, however a fact that this 
subsystem is specifically formed by the Higgs mechanism, is not yet experimentally confirmed. A decisive step in that 
direction would be experimental discovery of the Higgs boson by the LHC. 

The quark-gluon condensate is a system of mutually correlated non-perturbative fluctuations, arising during 
quantum-topological tunnel transitions between degenerated states of giuon vacuum. Energy density of the con- 
densate is 

Sqcd = ^■577(01 — G'^j.Gjf |0) , 

32 (XI.3) 

6 ~ 9 + 8Tg{mu + mrf + 0.8m^) ~ 9.6 , 

where 

(Ol^G^.Gf |0) ^u^^ (360 MeV)4 (XI.4) 

TT 

is the main energy-momentum parameter of the quark-gluon condensate; Tg ~ (1.5 GeV)~^ is a characteristic space- 
time scale of fluctuations; rUu, rud, rris are the masses of light quarks. According to the modern paradigm, the 
quark-gluon condensate has several phase states, in each of which the fluctuations have there own specific microstruc- 
ture. The value of the parameter, used in (|XI.4p . refers to the confinement phase of the out-of-hadron vacuum. In 
that case 



b_ 
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,Qc„(u) = ---u4 ^ „(265 MeV)-* (XI.5) 



As follows from (jXI.2p and (jXI.Sp . the scale of the electro-week vacuum is of the order of Mew ~ 120 GeV, and 

QCO 



the QCD scale is Mq^d — 265 MeV. The scale M^e ^ 10 "^"^ GeV, corresponding to the current density of the Dark 



Energy e^E ~ Af^^, is not on a pair with Mf.w, Mqco- It means that for the high-energy expansions of the Standard 
Model, the vacuum is of a more complicated structure. First, there should be vacuum subsystems with positive energy 
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density; Second, the theory should contain mechanisms of mutual compensation of the contributions from different 
vacuum subsystems to the full energy density of the equilibrium vacuum. 

The Standard Model, when applied to the description of the cosmological vacuum and cosmological plasma of 
elementary particles, forecasts that the vacuum subsystems (IXI.ip . (jXI.3[) are of evolutionary origin. The energy 
density of these subsystems is not constant by the definition. Higgs' and quark-gluon condensates appear in the 
early Universe in the processes of relativistic phase transitions. The energy density of the vacuum is a functional of 
non-equilibrium parameters V ^ v^lA ^ u: 

e.ac[VM) = Ao + eev.(V) + e^coiU) , (XI.6) 

where Aq = const is the energy density of the vacuum at the evolutionary stage before the electro-week transition. 
In the areas of phase transitions, functions eew(V) and SQcoiM) are changing from zero values to equilibrium values 
given in (jXI.2p and (jXI.Sp . The characteristic transition times and the times of relaxations of vacuum subsystems 
to the equilibrium states are incommensurably shorter than the characteristic times of the Universe evolution. By 
that reason, almost immediately after the phase transition, one might talk about the contributions of the respective 
condensates of the Standard Model to the cosmological constant. After the quark-hadron transition, we have the 
current value of the non-gravitational contributions to A-term: 

A = Ao + £e^v {v) + £qcd (m) ■ (XI.7) 

Thus, the theory of fundamental interactions, particularly, the Standard Model, leads to the conclusion that the 
treatment of the cosmological constant as an energy density of the equilibrium vacuum, is related to the late stages of 
the Universe evolution, situated on the cosmological scale after all the relativistic phase transitions in the vacuum and 
in the plasma of elementary particles are completed. 

Inevitability of the "fine tuning" of different vacuum subsystems, with the result that the asymptotical value of 
the cosmological constant (|XL7p turns out to be very small or even equal to zero, is actually the experimental fact, 
resulting from the very existence of the modern Universe. Actually, without the "fine tuning" any of the phase 
transitions from the Standard Model will lead either to the collapse of the Universe, or to the exponential inflation, 
preventing the formation of large-scale structures. At the same time, strong inequality |eew(t^)| 3> \£qcd{u)\ forces 
us to consider separately the problem of the "fine tuning" at the electro-week and QCD scales. 

The compensation of the energy of the Higgs condensate (jXI.ip at the fundamental level can, in principle, be 
provided by the supersymmetry of particles and interactions. Unfortunately, in the superstring scenarios of the 
"Theory of everything" , the effective algorithms to develop that kind of theory are absent at the present time. The 
supersymmetry, if it actually exists in Nature, is strongly broken, but the mechanisms of destruction, it seems, are not 
spontaneous and 4-D spacetime is not an arena of an action of these mechanisms. It is not excluded, that the reduction 
of the " Theory of Everything" to the low-energy Standard Model is provided by physical phenomena working in the 
extra spatial dimensions. Experimental test of the supersymmetry, as a realistic concept of the elementary particle 
physics, will be conducted at LHC. The discussion about the compensation mechanisms for the energy density of 
Higgs condensate (|XL1I) should continue after the results of the experiments start arriving. 

Supersymmetry has no relation to the problem of compensation of the quark-gluon condensate (|XI.3p . We think 
the problem should be solved in the framework of low-energy physics of strong interactions. One of the possible 
approaches is to assume that at the QCD scale there are additional contributions to (jXI.3[) . and the sum turns to zero 
due to the tuning of QCD parameters. The discussion of these scenarios is out of the scope of the current work. 

The hypothesis about the mutual compensations of non-gravitational contributions to the A-term do not remove 
the questions regarding the final value of A ^ 0, which is in agreement with observed data. At the order-of-value 
estimations, A-term, with its value equal to the modern Dark Energy density, satisfies Zeldovitch's relationship [4l| 
with the Kardashev's modification [9^ 

~ A = = 3 X 10-4^ GeV^ , (XI.8) 

where m^r — 138 MeV — pion mass; Mpi ~ 1.22 • 10^^ GcV — Planck mass. As it is known, the pion mass is an object 
of non-perturbative QCD. Therefore, in formula (|XI.8p the observed density of Dark Energy is expressed only via 
the combination of minimal Aqcd — 2mTr and maximal scales Mpi of the elementary particle physics. Adopting also 
the Sakharov's idea of Einstein equations modification by the effect of gravitational exchange interaction of identity 
particles [oTj . we may suggest that the energy density of vacuum (|XI.8p is formed by the gravitational exchange 
interaction of quantum-topological fiuctuations in the hadron vacuum. This idea leads to a formula 

ttT^ 7^-1 
- A = ^\^<^on fin' -r^ - 4 X 10-4^ GeV^ . (XI.9) 
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The difference of four orders between values (|XI.8|) and (jXI.9|) is no any longer a catastrophe for the theory (see 
formula (jXI.lOp below). 

The nature of the cosmological constant discussed above, (more precisely, the contributions to the energy density 
of the vacuum controlled by the modern theory of elementary particles), suggests a set of interesting analogies. 
First, because the vacuum subsystems of electromagnetic, week and strong interactions are known from experiments, 
it is appropriate to ask a question: what is the vacuum subsystem connected with the fundamental gravitational 
interaction? We hope that an answer is provided in our work: the existence of the graviton-ghost condensate, i.e. 
the Dark Energy, is a direct consequence of applying the first principles of quantum theory of gravitation. The 
second analogy is in the fact that all vacuum subsystems are the subjects of evolution as Universe progresses, and the 
asymptotical states appear as results of relativistic phase transitions. From the evolutional criterion, the difference 
between graviton-ghost vacuum subsystem and non-gravitational vacuum subsystems is of a quantitative character 
only. Non-gravitational vacuum subsystems relax on the respective characteristic scales of elementary particle physics, 
while due to the weakness of gravitational interaction, the graviton-ghost condensate relaxation takes place at the 
time intervals which are of the order the age of the Universe. Finally, the third general property of all vacuum 
subsystems is that the vacuum energy density in the asymptotical (equilibrium) state acquires a constant value. It is 
worth mentioning that the graviton-ghost condensate in its equilibrium state, as well the vacuum of QCD, is of the 
instanton microstructure. 

Of course, the problem of mutual compensation of different contributions to the energy density of non-gravitational 
vacuum subsystems remains unsolved. However, we believe that the consideration of the graviton-ghost component 
of the physical vacuum somewhat reduces the magnitude of the problem. Actually, as shown at the end of the section 
IVIBl the mathematical structure of the self-consistent theory of gravitons itself possesses an effective renormalization 
procedure applied to the energy density of non-gravitational vacuum by virtue of the graviton-ghost contributions. 
The possibility of such renormalization is here regardless of the sign of the A-term. Particularly, the renormalization 
of A-term (IXI.Qp . formed by the gravitational exchange interactions of non-perturbative quark-gluon fluctuations, is 
conducted by formula (|VI.34p and leads to the result: 



As follows from (|XI.10|1 . ehac 10"*^ GeV* with \Ng\ ^ 10^^®. Unfortunately, the question regarding the concordance 
of the graviton-ghost condensate's parameter |iVg| with the numerical values of parameters of a different physical 
nature (example (jXI.lOp . with QCD parameters), remains open. But in that case one always may appeal to the 
anthropic principle. 



The problem of calculating the anomalies in the energy-momentum tensor of gravitons (quantum field with spin 
J — 2) was discussed in [l8, 92, 93]. A large number of works are devoted to the study of graviton quantum field 
in De Sitter space. Here we provide only some links. In some works [13, [2l| the non-ghost models were considered; 
in other works [2^. the ghost fields with harmonic gauge were taken into account. Common feature of all known 
versions of quantum theory of gravitons in the isotropic Universe is the lack of one-loop finiteness off the mass shell. 

In Sections IXII Al IXIIBI we discuss a self-consistent theory of gravitons in isotropic Universe with the ghost sector 
not taken into account. As has been repeatedly stated, we believe that such a model is not mathematically sound. 
Gauges, completely removing the degeneracy, are absent in the theory of gravity. Thus, in the self-consistent theory of 
gravitons the ghost sector is inevitable present. Now, however, let us assume for the moment that the self-consistent 
theory of gravitons without ghosts is worth at least as a model of mathematical physics. The purpose of this Section 
is to get the properties of this model and to show that it is mathematically and physically internally inconsistent. 
Note also that the self-consistent theory of gravitons in the harmonic gauge (jIII.79P has qualitatively the same status 
as the non-ghost model: the one-loop non-renormalizability is hindering the attempts to obtain definitive results 
within this theory. 

A. Gravitons with no Ghosts. Vacuum Einstein Equations with Quantum Logarithmic Corrections 

It clear from the outset that in the non-ghost model the calculation of observables will be accompanied by the 
emergence of divergences. It is therefore necessary to formulate the theory in such a way that the regularization and 
renormalization operations are to be contained in its mathematical structure from the very beginning. We talk here 
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about changes in the mathematical formulation of the theory. The relevant operations should be introduced into the 
theory with care: first, in the amended theory, coexistence of classical and quantum equations should be ensured 
automatically; second, the enhanced theory should not contain objects initially missing from the theory of gravity. 

The dimensional regularization satisfies both above-mentioned conditions. Important, however, is the following fact: 
the use of dimensional regularization suggests that the self-consistent theory of gravitons in the isotropic Universe 
is originally formulated in a spacetime of dimension D = I + d, where 1 is the dimension of time; d = 3 — 2e is 
the dimension of space. The special status of the time is due to the two factors: (i) all the events in the Universe, 
regardless of its actual dimension, are ordered along the one-dimensional temporal axis; (ii) the canonical quantization 
of the graviton field in terms of the commutation relations for generalized coordinates and generalized momenta also 
presuppose the existence of the one-dimensional time. As for the space dimension, the limit transition to the true 
dimension d = 3 is implemented after the regularization and renormalization. 

Thus, we are working in a space with a metric 
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To avoid mathematical contradictions that could arise at the limit d — s- 3, Einstein equations in Z)-dimensional 
spacetime should be written down in exactly the form in which they were obtained from the variational principle: 
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Here >Cd is the Einstein gravitational constant in ZJ-dimensional spacetime. (Dimension [xdfi] = [l]^ ^0 The left 
hand sides of equations (|XII.2p satisfy the Bianchi identity: 
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In the right hand side of equations (jXII.2[) . the identity (|XII.4[) generates condition of the graviton EMT conservation 
that satisfies if the equations of motion (|XII.3|) are taken into account. Regarding the origin of the system of equations 
()XII.2|) and (|XII.3p . we should make the following comment. In this case it is inappropriate to invoke the reference 
to the path integral and factorization of its measures because the path integral inevitably leads to the theory of 
ghosts interacting with the macroscopic gravity. We can only mention a heuristic recipe: one should refer to the 
density of Einstein equations with mixed indices, define the exponential parameterization of the metric, and expand 
the equations into a series of metric fluctuations with an accuracy of the second-order terms. Deviations from this 
recipe (for example, linear parameterization gn- = gik + hik) lead to a system of inconsistent classical and quantum 
equations. To remove this sort of inconsistency, one is forced to use artificial transactions outside the formalism of 
the theory (see, for example, [13). 

While working with the system of equations (|XII.2[) . (jXII.3[) . we face with two mathematical problems. The first 
problem is that in the framework of that system of equations, except in very special cases, it is impossible to formulate 
the dynamics of operators on a given background that is to get the solution of the equation (jXII.3|) as an accurate 
operator function of time. This is due to the fact that formulae of (|XII.2p in reality are not yet specific equations. 
They are only a layout of Einstein equations with radiation corrections. These equations can only be obtained 
after regularization and renormalizations of the ultraviolet divergences. In addition, the functional form of equations 
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depends on which quantum gravitational effects are to be taken into account outside the sector of vacuum (i.e. zero) 
fluctuations of the graviton field. The only possible way to study the system of equations (|XII.2p and (|XII.3p is (i) 
to obtain the solution of operator equation (jXII.3|l in a form of a functional of the scale factor without specifying 
the dependence on a(r]) with a clear emphasis on zero fluctuations in this functional, (ii) to substitute the obtained 
functional in (jXII.2p under certain assumptions about the state vector; (iii) to regularize and renormalize and finally 
(iv) to solve the macroscopic Einstein equations, obtained after these operations. Implementation of the program, 
an essential element of which is the allocation of zero fluctuations generating ultraviolet divergences, is possible only 
when using the method of asymptotic expansions of solutions of operator equation in the square of wavelength of the 
graviton modes. Thus, the problem of the lack of macroscopic Einstein equations in the original formulation of this 
theory with divergences limits the methods of this theory to the short-wave approach. Note that this fact was clearly 
indicated by DeWitt [l3|- 

The second problem is related to the infrared instability of the theory, with the object of the theory being a 
conformal non-invariant massless quantum field. The problem is due to the fact that not every representation of 
the asymptotic series can be substituted into energy-momentum tensor to perform the summation over the wave 
numbers. For example, if in the explicit form, a term in the asymptotic series contains a large parameter k^" in the 
denominator, then starting from n = 2 in the integration over the wave numbers the infrared divergences will appear. 
Such an asymptotic series can not be used even for the renormalization of ultraviolet divergences, because when it is 
used in the space of the physical dimension rf = 3, the logarithmic divergences arise simultaneously at the ultraviolet 
and the infrared limits. In the method of dimensional regularization the problem is reduced to the fact that it is 
impossible to choose an interim dimension d in a way such that the integral exists at both limits. 

Formally, the technical problem described above is partly solved by reformatting the asymptotic series. In particular, 
the following method will be used, in which parameter of the asymptotic expansion is the effective frequency 
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In this method, the integrals over the wave numbers can be defined in terms of the principal value. Contributions of 
the poles at fc = ^/—p can not be mathematically verified if only because there are such contributions from each term 
of the infinite asymptotic series. The inability to describe infrared effects is the principal disadvantage of a theory 
with divergences, which uses only asymptotic expansions with respect to the wavelength. Meanwhile, as general 
considerations and the results of this work show, in the physics of conformal non-invariant massless field the most 
interesting and innovative effects occur in the infrared spectrum. The method of describing these effects, based on 
the exact BBGKY chain, can not be used in the theory with divergences, because a method regularizing the infinite 
chain of moments of the spectral function does not exist. 

The above problems automatically reduces the interest toward the theory with divergences. However, given that 
all previous works in this area have been implemented in the framework of regularization and renormalization, let 
us conduct our analysis to the end. In calculations, it is enough to consider the equation for the convolution. After 
identity transformations, using the equation of motion (|XIL3p . we get 
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where 



is the spectral function of gravitons. The calculation of the spectral function by the method of asymptotic expansion 
with respect to the square of wavelength was described in Section IIV Al Now we need to repeat this calculation 
excluding the ghosts, but with input from zero fluctuations in the spacetime of dimension D = d+1. The relevant 
calculations do not require additional comments. A spectral function is represented as: 
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where W, 



(vac) 



kcr vacuum component of the spectral function and W^f^^ is the spectral function of excitations. 

After passage to the limit d ^ 3, the contribution of to the EMT of short gravitons is exactly the same as 

pV.gp . (jlV.lOp . In the future, we discuss only the contribution from vacuum components of the spectral function. In 
the calculations, we must keep in mind that in the d-dimensional space the number of internal degrees of freedom of 
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transverse gravitons is Wg = [d + l){d— 2)/2. The solution for the vacuum spectral function is expressed in terms of 
the functional (|IV.4[) : 
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The powers of operator j| • 1 are defined by formulas (jIV.Sp . in which has the form (|XII.5[) . After substitution 
of (|XII.8[) into (|XII.6[) . the zero-term in the asymptotic expansion creates an integral, calculated by the rules of 
dimensional regularization: 
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The F— function in (|XII.9p diverges for d ^ i. Therefore, calculation of the integral (|XII.9[) and transformation of 
expressions with F— functions are carried out with those values of d which provide the existence of the integral and 
F— functions. At the final stage, the result of these calculations is analytically continued to the vicinity d = 2>. All other 
terms of the asymptotic expansion (jXII.Sp generate finite integrals and do not require a dimensional regularization. 
For reasons of heuristic rather than mathematical nature, it is considered that these terms are negligible compared 
to the contribution of the principal term of the asymptotic expansion (see below the effective Lagrangian (jXII.19|) ). 
Convolution of Z)-dimensional Einstein's equations (|XII.6p . containing the main term of the vacuum EMT of gravitons, 
has the form: 
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Other Einstein equations can be obtained using the Bianchi identities. A complete system of Einstein vacuum 
equations is written in Z?-covariant form: 
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Equation (jXII.ll|) are obtained by the variation of action 
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It is obvious from (|XII.lip . (|XII.12p that the method of dimensional regularization retains overall covariance of the 
theory. Of course, quantum corrections, appearing in (jXII.llj) . satisfy the condition of conservation. 

Renormalization and removal of regularization (limit d — > 3) are held at the level of action. A parameter with the 
dimension of length, which will eventually acquire the status of renormalization scale, is contained within the theory. 
This parameter, referred to as Lg, is appears in the Z3-dimensional constant of gravity: 
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The technique of removal the regularization assumes conservation of dimensionality for those objects in which the limit 
operation is performed. There are two such objects: the measure of integration d/i and the density of the Lagrangian 
£. As can be seen from (|XII.12p . (|XII.13p . the first (Einstein) term of the action is written down as 
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where D-dimensional objects C and have the same dimensions as the corresponding 4-dimensional objects. In this 
sector of the theory the Umit transition is trivial: R{d) ^ R, dfj, ^ ^y—gd'^x. In the sector of quantum corrections to 
the Einstein theory, we introduce the same measure and obtain the density of the Lagrangian: 
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It is necessary to emphasize that the operations of renormalizations and removal of regularization have to he mathemat- 
ically well-defined and generally-covariant. The condition of mathematical certainty assumes that the renormalization 
is conducted before the lifting of regularization. At the same time, the general-covariance of the procedure is auto- 
matically fulfilled if the counter-terms imposed in the Lagrangian are the D-dimensional invariants. Note also that 
if the mathematical value is finite at d = 3, then the above formulated conditions do not prevent the expansion of 
this quantity in a Taylor series over the parameter (3 — d)/2. In particular, we can write: 
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where /ig — 1/Lg\ ellipsis designate the terms which do not contribute to the final result. The substitution (|XII.16[) 
in (jXII.lSp provides: 
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According to (|XII.17p . the source Lagrangian of the theory requires a ZJ-invariant counter-term, which removes the 
contribution proportional to the diverging F-function: 
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In (|XII.18p . there is a new finite constant of the theory of gravity 1//^. The removal of the regularization in the 
renormalized Lagrangian is conducted by the regular transition: 
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where 
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is the renorm-invariant scale. There is a heuristic argument allowing to use the obtained expression: quantum 
corrections in the Lagrangian ( rcf (12.19)) dominate over all other neglected terms of the asymptotic series over The 
logarithmic parameter ln{X^/R) ^ 1. 

The renormalized Einstein vacuum equations with quantum corrections obtained from the Lagrangian (jXII.19|l are 
as follows: 
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Note that exactly the same equations are obtained from Z)-dimensional equations (jXILlip , provided that the opera- 
tions are performed in the same sequence: first a renormalization with the introduction of D-covariant counter-terms 
is conducted, and then a limit transition to the physical dimension is performed. 
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B. Intrinsic Contradiction of Theory with no Ghosts: Impossibility of One-Loop Renormahzation 



We are still discussing a formal model — self-consistent theory of gravitons with no ghosts. In the previous section 
it was shown that the renormahzation of divergences, that inevitably arise in this model, requires the imposition 
of an additional term quadratic in the curvature in the Lagrangian. It is now necessary to draw attention to two 
mathematical facts: (i) the need for a modification of Einstein theory is caused by quantum effects contained in the 
Lifshitz operator equation (jXII.3p : (ii) the original Lagrangian and operator equations of the modified theory have 
the form: 
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where is a covariant derivative in a space with the operator metric (jik- It is quite obvious that these facts 
contradict each other: the quantum effects in the Lifshitz equation lead to a theoretical model that contradicts the 
Lifshitz equation. Let us demonstrate that the contradiction is a direct consequence of the non-renormalizability of 
the model (|XII.2ip off the graviton mass shell. 

Equations (jXII.22p . after their linearization describe quantized waves of two types — tensor and scalar. It makes 
sense to discuss the problem of the scalar modes only in the event that at least preliminary criteria for consistency of 
modified theory will be obtained. Therefore, first of all, we should reveal properties of the tensor modes. Here is an 
expression for the Lagrangian of a system consisting of self-consistent cosmological field and tensor gravitons: 
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The equation for gravitons is produced either by the linearization of the equation (|XII.22[) . or from (|XII.23[) by the 
variation procedure: 
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Please note that the last term in (|XII.24p makes it impossible to retain the Lifshitz equation. After the transformation 
equation (jXII.24[) has a form 
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In (|XII.25p , the deviation from the Lifshitz equation is manifested in the effective frequency of gravitons — the latter 
contains an additional function of curvature's derivatives 
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When calculating quantum corrections to the macroscopic equations, the modification of the effective frequency 
leads to additional divergences. Averaged vacuum equations (jXII.22p . after their polynomial expansion in powers of 
curvature, look as follows (finite logarithmic corrections are omitted): 
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Here r(e) ~ 1/e is a divergent F-function obtained by dimensional regularization; I//0 is a seed constant of a theory 
with quadratic invariant. The complete quantum Lagrangian corresponding to equations (jXII.27[) has the form: 

-/ 

Renormalization of the second term in (IXII.28P is performed by selecting the seed constant: 

/2 2887r2 p ■ 

However, a divergent coefficient forms before the third term. To overcome this divergence, it is necessary to introduce 
a new seeding "fundamental" constant of the modified theory of gravity l/h^ with a renormalization rule: 
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Further actions are obvious and pointless: Lifshitz equation is the subject of the next modification; quantum 
corrections generate another new divergence; to renormalize the new divergence a new theory of gravity is introduced, 
etc. The only conclusion to be drawn from this procedure is that based on the criteria of quantum field theory, 
the one-loop self-consistent theory of gravitons in the isotropic Universe, and not possessing the property of one- 
loop finiteness outside of mass shell, does not exists as a mathematical model. In such a theory it is impossible to 
quantitatively analyze any physical effect. The theory of gravitons without ghosts is non-renormalizable even in the 
one-loop approximation. It is also important to stress that the correct alternative to a non-renormalizable theory is 
only a finite theory with the graviton-ghost compensation of divergences. 

In the future, from our perspective, the method of regularization and renormalization in general will be excluded 
from the arsenal of quantum theory of gravity, including one from the theory of one-loop quantum effects involving 
matter fields. Correct alternatives to existing methods of analysis of these effects to be found in extended supergrav- 
ities, finite at least in one-loop approximation. 

The situation prevailing in the scientific literature is a paradoxical one. On the one hand, inadequate nature of 
the regularization and renormalization methods in the quantum theory of gravity should be obvious from the latest 
development trends in the theories of supergravity and superstrings. On the other hand, however, in all works we 
know on cosmological applications of one-loop quantum gravity theoretical models are used, which, according to the 
criteria of quantum field theory, do not exist. We cannot comment on the specific results obtained in these models by 
the reasons clear from the content of this Section. Once again we should emphasize that the self-consistent theory of 
gravitons, if it exists as a theoretical model, must be finite outside the mass shell of gravitons. Effects arising in the 
finite theory are described in the main text of this work. 

C. Dimensional Transmutation of Finite Theory 

One-loop finiteness is the central important feature of the quantum theory of gravity, defining its mathematical 
structure and the algorithms of concrete calculations. Nevertheless, one should bear in mind that in the calculations 
we are dealing with, compensation of divergences comes from the graviton and ghost sectors. Therefore, the question 
arises: how are the structure and the results of the theory changed when an intermediate regularization of divergences 
is applied? Comparison of results obtained by different methods of computation will allow us to judge the stability 
or volatility of the results of one-loop quantum gravity with respect to the intermediate regularization. 

We will continue to use the method of dimensional regularization. The passage to the limit D — 4 will be applied 
only in the final expressions. General considerations show that after the transition to Z? = 4 in the expressions for 
the observed values, some terms may appear, whose mathematical structure has no analogues in the equations of the 
original theory. In such cases we are forced to talk about quantum anomalies that have arisen as a result of dimensional 
transmutation. If there are no such terms, the theory is stable with respect to the dimensional transmutation, and 
has no anomalies. Here we show that one-loop quantum gravity (without the matter fields) not only finite, but is 
void of anomalies. 

Let us turn to the convolutions of Einstein equations, assuming that the latter are written exactly in the form in 
which they were obtained from the path integral: 
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In the study of dimensional transmutation, one can not divide the left and right hand sides of equation (|XII.29p by 
the common multiplier y/\g(d)\, because the limit d ^ 3 on the left hand side is a regular one, but on the right hand 
side the same limit applies under conditions of compensation of divergences. In carrying out the operations in the 
right hand side, one must take into account the total dependence on the parameter d. 

The metric of rf-dimensional isotropic Universe and its scalar curvature are represented in the form of (jXII.ip . 
A theory of gravitons formulated in that spacetime has undergone a preliminary investigation. In doing so, it was 
established that (i) the ghost sector still consists of one complex Grassmann field, satisfying the Klein-Gordon-Fock 
equation in the space with a metric (jXII.ip . Thus, the number of internal degrees of freedom of the ghost field is 
Wgh — 2; (ii) the number of d— tensor gauge invariant degrees of freedom (the number of transverse polarizations) of 
the graviton field does not match the number of internal degrees of freedom of ghosts: 
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(iii) the equations for the wave functions of gravitons and ghosts have the same form: 
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(iv) the trace of the energy-momentum tensor of gravitons and ghosts in (|XII.29[) is as follows: 
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Each of two terms on the right hand side (|XII.32p contains the contribution of zero fluctuations of quantum fields, 
which, because (IXIL31I) . give rise to divergences of the same type. With d = 3 these divergences exactly cancel out 
each other out, because Wg = Wgh =2 — this comes from the nature of one-loop finiteness of the quantum theory of 
gravitons in the isotropic Universe. However, at d = 3 — 2£ there is no exact compensation: as seen from (IXII.30I) . 
Wg — Wgh = (d + 2)(d — 3)/2 ~ 5e. The difference in the numbers of internal degrees of freedom is multiplied by a 
divergent coefficient proportional to l/e. As a result, in the limit of d — *■ 3, a conformal anomaly arises, caused by 
the spontaneous dimensional transmutation in the theory of gravitons and ghosts. 

The goal is to calculate the conformal anomaly as a functional of the spacetime metric. Clearly, this requires 
solutions of operator equations (|XII.32|) in the form functionals of the same-type. As we know, when examining 
the zero fluctuations by the methods of regularization and renormalization, it is enough to have the solutions in the 
form of asymptotic expansion in powers of curvature. Virtually all the computations coincide with the computations 
already described in Section IXII Al down to the formula (jXILlOp — with the only difference being that it is now 
necessary to take into account the additive contributions of gravitons and ghosts. The trace of the energy-momentum 
tensor (|XII.32[) is divided into two terms: 
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The first of these contributors describes zero vacuum fluctuations, "deformed" by the self-consistent gravitational 
field, with the graviton-ghost compensation (which is incomplete with d 7^ 3) taken into account: 
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Note that the relation (|XII.30p was used to obtain (|XII.34[) . The second term in (jXII.33p is a trace of energy- 
momentum tensor of the excitations. Because of the obvious limitations on the range of excitations, the term has no 
divergent integrals, so the limit is taken in the regular manner. In doing so, of course, an expression is obtained that 
exactly coincides with the trace of the energy-momentum tensor in 4-dimensional finite theory not containing the 
contribution of zero fluctuations: 



:^(vI;|T(-^)|vI/) = _±^Z?, 

OK 



(XII.35) 



where D is a function appearing in the 4-dimensional finite BBGKY chain (|V.7|) — (|V.9[) . 
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Note that previous calculations in this section and further operations with the expression (|XII.34p are formally 
accurate. As for the (|IV.5[) . we should notice that all integrals of expression 
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have regular limit at d ^ 3, so the multiplication of these integrals by the factor d—3, would provide zero. The effect 
of dimensional transmutation is contained entirely in the expression 



-gT 



(vac) 



lim — (n- l)(n + 2)(n-3) 

n-*3 16 



(XII. 36) 



The integral in dimensional space of wave numbers is defined in (|XII.9[) . Substituting this expression into (jXII.36p . 
we get the final result: 
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From (|XII.37p . using the Bianchi identities, one can reconstruct all components of the anomalous energy-momentum 
tensor of gravitons and ghosts: 
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Expression (jXII.38|) describes a concrete quantum effects, which is the deformation of the spectrum of zero fluctuations 
of gravitons and ghosts in the self-consistent classical gravitational field. In obtaining (|XII.38p . the renormalization 
with the removal of divergences by introducing the counter-terms was not used. The existence of this effect does not 
require modification of the original quantum Lagrangian. 

In an isotropic space with the metric (jXII.ip , two types of conformal anomalies of the energy-momentum tensor of 
quantum fields were repeatedly discussed: 
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(Numerical coefficients C^P and cf^ for the fields for the fields with spin J = 0, 1/2, 1 are given, for example, in the 
monograph [6].) The anomaly of second type has exited an interest, in particular because the quantum corrections 
(jXIL40p , added to the Einstein equation, are capable of providing a self-consistent De Sitter solution in the vicinity 
of the Plank's values of the curvature [93|. As we can see from (|XII.38p . the dimensional transmutation of the finite 
theory of gravitons and ghosts generates only the first type anomaly in the energy-momentum tensor. The anomaly of 
the second type simply does not arise under the dimensional transmutation. It means that in the one-loop quantum 
gravity (without matter fields) the effects of zero fiuctuations are not able to sustain an infiationary expansion with the 
constant parameter of inflation. In the finite theory, the De Sitter solution can be formed only by the graviton- ghost 
condensate. 

The last question is this: does the emergence of anomalies (IXIL38I) change of the mathematical structure of the 
system of equations of the theory consisting of the BBGKY chain (|V.7p — (|V.9p and the macroscopic Einstein's 
equations (jV.lOp ? The answer is negative. The fact is that the BBGKY chain (|V.7p — (|V.9p is form-invariant with 
respect to the additive transformation of the moments of the spectral function: 

Wn^W^ + ^A_1_K-+' . 1 , n = 0, 1, 2, oo , (XII.41) 



where b is an arbitrary numerical parameter; Kp is an integral-differential operator, functionally dependent on p = 
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—a" /a. The operator is defined as follows: 

k;+' • 1 = • (i^; • i) . 

The transformation (|XII.41|) can be seen as a trace of renorm-group symmetry of a theory with the quadratic invariant 
— the divergences have disappeared, but their existence in the graviton and ghost sectors was separately recorded 
in the symmetry properties of the BBGKY chain. This transformation leaves the chain (jV.7[) — (jV.9[) unchanged, 
but the same transformation D and Wi with a coefficient 6 = — 5 eliminates the conformal anomaly from Einstein 
equations (jV.lOp . Thus, a simple renaming of the moments of the spectral function returns the system of equations 
to its former view. The result means that the one-loop quantum gravity not only finite, but anomaly-free as well. 

The anomaly (jXII.38[) is contained within the BBGKY chain along with all other quantum gravitational effects. In 
the general solution, the relative role of the anomaly is exclusively governed by the initial conditions. Of course, the 
particular solution of the chain (i.e. a solution containing only the anomaly) also has a meaning. In that particular 
solution, the value 6 = 5, obtained by the method of dimensional transmutation, is not a special one. In examining 
particular solutions of this type in the vacuum energy-momentum tensor (jXII.38p , the coefficient 5 should be replaced 
by an arbitrary constant b^. Different values of the constant bp physically correspond to different versions of incomplete 
compensation of energy of zero-fluctuations of gravitons and ghosts. The solution of such vacuum equations at 6o < 
("physics of the scalarons" ) was discussed in [9J]. 
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